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FOREWORD 


The  U.  S.  Environmental  Protection  Agency  is  charged  by  Congress  with  pro¬ 
tecting  the  Nation's  land,  air,  and  water  resources.  Under  a  mandate  of  national 
environmental  laws,  the  Agency  strives  to  formulate  and  implement  actions  lead¬ 
ing  to  a  compatible  balance  between  human  activities  and  the  ability  of  natural 
systems  to  support  and  nurture  life.  To  meet  this  mandate,  EPA's  research 
program  is  providing  data  and  technical  support  for  solving  environmental  pro¬ 
blems  today  and  building  a  science  knowledge  base  necessary  to  manage  our  eco¬ 
logical  resources  wisely,  understand  how  pollutants  affect  our  health,  and  pre¬ 
vent  or  reduce  environmental  risks  in  the  future. 

The  National  Risk  Management  Research  Laboratory  is  the  Agency's  center  for 
investigation  of  technological  and  management  approaches  for  reducing  risks 
from  threats  to  human  health  and  the  environment.  The  focus  of  the  Laboratory's 
research  program  is  on  methods  for  the  prevention  and  control  of  pollution  to  air , 
land,  water,  and  subsurface  resources;  protection  of  water  quality  in  public  water 
systems;  remediation  of  contaminated  sites  and  groundwater;  and  prevention  and 
control  of  indoor  air  pollution.  The  goal  of  this  research  effort  is  to  catalyze 
development  and  implementation  of  innovative,  cost-effective  environmental 
technologies;  develop  scientific  and  engineering  information  needed  by  EPA  to 
support  regulatory  and  policy  decisions;  and  provide  technical  support  and  infor¬ 
mation  transfer  to  ensure  effective  implementation  of  environmental  regulations 
and  strategies. 

This  publication  has  been  produced  as  part  of  the  Laboratory's  strategic  long¬ 
term  research  plan.  It  is  published  and  made  available  by  EPA's  Office  of  Re¬ 
search  and  Development  to  assist  the  user  community  and  to  link  researchers 
with  their  clients. 


E.  Timothy  Oppelt,  Director 

National  Risk  Management  Research  Laboratory 
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ABSTRACT 


CFC-1 14  has  been  used  as  a  working  refrigerant  in  shipboard  and  submarine  chiller  units.  Because 
CFC-114  is  one  of  the  designated  phase-out  refrigerants,  alternatives  for  CFC-1 14  in  shipboard  and  submarine 
chiller  units  are  being  sought  by  the  U.S.  Navy  with  the  assistance  of  the  Environmental  Protection  Agency  (EPA). 
Presently,  one  of  the  leading  candidates  to  replace  CFC-114  isHFC-236ea  (hexafluoropropane),  which  was  identified 
by  EPA.  This  study,  one  part  of  the  EP A/Navy  program,  is  an  evaluation  of  transport  properties  of  HFC-236ea. 

Liquid  viscosity  and  thermal  conductivity  of  HFC-236ea  are  the  two  main  transport  properties  of  interest  in  this  study. 

This  study  used  a  novel  method  for  simultaneously  measuring  viscosity  and  thermal  conductivity  by  using 
inline  property  sensors  in  series  with  a  heat  transfer  measurement  system.  For  this  purpose,  a  test  system  with 
inline  viscosity  measurements  and  a  single-phase  heat  transfer  test  facility  was  established  in  the  Refrigeration 
Laboratory'  at  Iowa  State  University.  In  this  approach,  viscosity  was  measured  with  an  inline  torsional  oscillation 
viscometer,  while  the  thermal  conductivity  was  determined  from  the  knowledge  of  the  single-phase  heat  transfer 
characteristics  of  a  heated  test  section. 

The  resulting  viscosity  and  thermal  conductivity  of  CFC-1 14  by  the  current  measuring  method  were 
compared  with  ASHRAE  (American  Society  of  Heating,  Refrigerating  and  Air-Conditioning  Engineers)  data  with  a 
difference  within  ±5%  in  thermal  conductivity  and  ±2%  in  viscosity.  For  HFC-236ea,  the  measured  data  were 
compared  with  predictions  of  REFPROP,  a  theoretical  package  developed  by  the  NIST  (National  Institute  of 
Standards  and  Technology),  with  an  average  deviation  within  +15%  in  thermal  conductivity  and  -5%  in  viscosity. 

The  properties  of  HFC-236ea  mixed  with  a  lubricant  (Castrol  oil  SW68)  were  also  investigated.  The 
results  showed  that  the  thermal  conductivity  increased  with  lubricant  concentration  in  a  low  temperature  range,  but 
decreased  with  the  mixture  in  a  high  temperature  range.  However,  there  was  no  significant  difference  of  thermal 
conductivity  of  refrigerant/lubricant  mixtures  for  mixtures  of  various  lubricant  concentrations  at  high  temperatures. 

For  the  viscosity  of  HFC-236ea/lubricant  mixtures,  the  results  showed  a  significant  increase  at  the  higher  lubricant 
concentrations  and  lower  temperature  ranges.  In  addition,  the  specific  heat  and  density  of  refrigerant/lubricant 
mixtures  were  also  determined  in  this  study. 

This  report  was  submitted  in  partial  fulfillment  of  CR  820755-01-4  by  Iowa  State  University  under  the  sponsorship 

of  the  U.S.  Environmental  Protection  Agency  with  funding  from  the  Department  of  Defense's  Strategic  Environmental 
Research  and  Development  Program  (SERDP).  This  report  covers  a  period  from  October  1992  to  May  1995. 
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M  —  molecular  weight,  g/mole 

m  -  mass  flow  rate,  kg/s 

m  -  mass  of  molecule,  g 

T/ud  —  Nusselt  number  based  on  diameter  D 

n  —  refractive  index  or  number  of  molecules  in  a  unit  volume 

P  —  pressure.  Pa 

Pc  -  critical  pressure.  Pa 

A P  —  differential  pressure,  Pa 

P^p  —  vapor  pressure,  Pa 

A Pr  —  differential  reduced  pressure,  (P  -  P^,)/Pc 

Pr  —  Prandtl  number 

Pr  -  reduced  pressure,  P/Pe 

Q,q  —  heat  rate,  W 

QvQr  -  total,  radiation  heat  rate,  W 

Qgon,  —  heat  transfer  by  convection,  W 

Qcond  —  heat  transfer  by  conduction,  W 

(f  —  heat  flux,  W/m1 

q£ss  -  loss  heat  flux,  W/m7 

q£t  —  net  heat  flux,  W/m2 

qt"ot  -  total  heat  flux,  W/m1 

qhss  -  heat  loss  rate  of  test  section,  W 

qnet  —  net  heat  input  rate  to  test  section,  W 

qlol  -  total  heat  loss  rate  of  test  section,  W 

R  —  radiation  parameter,  r  •  ka.  Equation  4.19 

Ra  —  Rayleigh  number,  Gr  •  Pr 

r  —radius  of  line  source,  m 

r,  -  radial  distance  of  point  i  from  a  line  source,  i  *=  0, 1, 2, ... 

ReD  —  Reynolds  number  based  on  diameter  D 

r0  —  wire  radius,  \un 

r,  —  inside  radius  of  tube,  m 

r2  —  radius  at  the  point  where  thermocouple  heads  are  buried. 


Abbreviations  (continued) 

SS  —  sum  of  square  quantity 


Tr,l 

Tra 
t  or  /, 

Urdk 


—  ambient  temperature,  °C 

—  boiling  temperature,  °C 

—  reduced  temperature 

—  static  temperature  or  outer  insulation  temperature,  °C 

—  total  temperature,  °C 

—  subsection  wall  temperature,  °C 

—  equilibrium  temperature  of  the  test  fluid,  °C 

—  average  inlet  temperature  of  the  test  section,  °C 

—  average  outlet  temperature  of  the  test  section,  °C 

—  average  fluid  temperature  of  the  test  section,  °C 

—  average  outer  insulation  surface  temperature,  °C 

—  average  wall  temperature  of  the  test  section,  °C 

—  average  ambient  temperature,  °C 

—  average  temperature  difference  between ,  Tw  and  T}t  °C 

—  ideal  temperature  rise  in  transient  hot-wire  method,  °C 

—  the  i-th  subsection  local  average  temperature  difference,  Tw,  -  °C 

—  dimensionless  temperature  group  in  test-tube,  (  Ta  -  T{ )  /  (  Tw-T^) 

—  temperature  at  a  radial  distance  i  from  a  line  source,  °G 

—  temperature  at  r1?  °C 

—  temperature  at  r2,  °C 

—  time  or  time  i,  sec 

—  uncertainty  percentage  of  k  by  the  Dittus-Boelter  correlation 

—  uncertainty  percentage  of  k  by  the  Petukhov  and  Popov  correlation 

—  uncertainty  percentage  of  k  by  the  Gnielinski  correlation 

—  velocity,  m/s 


Vs  —  speed  of  sound,  m/s 

Ax  —  distance  between  two  parallel  plates,  m 


Superscripts 

e  —  experimental 

c  —  correlation 

*  -  nondimensional  quantity 


Subscripts 


a 

app 

D 

f 

i 

net 

o 

oo 

0 

r 

s 

sat 

t  or  tot 
w 


—  air  or  absorption 

—  apparent  quantity 

—  diameter  based 

—  fluid 

—  inside  or  inlet 

—  net  quantity 

—  outside  or  outlet 

—  uniform  stream 

—  equilibrium  quantity 

—  radiation  or  reduced  property 

—  surface 

—  saturation  state 

—  total  quantity 

—  wall 


Greek  Symbols 


a 

£ 


Eu 

n 

p 

<*s 

© 


—  volumetric  thermal  expansion  coefficient  of  fluid,  K~ or  thermal  diffiisivity,  rrP-ls 
~  emissivity  of  surface 

—  emissivity  of  surface  /,  /  =  1,2, ...  etc. 

—  emissivity  of  surface  1  to  2 

—  viscosity,  N  •  slni2  or  cp  (centipoises  =  10-3  N  *  s/m 2 ) 

—  density,  kg/m3 

—  Stefan-Boltzmann  constant,  5.6697  x  10~8,  W/m2KA 

—  nondimensional  apparent  thermal  conductivity  group.  Equation  4.18 

—  acentric  factor 
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CHAPTER  1 


INTRODUCTION 


Chlorofluorocarbon  refrigerants  (CFCs)  have  been  widely  used  as  working  fluids  in  refrigeration  and  air- 
conditioning  applications.  In  the  past  few  decades,  CFC  refrigerants  have  been  used  as  refrigerants  due  to  their 
excellent  thermal  performance,  stable  chemical  properties,  and  low  toxicity.  However,  because  of  the  problems 
with  CFCs  that  were  exposed  in  the  early  1970's,  CFC  refrigerant  production  ends  on  December  3 1,  1995  under 
the  Montreal  protocol  and  its  amendments.  In  addition,  through  usage  regulation  and  prohibition,  the  CFC 
refrigerants  are  to  be  totally  phased-out  in  the  next  few  years  [1,2].  Since  the  CFC  problems  surfaced,  people  have 
sought  CFC  alternatives  to  take  the  place  of  existing  refrigerants.  So  far,  several  alternatives  have  been  found  for 
some  banned  refrigerants,  such  as  HFC-134a  for  CFC-12  and  HCFC-123  for  CFC-1I. 

A  cooperative  project  has  been  initiated  by  the  United  States  Navy  and  the  Environmental  Protection 
Agency  (EPA).  The  U.S.  Navy  currently  uses  CFC-114  ( CCIF2CCIF2)  as  the  working  refrigerant  in  shipboard 
chiller  units.  Because  CFC-1 14  performs  well  in  shipboard  chillers,  which  use  sea  water  as  the  heat  sink  fluid  on 
the  condenser  side,  CFC-1 14  has  served  as  a  good  working  fluid  in  the  past  However,  CFC-1 14  is  one  of  the 
banned  refrigerants.  With  the  mandatoiy  phase-out  program  of  CFCs  as  stipulated  by  the  Montreal  Protocol,  it  is 
an  urgent  task  for  the  U.S.  Navy  under  national  policy  to  find  a  new  refrigerant  alternative  which  is  non-ozone 
depleting  and  has  acceptable  thermal  characteristics  for  replacing  CFC-1 14.  Some  points  of  concern  for  the  Navy 
in  switching  refrigerant  to  an  acceptable  alternative  are:  1.  safety  (non-flammable,  non-toxic)  in  application;  and  2. 
easy  to  exchange  with  the  present  refrigerant  (CFC-114)  in  the  field  with  only  minor  equipment  modifications. 

HFC-236ea  (CF3CFHCF^f),  a  new  refrigerant,  has  already  been  selected  as  one  of  the  candidates  for 

CFC-114  alternatives.  There  are  several  reasons  that  this  refrigerant  shows  promise  to  replace  CFC-114.  Firstly, 
the  thermophysical  properties  were  investigated  and  found  to  be  similar  to  those  of  CFC-1 14  [3,4].  Secondly, 
according  to  simulation  results,  the  performance  of  HFC-236ea  is  similar  in  the  performance  to  CFC-1 14.  Thus,  it 
is  a  retrofit  candidate.  This  latter  characteristic  is  attractive  for  the  shipboard  or  submarine  chiller's  refrigerant 
transition.  Other  preliminary  test  results,  such  as  flammability,  material  compatibility,  and  toxicity,  are  also 
favorable.  The  EPA,  which  identified  HFC-236  isomers  as  promising  non-ozone  depleting  refrigerants,  has  been 
cooperating  with  the  Navy  in  evaluating  HFC-236ea  (and  HFC-236fa  as  well)  as  an  alternative  to  CFC-114  for 
shipboard  chillers.  The  EPA  effort  has  included  property  measurements,  performance  modeling,  and  small-scale 
engineering  tests. 
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PROPERTIES  OF  INTEREST 


Some  of  the  general  steps  for  refrigerant  replacements  include  property  evaluations,  heat  transfer  characterization 
studies,  design  evaluation  for  retrofit  application,  and  field  tests.  The  property  evaluations  of  interest  include 
thermodynamic  properties,  solubility  and  miscibility  characteristics,  and  transport  properties.  Some  of  the  properties, 
such  as  viscosity,  thermal  conductivity,  and  specific  heat,  are  important  properties  in  any  heat  transfer  study. 

Today  these  properties  are  known  for  numerous  fluids;  however,  there  are  still  potentially  useful  fluids 
whose  properties  need  to  be  measured.  Among  these  fluids  are  new  refrigerants,  refrigerant  mixtures,  and  the 
refrigerant/lubricant  mixtures.  This  study  focuses  on  the  transport  property  evaluations  of  the  alternative 
refrigerant  HFC-236ea.  The  transport  properties  of  interest  in  this  study  include  thermal  conductivity  and  viscosity 
of  pure  HFC-236ea  and  mixtures  of  HFC-236ea  with  lubricant. 


LUBRICANTS  OF  INTEREST 

Lubricants  are  required  for  the  proper  operation  of  the  mechanical  compressor  in  refrigeration  systems. 
However,  the  use  of  lubricants  will  cause  changes  in  a  refrigerant's  properties,  such  as  viscosity,  miscibility,  and 
solubility,  and  eventually  affect  heat  transfer  characteristics  [5].  Finally,  mixing  lubricants  with  refrigerants  could 
affect  material  compatibilities  [6]. 

The  mineral,  naphthenic,  paraffinic,  and  synthetic  lubricants  usually  used  with  CFC  refrigerants  are  not 
miscible  in  the  HFC  refrigerants  [7,8] .  Several  lubricant  alternatives  are  presently  considered  as  suitable  with  HFCs. 
These  lubricants,  which  include  polyol  ester  lubricants  (POEs),  alkylbenzenes,  and  polyalkylene  glycols  (PAGs), 
are  currently  being  tested  by  compressor  manufacturers  and  lubricant  suppliers  for  a  wide  range  of  characteristics, 
such  as  miscibility,  lubricity,  and  material  compatibility.  Recent  research  shows  that  the  POEs  and  PAGs  are  more 
favorable,  due  to  good  miscibility  and  stability  when  subjected  to  a  wide  range  of  operating  conditions,  for  new 
HFC  refrigerants  such  as  HFC-236ea  [9].  According  to  the  solubility  and  miscibility  test  results,  Icematic  Castrol 
SW68  lubricant  was  shown  to  be  compatible  with  HFC-236ea  [9]  and  was  selected  for  property  measurements  in 
this  study.  Various  concentrations  of  lubricant  which  cover  a  possible  range  of  concentrations  in  a  vapor 
compression  refrigeration  system  were  mixed  with  HFC-236ea.  Measurements  were  performed  on  the  refrigerant/ 
lubricant  mixtures  in  order  to  understand  the  effect  of  lubricant  concentration  on  selected  properties. 

OBJECTIVES 

Refrigerant  properties,  such  as  thermodynamic  properties  and  transport  properties,  provide  important 
information  as  to  the  acceptability  of  a  refrigerant  for  use  as  a  working  fluid  in  a  refrigeration  cycle.  Thermal 
conductivity  and  viscosity  are  especially  important  in  any  heat  transfer  analysis  as  applied  to  refrigeration  cycle 
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performance  studies.  Because  HFC-236ea  is  a  new  refrigerant,  there  are  very  limited  property  data  published  for 
this  potential  refrigerant.  As  a  result,  the  objectives  of  this  research  focus  on  measuring  the  transport  properties  for 
HFC-236ea.  The  transport  properties  of  interest  are  liquid  thermal  conductivity  and  viscosity.  The  lubricant's 
effects  on  the  transport  properties  of  HFC-236ea  were  also  investigated. 

METHODOLOGY 

Because  of  the  limitation  in  accurately  predicting  transport  properties  (such  as  thermal  conductivity)  by 
using  theory,  experimental  measurement  is  the  only  reliable  method  of  obtaining  viscosity  and  thermal 
conductivity.  There  are  a  number  of  ways  to  measure  viscosity  such  as:  the  capillary  viscometer  [10,11],  the 
vibrating  principle  viscometer  [12],  and  the  torsional  oscillation  viscometer  [13-16],  For  thermal  conductivity,  the 
vertical  cnavial  cylinder  method  [17,18],  transient  hot-wire  method  [19-23],  and  transient  hot-stnp  method  [24]  are 
the  usual  measurement  methods  employed. 

As  can  be  seen  from  the  above  list  of  methods,  each  property  is  measured  by  using  a  method  independent 
of  the  other  properties.  In  this  study,  a  new  approach  is  adopted  for  simultaneously  measuring  these  properties 
[25,26].  This  approach  uses  an  inline  viscometer  and  a  heat  transfer  test  section  to  dynamically  measure  viscosity 
and  thermal  conductivity  at  the  same  time.  Viscosity  is  measured  by  a  torsional  oscillation  inline  viscometer, 
while  thermal  conductivity  is  measured  from  knowledge  of  the  single-phase  heat  transfer  characteristics  of  a  heated 
test  section.  In  addition,  by  placing  a  mass  flow  meter  and  densimeter  in  series  with  the  viscometer  and  heat 
transfer  test  section,  density  and  specific  heat  can  also  be  determined. 

RESEARCH  PROGRAM 

The  main  goal  of  this  study  is  to  measure  transport  properties  of  HFC-236ea  with  and  without  lubricant  (Castrol 
oil  SW68)  by  using  a  new  methodology  [25,26].  Measured  data  for  pure  refrigerant  were  compared  with  other  theoretical 
data  sources  (e.g.,  REFPROP  [27] ).  The  test  facility  was  calibrated  by  using  the  well-known  properties  of  HCFC-22, 
CFC-12,  and  CFC-1 13,  and  the  methodology  was  verified  by  using  the  properties  of  CFC-1 14. 

# 

SCOPE 

The  scope  of  this  study  was  as  follows: 

•  Design  and  construct  a  test  rig  for  measuring  transport  properties,  such  as  thermal  conductivity  by 
using  an  approach  based  on  the  knowledge  of  single-phase  heat  transfer. 

•  Install  a  viscometer  in  series  with  a  heat  transfer  test  section  for  simultaneously  measuring  viscosity 
and  thermal  conductivity. 
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•  Calibrate  heat  transfer  and  heat  loss  characteristics  of  the  test  section  by  using  refrigerants  of  known 
properties. 

•  Verify  the  accuracy  of  the  viscometer  with  fluids  of  known  viscosity. 

•  Measure  the  viscosity  and  thermal  conductivity  of  CFC-1 14  and  compare  them  with  the  ASHRAE 
(American  Society  of  Heating,  Refrigerating  and  Air-Conditioning  Engineers)  standard  handbook  data. 

•  Measure  the  viscosity  and  thermal  conductivity  of  HFC-236ea  and  compare  them  with  theoretically 
predicted  data  (i.e.,  REFPROP). 

•  Measure  the  viscosity  and  thermal  conductivity  of  HFC»236ea  with  lubricant. 

•  Develop  prediction  equations  from  measured  properties. 

In  this  study,  density  (p)  and  specific  heat  (Cp)  of  the  refrigerants  of  interest  are  also  measured  and  discussed. 


CHAPTER  2 


CONCLUSIONS 


This  study  uses  a  new  approach  for  simultaneously  measuring  several  thermophysical  properties,  such  as 
thermal  conductivity,  viscosity,  specific  heat,  and  density.  This  approach  uses  single-phase  in-tube  heat  transfer 
knowledge  to  obtain  thermal  conductivity.  Viscosity  is  measured  by  a  viscometer  placed  in-line  with  the  heat 
transfer  test  section.  There  are  two  approaches;  Approach  1:  the  Nusselt  number  (Nu)  method  and  Approach  2:  the 
Prandtl  number  (Pr)  method.  The  uncertainty  analysis  was  presented  in  this  study.  Approach  2  seems  to  have  less 
uncertainty  than  Approach  1  and  T,  and  A T  do  not  significantly  affect  the  uncertainties.  However,  AT^  is  a 
significant  parameter  that  affects  the  uncertainties. 

In  Approach  1,  the  determination  of  a  calibration  function  (CF)  by  experiments  using  fluids  with  known 
properties  is  shown  to  be  important  for  accurate  thermal  conductivity  measurements.  Three  different  Nusselt 
number  correlations  were  used  for  calculating  thermal  conductivity  in  this  study,  and  they  have  been  examined  and 
discussed  herein.  Four  refrigerants,  HCFC-22,  CFC-12,  CFC-113,  and  CFC-114,  were  used  for  calibration  and 
verification  purposes.  These  refrigerants  cover  the  Pr  range  from  3  to  9  and  ReD  range  from  8,000  to  180,000. 
Based  on  the  calibration  results,  the  CF  functions  were  found  for  three  different  correlations  which  were  examined 
in  this  study.  In  Approach  2,  bypassing  the  Nusselt  number,  the  thermal  conductivity  was  found  from  Pr  which  is 
directly  related  to  ReD  and  non-dimensional  temperature,  AT*.  This  approach  was  shown  to  be  more  accurate  and 
convenient  to  use  because  less  variables  were  involved.  A  theoretical  uncertainty  analysis  also  showed  this 
approach  to  have  less  uncertainty  than  Approach  1.  The  measured  results  were  also  compared  and  discussed  for 
both  approaches,  and  consistency  was  shown  between  them. 

Viscosity  was  measured  by  a  torsional  oscillation  inline  viscometer.  The  accuracy  of  the  viscosity 
measurement  was  verified  with  CFC-113,  CFC-12,  and  pure  water,  and  shown  to  be  within  ±2%  when  compared 
with  the  ASHRAE  data.  The  measured  properties  included  specific  heat,  density,  viscosity,  and  thermal 
conductivity.  They  were  also  examined  for  CFC-1 14,  compared  with  ASHRAE  data,  and  shown  to  be  matched 
closely  within  ±5%  for  thermal  conductivity,  ±3%  for  specific  heat,  and  within  ±1%  for  density.  For  HFC-236ea 
property  m^r^ntc  REFPROP-4.0  data  were  used  as  a  comparison  with  the  measured  data.  The  deviations  of 
measured  properties  from  REFPROP-4.0  arc  +4.8%  for  specific  heat,  -5.0%  for  viscosity,  ±1%  for  density,  and 
+15%  for  thermal  conductivity. 
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A  lubricant  (Castrol  oil  SW68)  was  selected  to  mix  with  HFC-236ea.  Properties  were  measured  for  five 
lubricant  concentrations  over  a  temperature  range  of  -10  to  50°C.  Thermal  conductivity  effects  due  to  the  lubricant 
seemed  to  be  more  significant  at  low  temperatures  than  at  high  temperatures.  For  example,  at  a  low  temperature 
(-10°C)  the  thermal  conductivity  of  the  mixture  was  found  to  increase  by  40%  as  the  mixture  concentration 
changed  from  0.0%  oil  concentration  to  7.4%.  At  a  high  temperature  (50°C),  this  decrease  of  thermal  conductivity 
was  found  to  be  less  than  15%  at  a  for  oil  concentration  changing  from  0%  to  7.4%.  Viscosity  was  obviously 
affected  by  lubricant  concentration,  especially  at  a  low  temperature  and  high  lubricant  concentration  at  the  low 
temperature.  It  increased  over  50%  as  the  lubricant  concentration  went  from  0%  to  7.4%  at  the  low  temperature 
while  it  increased  nearly  100%  at  the  high  temperature.  Curve  fit  equations  for  both  one  variable  (temperature) 
and  two  variables  (temperature  and  lubricant  concentration)  were  provided  for  convenient  use.  Other  properties, 
such  as  density,  specific  heat,  thermal  diffusivity,  and  Prandtl  number,  were  also  calculated. 

Transport  properties  are  important  for  the  evaluation  of  refrigerant  alternatives.  The  method  developed  in 
this  study  shows  that  several  properties  can  be  quickly  and  simultaneously  measured.  Today  refrigerant 
alternatives  are  being  screened  and  tested  for  application.  The  methods  developed  in  this  study  can  provide  a  quick 
way  for  obtaining  these  properties. 
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CHAPTER  3 


RECOMMENDATIONS 


The  following  are  suggestions  for  further  improvement  of  the  current  study  and  possible  extended  future 
research  using  the  current  test  facility  and  methodology: 

1.  The  Prandtl  number  (Pr)  and  Reynolds  number  (Re)  ranges  of  calibration  can  be  further  extended  by 
testing  more  fluids  of  known  properties  to  extend  the  ranges  of  application  for  the  adopted 
methodology  in  this  study  so  that  wider  range  of  properties  can  be  further  evaluated. 

2.  The  accuracy  of  the  thermal  conductivity  and  viscosity  measurements  can  be  further  improved  by  using 
more  accurate  sensors. 

3.  Property  correlations  for  HFC-236ea  with  lubricant  can  be  further  developed  if  more  data  are  taken 
and  accumulated. 

4.  Although  properties  are  the  primary  interest  in  this  study,  the  single-phase  heat  transfer  and  flow 
characteristics  of  lubricant/refrigerant  mixtures  can  be  further  investigated  by  the  current  test  facility. 

5.  An  obvious  extension  of  this  study  would  be  to  use  other  potential  lubricants  mixed  with  HFC-236ea, 
subjected  to  a  wide  range  of  temperatures.  A  study  of  partly-miscible  lubricant-refrigerant  mixtures 
would  be  another  interesting  study. 
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CHAPTER  4 


REVIEW  OF  EXPERIMENTAL  METHODS 
FOR  THE  DETERMINATION  OF  TRANSPORT  PROPERTIES 


This  chapter  will  review  the  theory  surrounding  transport  properties  and  the  experimental  methods  used 
to  determine  transport  properties.  Several  traditional  methods  of  measuring  transport  properties  will  be  reviewed. 
Existing  error  sources  of  each  method  will  also  be  discussed.  Finally,  the  advantages  and  disadvantages  of  each 
traditional  method  will  also  be  compared. 

REVIEW  OF  THEORY  OF  TRANSPORT  PROPERTIES 

The  elementary  model  of  a  molecule  assumes  that  the  shape  of  the  molecule  is  like  a  rigid  sphere  of  diameter,  d. 
When  molecules  move  randomly  at  a  mean  velocity,  v,  they  collide  and  transfer  momentum  or  energy  assuming 
velocity  or  temperature  gradients  exist  These  momentum  or  energy  transfers  will  cause  an  exchange  of  momentum  or 
energy  flux  as  the  molecules  try  to  reach  an  equilibrium  status.  The  momentum  or  energy  flux  was  found  to  be 
proportional  to  the  velocity  or  temperature  gradient  The  corresponding  constants  of  proportionality  for  the  momentum 
and  energy  flux  equations  are  known  as  transport  properties.  Viscosity  is  associated  with  the  velocity  gradient 
while  thermal  conductivity  is  proportional  to  the  temperature  gradient  These  coefficients  are  also  proportional  to 
vLt 3  according  to  kinetic  theory  [28-3 1],  where  v  is  the  average  molecular  speed  and  L  is  the  mean  free  path  between 
molecular  collisions.  The  relations  between  these  properties  are  usually  described  by  the  following  equations  for 
low  density  gas: 


Momentum  flux 


dv  vL  dv 

=  ~Tnnrd2 


Energy  flux 


vL  JT 


3^n  dz 


(4.1) 

(4.2) 


where  n  is  the  number  of  molecules  in  a  unit  volume,  m  is  the  mass  of  the  molecule,  and  Cv  is  the  specific  heat 
under  constant  volume.  The  term  nmv  is  often  called  the  momentum  density  while  CvnT  is  the  energy  density. 
The  above  equations  define  the  transport  properties  in  the  momentum  and  energy  flux  equations,  which  are 
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viscosity,  p,  and  thennal  conductivity,  k,  respectively.  However,  if  the  average  speed  is  proportional  to  (RT!M)m 
and  the  free  path  to  (mf2)*1,  where  R  is  the  gas  constant,  M  is  the  molecular  weight,  and  T  is  the  temperature, 

then  the  viscosity  and  thermal  conductivity  can  be  rewritten  as: 


mpvL  TxriMm  .. .. 

H  ^ —  «  ( constant )  ^  (4  3) 

k  =  =  (constant) 

where  p  is  the  mag*  density.  Another  theoretical  equation  for  the  thennal  conductivity  of  liquids,  which  was  also 
basfd  on  a  rigid  sphere  model,  was  proposed  by  Bridgman  in  1923  and  later  modified  by  Powell  et  al.  [31].  The 
resulting  equation,  which  has  been  used  for  pure  substances,  is: 

k  =  l.&^yiaVs  (4.5) 


where  N  is  Avogadro's  number,  V  is  the  molar  volume,  o  is  the  Boltzmann's  constant,  and  IS  is  the  speed  of  sound. 

The  above  equations  are  based  on  ideal  models  which  neglect  an  interaction  force  that  might  exist 
between  molecules.  As  a  result,  a  rigid  sphere  model  which  neglects  forces  between  molecules  is  no  longer 
acceptable  and  must  be  modified.  For  different  molecules,  the  modification  will  be  different. 

Although  a  number  of  prediction  rules  were  developed,  there  is  still  no  universal  principle  which  can  be 
used  to  cover  every  group  of  fluids,  especially  when  applied  to  mixtures.  Baroncini  et  al.  (1981)  developed  a 
general  correlation  for  organic  liquid  thennal  conductivity  [32].  Although  it  provided  useful  information  for  some 
existing  materials,  it  still  needs  further  study  for  new  fluids  or  refrigerant  mixtures.  Kandlikar  et  al.  (1975)  also 
published  a  paper  to  theoretically  predict  viscosity  and  thermal  conductivity  for  HCFC-22  and  CFC-12  mixtures 
[33].  Later,  Levy  (1981)  used  a  modified  Maxwell-Eucken  equation  to  calculate  thermal  conductivity  for  two 
component  solution  mixtures  [34],  Ely  and  Hanley  (1983)  used  a  corresponding  state  method  to  predict  the 
viscosity  and  thermal  conductivity  of  mixtures  [35]. 

Although  the  above  studies  provide  useful  information,  they  still  have  some  limitations  in  actual  use, 
especially  for  multi-component  mixtures.  This  is  because  the  interaction  force  between  molecules  is  still  not  fully 
understood.  and  some  correlations  need  semi-empirical  data  from  experiments  to  modify  the  models.  However,  one  of 
the  mixing  rules  has  been  widely  used  in  past  studies  [36,37].  This  rule  is  described  by  the  equation  shown  below: 

km  =  kr  (1  -  Q  +  kaC  -  0.72  (kD  -kr)(\-QC  (4.6) 

where  C  represents  the  lubricant  concentration. 
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Lin  and  Pate  (1992)  used  thermodynamic  properly  data  and  the  speed  of  sound  of  the  mixture  to  predict 
the  thermal  conductivity  of  refrigerant  mixtures  by  using  Bridgman  and  Powell's  equation  (Equation  4.5)  to 
develop  a  mixing  rule  [38].  However,  this  method  needs  modified  coefficients  which  requires  calibrations  and 
numerous  tests  with  known  fluid  properties. 

Although  some  theoretical  prediction  methods  and  limited  experimental  data  for  obtaining  properties  of 
new  pure  refrigerants  and  refrigerant  mixtures  have  been  published  [39-43],  there  are  still  very  limited  experimental 
data  for  verification,  especially  for  lubricant/refrigerant  mixtures.  Therefore,  experimental  methods  for  obtaining 
the  required  data  become  vital  in  the  current  CFC  transition. 

EXPERIMENTAL  METHODS  FOR  THE  DETERMINATION  OF  THERMAL  CONDUCTIVITY 

Thermal  conductivity  can  be  measured  by  experimental  methods  [44,45]  by  using  Fourier's  law.  The 
methods  developed  for  measuring  thermal  conductivity  can  be  divided  into  the  following  two  categories:  steady- 
state  methods  and  unsteady-state  methods.  The  following  is  a  brief  review  of  these  methods. 

Steadv-State  Methods 

Steady-state  methods  differ  primarily  in  the  different  geometries  used  in  the  apparatus.  In  addition,  all 
steady-state  methods  require  careful  experimentation  to  minimize  errors  for  heat  losses  and  to  prevent  the  onset  of 
any  convective  motion  in  the  fluid. 

1.  Horizontal  Flat-Plate  Method- 

In  this  method,  the  heat  passes  through  a  layer  of  fluid  located  between  two  circular  horizontal  metal 
plates.  With  a  suitable  design  along  with  radiation  and  convection  corrections  in  the  measurement,  the  thermal 
conductivity  accuracy  can  be  optimized.  This  method  measures  the  temperature  difference  across  the  fluid  layer 
and  the  heat  flow  rate  with  known  geometric  dimensions.  The  accuracy  of  thermal  conductivity  measurement  by 
using  this  method  depends  on  the  precision  of  the  temperature  measurement  and  the  accuracy  of  the  geometric 
dimensions  of  the  cell. 

2.  The  Vertical  Coaxial  Cylinder  Method- 

In  this  method,  two  coaxial,  vertical  cylinders  form  a  test  cell.  Fluid  fills  the  annulus  of  the  cylinders,  and 
heat  is  generated  in  the  center  of  the  inner  cylinder.  By  measuring  the  heat  passing  through  the  test  fluid  and  the 
temperatures  at  the  inner  and  outer  surfaces  of  the  two  cylinders,  the  thermal  conductivity  can  be  determined. 
Again,  the  experimental  uncertainty  needs  to  be  minimized  for  improving  accuracy. 
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3.  Hot-Wire  Method— 

The  hot-wire  method  is  in  principle  a  special  case  of  the  coaxial  cylinder  method  in  which  the  inner 
cylinder  is  replaced  by  a  thin  wire.  The  wire  not  only  performs  as  a  heater,  but  also  performs  as  an  internal 
thermometer.  The  temperature  of  the  outer  surface  of  the  fluid-filled  annulus  can  be  evaluated  from  the 
dimensions  heat  flow  rate,  thermal  conductivity  of  test  fluid,  and  the  outer  surface  temperature  of  the  outer 
cylinder.  The  heat  flow  generated  in  the  wire  is  determined  from  the  electric  power  added  to  wire.  Correction  for 
the  heat  loss  is  required  for  an  accurate  heat  flow  determination.  The  thermal  conductivity  is  determined  from  the 
measured  quantities  and  the  dimensions  of  the  cell. 

4.  Concentric  Sphere  Method- 

In  this  method,  two  concentric  spheres  make' up  the  test  cell.  The  test  fluid  is  contained  in  the  volume 
between  the  spheres.  A  heater  located  in  the  cavity  of  the  inner  sphere  generates  heat  that  is  conducted  radially 
outward  to  a  concentric  spherical  shell.  Temperatures  on  the  inner  and  outer  surfaces  of  the  spheres  are  measured 
by  sensors  embedded  in  the  surface  walls.  With  known  dimensions,  the  heat  generated  in  the  inner  sphere,  and 
temperature  measurements,  the  thermal  conductivity  can  be  obtained. 

Unsteady  Methods 

The  establish™^  of  steady  state  in  an  apparatus  of  considerable  mass  is  a  slow  process,  and  measurements 
may  be  rather  time  consuming.  It  is  therefore  natural  to  explore  the  possibility  of  employing  transient  techniques 
in  which  the  same  information  may  be  obtained  while  the  fluid  system  is  subjected  for  a  short  time  to  the  influence 
of  a  transient  temperature  field. 


1.  Continuous  Line  Source- 

A  transient  heat  transfer  characteristic  of  a  line  source  with  an  infinite  length  and  a  constant  energy 
generation  in  an  extended  isotropic  fluid  can  be  used  to  determine  the  fluid  thermal  conductivity.  From  the 
measurement  of  the  temperatures  at  a  fixed  radial  distance  at  two  different  times,  the  thermal  conductivity  can  be 
determined  by  the  following  equation: 


Q  /"( W 

4*  (^2,oo  _  V 


(4-7) 


where  Q  is  the  total  heat  added  to  the  line  source  between  times  /j  and  t2 ,  and  7"1>o0  and  T2ko  are  the  measured 
fluid  temperatures  at  a  fixed  radial  distance  from  the  heat-dissipated  line  source  at  two  different  times,  fj  and  t2, 
respectively.  A  plot  of  temperature  versus  the  natural  logarithm  of  time  should  result  in  a  straight  line.  The 
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unique  feature  of  the  solution  from  this  method  is  that  it  yields  thermal  conductivity,  k,  directly.  Recently,  the 
transient  hot-wire  method  has  been  applied  to  the  measurement  of  liquid  thermal  conductivity  [46-50],  This 
method  uses  a  veiy  thin  wire  as  a  line  source  in  which  heat  is  generated  electrically.  The  thermal  conductivity,  k, 
of  the  test  liquid  is  determined  by  measuring  heat  input,  temperature  rise  of  the  wire,  and  the  time  of  heating.  The 
equation  [51,52]  is  given  as  follows: 


in  which  q  is  the  heat  input  per  unit  length  of  the  wire,  k  is  the  liquid  thermal  conductivity  surrounding  the  wire,  t  is  the 
time  of  hearing,  r0  the  wire  radius,  and  C  a  numerical  constant  (C » ef),  y = 0.5772157-,  Euler's  constant).  The  symbol  a 
represents  the  thermal  diffusivity  of  the  liquid  surrounding  the  wires. 

2.  Cylindrical,  Spherical,  and  Plane  Sources- 

A  different  approach  from  that  used  in  the  transient  hot-wire  method  is  based  on  a  cylindrical,  spherical, 
or  plane  source  taking  the  place  of  the  line  source  (wire).  In  these  methods,  heat  is  generated  electrically  for  a 
short  period  of  time  in  a  relatively  large  solid  body  of  a  simple  geometry  (cylinder,  sphere,  plane).  After  the 
temperature  of  this  body  has  been  raised  by  a  few  degrees  above  that  of  the  heat  receiving  body,  which  is  assumed 
to  remain  at  a  constant  temperature  in  its  isothermal  surroundings,  heating  is  discontinued.  Heat  is  now  conducted 
through  the  narrow,  fluid-filled  gap  between  the  two  geometrically  similar  bodies  and  the  change  of  temperature 
with  time  is  recorded.  The  thermal  conductivity  is  then  calculated  using  the  equations  which  describe  the  transient 
heat  transfer  characteristics  for  different  geometries. 

ERROR  SOURCES  FOR  EXPERIMENTAL  METHODS 
FOR  THE  DETERMINATION  OF  THERMAL  CONDUCTIVITY 

Despite  precautions  and  sophistication  in  the  test  system's  design  and  of  the  sensors,  it  is  impossible  to 
completely  eliminate  all  sources  of  error  for  the  measured  quantities.  There  are  several  error  sources  which  exist 
in  the  methods  mentioned  above;  however,  two  sources  of  error,  namely  radiation  and  convection,  are  the  major 
contributors  affecting  the  energy  transfer  through  the  fluid  film. 

Radiation  Effect  in  Parallel  Plates  or  Cylinders 

Thermal  radiation  in  partly  transparent  fluids  affects  the  heat  transfer  in  two  ways:  first,  the  energy  is 
directly  absorbed  and  emitted  by  the  test  fluids,  and  second,  the  energy  is  absorbed  and  emitted  by  the  walls  or 
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surfaces  of  the  test  cell  or  gap.  Thus,  the  radiation  heat  transfer  affects  the  conduction  heat  exchange  determined 
from  Fourier's  Law  which  is  the  basis  for  the  thermal  conductivity  measurement  methodology.  Since  Fourier's  Law 
does  not  account  for  radiation,  the  thermal  conductivity  measurement  needs  to  be  corrected  for  radiation  heat 
transfer.  The  temperature  distribution  depends  not  only  on  the  geometry  of  the  cell  system,  but  also  on  the 

temperature  of  the  surfaces  surrounding  the  test  fluid. 

Leidenfrost  analyzed  the  thermal  radiation  effects  on  the  thermal  conductivity  measurement  in  both 
vertical  and  horizontal  parallel  plate  cells  [53].  Leidenfrost  calculated  the  ratio  of  the  heat  transferred  by  radiation 
to  the  total  heat  transfer  and  found  that  this  ratio  was  strongly  related  to  emissivity.  Table  4.1  was  obtained  from 
the  calculated  results  shown  in  Figure  10  of  Leidenfrost's  publication  [53]. 

TABLE  4.1.  RADIANT  PERCENTAGE  BETWEEN  PARALLEL  WALLS 


Emlssivitv',  e*  0,1 

Emlssivitv.  e«=i.O 

273 

0.1% 

1.0% 

383 

0.2% 

2.0% 

483 

0.5% 

4.1% 

600 

1.0% 

8.8% 

example  for  toluene  at  AT  =  5  K 

From  the  radiant  heat  transfer  percentages  shown  in  Table  4. 1,  it  was  noted  that  the  radiation  heat  flow 
affects,  the  heat  conduction  which,  in  turn,  affects  the  accuracy  of  the  thermal  conductivity  measurement.  The 
higher  the  temperature  and  emissivity,  the  larger  the  radiant  heat  exchange.  At  a  temperature  of  600  K  for 
toluene,  the  radiation  heat  transfer  ratio  can  even  reach  8.8%  if  black  bodies  are  used.  This  result  suggests  that  if 
no  corrections  for  radiant  heat  exchange  are  applied  to  the  measurements,  then  the  values  of  thermal  conductivity 
will  be  too  high  by  the  same  percentages. 

As  a  further  step  in  understanding  how  radiation  affects  conduction  in  flat  plate  and  concentric  cylinder 
cells,  a  theoretical  analysis  is  performed  below.  The  situation  is  one  in  which  the  test  fluid  is  completely  transparent  to 
radiation  of  all  wavelengths.  For  this  situation,  the  heat  transferred  by  radiation  can  be  calculated  as  follows. 

1.  Parallel  Flat  Plate  Case- 

The  equations  governing  the  radiation  and  conduction  heat  transfer  between  two  parallel  flat  plates  are 
given  as: 

*5  (4-9) 
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where  e1>2  is  the  radiation  interchange  factor,  oJ  is  the  Stefan-Boltzmann  constant,  A  is  the  heat  transfer  area,  and 
7j,  T2  represent  the  surface  temperatures  of  these  two  flat  plates.  For  two  parallel  walls,  e12  is  given  as: 


(4.10) 


However,  the  total  heat  transferred  by  conduction,  which  is  the  quantity  used  to  calculate  thermal  conductivity  in 
this  apparatus,  is  based  on  Fourier's  Law  and  expressed  as: 


_kA{Tx-  T2) 

Qt  Ax 


(4.11) 


where  k  is  thermal  conductivity  of  test  fluid  and  Ax  is  the  normal  distance  between  emitter  and  receiver  plates  or 
the  distance  of  test  cell  gap.  The  ratio  of  the  radiation  heat  transfer  to  the  total  heat  transfer  is  then  calculated  by 
dividing  these  two  equations,  which  is: 

| - (4,2) 

a.  +  — -i) 

It  is  of  interest  to  calculate  this  ratio  in  order  to  determinate  the  significance  of  the  radiation  heat  transfer 
in  relation  to  the  total  heat  transfer.  For  this  example  CFC-114  was  selected  as  the  test  fluid.  Figure  4.1  plots  the 
radiant  heat  transfer  ratio  versus  temperature  for  this  case.  As  shown  in  the  figure,  the  radiation  heat  transfer  is 
quite  significant  at  high  temperatures  and  high  emissivities  of  the  plates.  This  ratio,  for  example,  can  even  reach 
+15%.  In  other  words,  the  thermal  conductivity  measured  under  this  situation  will  be  15%  higher  than  the  true 
value.  Viskanta  and  Grosh(1962)  calculated  the  simultaneous  conduction  and  radiation  in  an  absorbing  medium 
between  parallel  flat  plates  (54],  The  calculated  results  showed  that  the  radiation  effect,  which  depend  on  the 
distance  between  the  two  plates,  was  of  the  same  order  of  magnitude  as  the  conduction. 

2.  Concentric  Cylinder  Case— 

The  equation  which  governs  the  radiation  heat  transfer  in  between  two  long  concentric  cylinders  is  given  as: 

Cr-«U»^W-yJ)  (413) 

where,  again,  e1>2  is  the  radiation  interchange  factor,  ct  is  the  Stefan-Boltzmann  constant,  /  is  the  length  of  each 
cylinder,  r,  is  the  center  radius  of  the  inner  cylinder,  and  r2  is  the  inner  radius  of  the  outer  cylinder  (r,  <  r2).  For 
this  case,  e1|2  is  given  as: 
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(4.14) 


In  this  case,  the  total  heat  transferred  by  conduction  is  expressed  as: 


(4.15) 


The  ratio  of  radiation  to  the  total  conduction  heat  transfer  is  then  calculated  by  dividing  these  two  equations 
(Equations  4. 13  and  4. 15),  resulting  in 


4 

Q, 


(7?-72)/n(r2/rl> 


(4.16) 


This  ratio  was  calculated  by  using  CFC-1 14  as  a  test  fluid.  Figure  4.2  shows  the  percentage  of  the  ratio  erf  radiant 
heat  transfer  to  total  heat  transfer  versus  temperature  for  CFC-1 14  in  the  concentric  cylinder  case.  Compared  with 
the  parallel  flat  plate  case  shown  in  Figure  4.1,  the  radiant  ratio  for  the  concentric  cylinder  case  is  quite  similar. 
Again,  the  radiant  heat  transfer  is  quite  significant  at  high  emissivities  and  high  temperatures. 

In  summary,  the  errors  in  thermal  conductivity  measurements  introduced  by  radiant  heat  exchange  in  both 
parallel  flat  plate  and  concentric  cylinder  cells  should  not  be  assumed  to  be  negligibly  small  until  verified.  In  order 
to  achieve  the  desired  accuracy,  corrections  should  be  made  and  the  material  selected  for  the  test  cell  walls  should 
be  appropriate  so  that  the  errors  can  be  minimized. 


Radiation  Effect  on  Hot-Wire  Cell 

In  the  transient  hot-wire  method,  the  essential  relationship  is  that  of  the  temperature  rise  of  the  thin  wire 
imnipro-H  in  the  fluid  as  a  function  of  time  following  the  stepwise  initiation  of  a  heat  flux  within  it  There  have 
been  a  number  of  attempts  to  cany  out  an  analysis  of  the  process  of  simultaneous  conduction  and  radiation  in  an 
absorbing  fluid  for  the  transient  hot-wire  instrument  [55-57],  Menashe  and  Wakeham  analyzed  the  combined 
radiation  and  conduction  heat  exchange  in  a  transient  hot-wire  cell  in  1982  [56].  They  theoretically  calculated  the 
apparent  thermal  conductivity  measured  by  this  method  and  drew  the  conclusion  that  the  apparent  thermal 
conductivity  ratio,  kapp  /  k,  is  proportional  to  the  equilibrium  temperature,  T0, 

kapp/k  oc  T\  (4.17) 

They  also  calculated  the  kapp  /  k  for  n-undecane  at  348  K  and  found  that  the  apparent  thermal  conductivity  of  the 
fluid  was  as  much  as  2.5%  greater  than  the  true  value. 
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Figure  4.1:  Radiant  heat  transfer  percentage  in  two  parallel  flat  plates  using  CFC-114  as  the  test  fluid. 
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Figure  4.2:  Radiant  heat  transfer  percentage  in  concentric  cylinder  cell  using  CFC-114  as  the  test  fluid. 


Saito  and  Vcnart  analyzed  the  radiation  effect  on  a  transient  line  source  measurement  by  a  theoretical 
approach  157].  They  analyzed  the  ratio  or  apparent  thermal  conductivity  (radiation  involved)  to  true  thermal 
conductivity  and  found  it  was  a  function  of  radiation  properties  such  as  emissivity  and  refractivity  plus  equilibrium 
temperature  and  geometric  parameters,  such  as  the  radius  of  the  line  source.  The  results  are  tabulated  in  Table  4.2 
(values  were  cited  from  Figure  6  of  Reference  [57] ). 

TABLE  4.2:  +  VALUES0  FOR  APPARENT  THERMAL  CONDUCTIVITY  RATIO 


II  III  ‘ 

SliWoIil# 

R*O.0OJ 

e*0.2 

2.750 

2.125 

0.750 

e  =  0.4 

3.188  ‘ 

2.500 

1.250 

e  =  0.6 

3.625 

3.000 

1.750 

O 

II 

O 

oo 

4.063 

3.500 

2.313 

s  —  1.0 

4.500 

4.000 

2.813 

<|>  is  defined  in  Equation  4.18  and  R  is  defined  in  Equation  4.19 
a  values  were  cited  from  Figure  6  of  Reference  [57] 

In  Table  4.2,  the  non-dimensional  group  (<t»)  as  a  function  of  apparent  thermal  conductivity  ratio  is  defined  as: 

$  =  (kappfk- l)16ri^57J  (418) 

R  =  rl’ka  (4.19) 

where  rj  is  the  radius  of  line  source,  n  is  the  index  of  refraction,  and  ka  is  absorption  coefficient,  m*1. 

Table  4.2  shows  the  <j>  values  (read  from  Figure  6  of  Ref.  [57] )  for  various  e  and  R.  Because  4>  is 
proportional  to  kapp  /  it,  an  increasing  tendency  of  kapp  /  k  with  larger  e  and  R  was  obtained  from  Table  4.2.  These 
results,  as  expected,  have  the  same  trends  for  both  the  parallel  flat  plate  and  the  concentric  cylinder  arrangements. 

In  Equation  4.18,  the  equilibrium  temperature,  T0,  is  another  quantity  that  affects  the  radiation  strength. 

A  plot  of  kapp  /  k  versus  T0  is  interesting  to  view.  Applying  Equation  4. 18  to  a  test  condition,  the  kapp  /  k  ratio 
can  be  shown.  Figures  4.3  through  4.5  show  plots  of  kapp  /  k  versus  k  for  various  T0  values  at  e  =  0.2,  e  =  0.6,  and 
e  =  1.0,  respectively. 

As  shown  in  these  figures,  the  ratio  of  kapp  /  k  generally  increases  with  increasing  temperature,  T0,  at  the 
smaller  k  values.  In  addition,  in  the  lower  thermal  conductivity  range,  this  ratio  becomes  more  significant  than  it 
is  at  the  higher  thermal  conductivity  ranges.  Figure  4.6  shows  the  kapp  /  k  versus  fcata  high  temperature  of  T0  = 
120°C  for  various  emissivities.  As  expected  in  this  case  of  higher  emissivity,  radiation  heat  transfer  contributes  more 
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to  the  total  conduction  heat  transfer  so  that  the  ratio,  kapp /k,  becomes  higher.  Figure  4.7  shows  the  kapp /k 
versus  temperature  at  e  =  1.0  for  different  thermal  conductivity  ranges.  It  indicates  that  the  higher  the 
temperature,  the  higher  the  ratio,  and  for  low  thermal  conductivity  fluids,  this  tendency  would  be  more  significant. 
These  results  show  that  in  the  transient  hot-wire  method  the  temperature  and  emissivity  of  the  wire  are  two 
important  factors  affecting  the  thermal  conductivity  measurement.  C.A.  Nieto  de  Castro  et  al.  (1991)  [58]  stated 
that  the  contribution  of  radiative  heat  transfer  to  thermal  conductivity  was  found  to  range  between  2  and  50%. 

C.A.  Nieto  de  Castro  et  al.  investigated  toluene  and  dimethylphthalate  liquid  thermal  conductivity  [59]. 
They  indicated  that  the  radiant  effect  in  determining  the  thermal  conductivity  using  the  transient  hot-wire  method 
is  quite  significant  without  any  corrections.  Gross  et  al.  (1992)  [60]  measured  the  thermal  conductivities  of 
refrigerants  HFC-134a,  HFC-152a,  and  CFC-123  using  the  hot-wire  transient  method.  Although  the  accuracy  was 
claimed  within  ±1.6  to  -2.0%  with  proper  corrections,  the  uncertainties  were  attributed  to  limitations  in  the  ideal 
mathematical  model  of  a  line  source  in  infinite  space  and  the  effects  of  convection  and  radiation. 

Another  factor  which  affects  the  thermal  conductivity  measurement  error  is  the  wire  diameter.  From 
Saito's  apparent  thermal  conductivity  ratio  equation,  Equation  4.18,  the  apparent  thermal  conductivity  ratios  were 
dominated  by  the  diameter  (2 r  j)  of  the  thin  wire  used.  Table  4.3  indicates  how  the  ratio  was  affected  by  the  wire 
diameter  used. 


TABLE  4.3:  RATIO  OF  APPARENT  THERMAL  CONDUCTIVITY 
TO  ACTUAL  THERMAL  CONDUCTIVITY 


1 1  diameter  (2rj ),  pm  llli§|f 

5.0 

1.003 

10.0 

1.006 

50.0 

1.030 

100.0 

1.060 

1 .  kapp  /  k  was  calculated  by  Equation  4.18 

2.  Cases  for  toluene  with  the  following  parameters: 
equilibrium  temperature,  Tn  =  360  K 

R  =  0.01 

thermal  conductivity,  k  =  0. 1 137  Wl(m  •  K ) 
refractive  index,  n  =  1.4961 

Stefan-Boltzmann  constant,  crs  =  5.6697  x  10*8  WKm2  *  K4  * 

As  shown  in  Table  4.3,  the  radiation  effect  is  affected  by  the  wire  diameter  and  is  never  completely  eliminated 
although  minimized  by  using  a  thin  wire.  However,  the  experimental  uncertainty  of  the  smaller  diameter  wire 
would  be  larger  than  that  for  a  larger  diameter  wire.  Therefore,  an  appropriate  wire  diameter  must  be  selected  in 
order  to  achieve  optimal  accuracy  by  this  method. 
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Figure  4.6:  kapp  I  k  ratio  as  a  function  of  k  for  various  e  at  high  temperature  case 
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Figure  4.7:  kapp  /  k  ratio  as  a  function  of  temperature  for  different  thermal  conductivity  ranges. 


Convection  Effects  on  Thermal  Conductivity  Measurement 


The  design  of  an  experimental  apparatus  for  the  determination  of  the  thermal  conductivity  of  a  liquid 
should  preclude  the  possibility  of  free  convection  [61,62].  Tsedeiberg  [45]  mentioned  Kraussold's  study  of  the 
conditions  at  which  free  convection  occurs  for  various  liquids  located  between  two  concentric  cylinders  when  the 
thermal  conductivity  of  the  fluids  was  known.  Kraussold  found  that  free  convective  heat  transfer  arises  when  (Gr  * 
Pr)  £  1000.  Kraussold,  indicating  that  both  horizontal  and  vertical  placements  of  the  tube  were  similar,  analyzed 
the  apparent  thermal  conductivity  ratio  for  concentric  cylinders  and  developed  the  following  correlation: 

*2EEl±  =  0.11  (Gr  *  Pr)0-29  (4.20) 

This  equation  is  only  valid  for  3.8  £  (Gr  •  Pr)  <;  6.0.  Instead  of  Kraussold’s  equation.  Equation  4.20,  Wakeham  et 
al.  [51]  used  a  different  equation  to  estimate  the  heat  transfer  by  convection,  which  for  a  cylindrical  geometry  is: 

Qcohv  _  dRa  ..  ... 

Qw=  720/ 

where  d  is  the  thickness  of  the  fluid  layer,  /  the  length  of  the  internal  cylinder,  and  Ra  =  (Gr  *  Pr)  is  the  Rayleigh 
number.  Results  from  using  Equation  4.21  are  shown  in  Table  4.4. 

TABLE  4.4:  AN  EXAMPLE  OF  CONVECTION  EFFECTS  ON 
CONDUCTION  IN  CYLINDRICAL  GEOMETRY 


tWffiiiOTKi'rent  thermal  conductivity  ratio  (based 

§j|  -Pr'"  (§ 

$«*• p: 

(d  mm} 

1000 

0.0023 

0.0058 

5000 

0.0116 

0.0289 

10000 

0.0231 

0.0579 

From  Table  4.4,  it  is  noted  that  the  significance  of  free  convection  depends  on  the  Rayleigh  number,  Ra.  For  large 
Ra,  the  convection  effect  would  be  large  if  the  fluid  layer,  d,  is  large.  Therefore,  control  of  Ra  is  necessary  in  order 
to  decrease  the  effect  of  free  convection  on  total  heat  transfer.  However,  the  measured  thermal  conductivity  still 
needs  to  be  corrected  for  free  convection  effects.  For  a  non-steady  method,  such  as  the  transient  hot-wire  method, 
Tye  [44]  mentioned  and  suggested  that  there  might  be  a  corresponding  pseudo-steady  state,  and  convection  would  be 
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expected  to  occur  when  the  pseudo-steady  state  temperature  reaches  a  critical  value  defined  by  RaEGr  ■  PrZ.  1000 
[44].  A  criterion  for  the  onset  of  convection  in  a  non-steady  hot-wire  cell  was  suggested: 

(  )  7*  (  e0  >T>  -  1  y  <Z  1000  (4.22) 

where  T*  =  4nkTga  /  Q,  p  is  dynamic  viscosity  of  liquid,  Cw  is  the  wire  heat  capacity,  is  fluid  temperature,  Q  is 
the  heat  input,  and  a  is  thermal  diffusivity.  A  plot  of  t*  (  =  A  Cw/r^)  against  T*  ( e0ST*  -  1 )  can  be  used  to 
determine  the  true  time,  f,  for  the  onset  of  convection,  and  the  necessary  error  correction  can  be  made. 

Summary  of  Thermal  Conductivity  Measuring  Methods 

The  methods  mentioned  above  for  experimental  thermal  conductivity  determination  are  summarized  as: 
the  steady-state  method  and  the  non-steady-state  method.  The  main  difference  between  the  steady-state  and 
unsteady  method  is  the  amount  of  time  consumed.  The  steady-state  method  requires  more  time  to  reach  steady 
state,  which  affects  the  accuracy  of  measurement  because  the  mathematical  model  is  based  on  the  real  steady-state 
conditions,  whereas  the  unsteady-state  method  allows  prompt  measurements  to  be  taken.  Both  methods  demand 
careful  experimentation  to  minimize  corrections  for  heat  losses  and  to  prevent  the  onset  of  convective  motion  in  the 
fluids.  However,  other  effects  such  as  time  counting  deviation,  geometric  errors,  temperature  sensor  errors,  and 
constant  surface  temperature  controlling  scheme  also  contribute  to  the  overall  uncertainty  in  values  of  thermal 
conductivity  for  both  methods.  Table  4.5  is  a  brief  discussion  of  those  steady-state  methods  mentioned  above,  and 
Table  4.6  describes  the  advantages  and  disadvantages  for  unsteady-state  methods. 

The  uncertainty  sources  are  generally  known  for  all  of  these  methods.  Typically,  the  main  uncertainties  have 
been  catalogued  into  two  sources:  radiation  error  and  convection  error.  Table  4.7  shows  the  errors  of  these  main  sources. 
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TABLE  4.5:  SUMMARY  OF  ADVANTAGES  AND  DISADVANTAGES  OF  STEADY-STATE 

METHODS  OF  THERMAL  CONDUCTIVITY  MEASUREMENTS 


horizontal,  flat-plate  gg 
method  and  vertical,  coaxial 

cylinder  method  m 

hot-wire  method  g| 

concentric  sphere 
method 

advantages  ; 

*  static  fluid  -  no  flow 

•  less  fluid  properties  related 

•  static  fluid  -  no  flow 

•  simplicity  and  ease  of 
construction 

•  less  fluid  properties  related 

•  static  fluid  -  no  flow 

*  no  other  properties  related 

disadvantage* 

•  difficulties  in  eliminating 
convection  pertuibation 
along  the  edges  of  the 
heated  plates 

•  difficulties  in  controlling 
isothermal  conditions  on 
the  surfaces  of  the  plates 

•  radiation  effect  not  negligible 

•  real  steady  state  attainment 

•  time  consuming 

•  axial  conduction  effect  of 
the  wire 

•  precise  determination  of 
the  wire  and  tube  diameters 

«  accurate  coaxial  alignment 
of  the  wire  in  the  tube 

•  difficulty  in  apparatus's 
construction  and  accurate 
alignment 

•  effect  of  convection 

•  difficulty  in  real  surface 
temperature  measurement 

approximate 

±1.0% -±15% 

±1.0% -±18% 

±1.0%  ~  ±15% 

TABLE  4.6:  SUMMARY  OF  ADVANTAGES  AND  DISADVANTAGES  OF  UNSTEADY-STATE 
METHODS  OF  THERMAL  CONDUCTIVITY  MEASUREMENTS 


• "  source  •• 

*  less  time  consuming 

*  no  thermal-resistance 
layer 

•  less  time  consuming 

*  simple  construction 

disadvantages  . 

•  effect  of  finite  length 
(originally  assumed  as 
infinitely  long) 

•  effect  of  heat  capacity 
depending  on  temperature 

•  effect  of  bounding  wall 

•  onset  of  convection  effect 
(Ra  £  1000)  -  when  time 
going 

•  radiation  effect  never 
completely  eliminated  - 
correction  is  difficult 

•  difficulty  in  apparatus's 
construction  and  accurate 
alignment 

•  effect  of  heat  capacity 

•  effect  of  thermal  boundary 
layer 

•  onset  of  convection  effect 
(Ra Z  1000) 

•  radiation  effect  never 
completely  eliminated  - 
correction  is  difficult 

|||^roxiniate  j 
errorrange 

±1.0% -±18% 

±1.0% -±15% 
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TABLE  4.7:  SUMMARY  OF  TYPICAL  UNCERTAINTIES  FOR  EXPERIMENTAL 
THERMAL  CONDUCTIVITY  MEASUREMENTS 


method 

radiation 

H  ■  convection  | 

^  - -  total  error 

steadystaie 

concentric  cylinders 

emissivity  and  temperature 
dependent  -  low  temp.: 

±0.1%  -±1.0% 

(s*  0.1 -1.0) 

high  temp.:  ±1.0% -±8.8% 

(e  =  0.1 -1.0) 

negligible  for  Ra  <  1000 
±0.3% -±4.2% 

Ra:  1000-10,000 

±0.4% -±13% 

;^p:non4*eady  gitatK^ii! 
transient 
' M  method  • 

emissivity,  wire  size,  and 
temp,  depending  on  - 
e  effect:  ±1.0% -±10% 
wire  diameter:  ±0.3%  -  ±6% 
temp,  effect:  ±1%  —  ±10% 

less  than  ±1% 

±3% -±26% 

EXPERIMENTAL  METHODS  FOR  VISCOSITY  MEASUREMENT 


Some  methods  usually  used  for  the  viscosity  measurements  include  the  capillary  tube  [10,1 1,63,64,65], 
falling  cylinder  [66,67],  rolling  ball  [68,69],  vibrating-wire  method  [12],  and  torsional  oscillation  method  [14-16], 
The  following  is  a  brief  introductory  review  of  these  methods. 


Capillary  Viscometer 

The  capillary  viscometer  [10,11]  uses  the  flowing  fluid  characteristics  in  a  capillary  tube  of  a  certain 
length.  The  viscosity  is  related  to  the  following  quantities:  volumetric  flow  rate,  pressure  drop,  density  of  fluid, 
thermal  expansion  coefficient  of  the  capillary,  kinetic  energy  correction  factor,  and  the  capillary  constant.  The 
capillary  constant  is  determined  using  a  known  viscosity  fluid,  usually  pure  water.  The  viscosity  is  calculated  from 
a  modified  Hagen-Poiseuille  equation. 

Vibrating-Wire  Viscometer 

The  vibrating-wire  viscometer  [12]  is  used  for  Newtonian  fluids.  A  circular-section  beam  of  constant 
length  with  both  ends  clamped  at  fixed  points  is  subjected  to  a  working  fluid.  The  beam  is  set  into  oscillation  in  a 
single  mode  within  a  plane  containing  the  axis  of  the  beam  and  perpendicular  to  it  by  means  of  an  initial 
displacement  at  time  t  =  0  in  an  initially  stationary  fluid  of  constant  viscosity  and  density.  The  viscosity  is 
determined  tty  the  relationship  between  the  frequency  and  the  magnitude  of  the  oscillation  of  the  beam  or  wire  in  a 
fluid.  A  correction  is  required  for  the  specific  wire. 
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Torsional  Oscillation  Viscometer 


This  kind  of  viscometer  uses  the  principle  of  surface  loading,  that  is,  a  vibrating  surface  in  contact  with  a 
liquid  experiences  a  force  which  is  a  function  of  viscosity  [14].  The  oscillatory  surface  may  be  spherical,  cylindrical, 
or  planar.  The  wave  propagation  is  created  outward  by  the  oscillating  surface.  To  eliminate  the  complexity  of  the 
counter-reflected  wave  effects,  the  container  size  must  be  large  enough  so  that  no  disturbance  reaches  the  walls. 
The  viscosity  can  be  found  from  the  relationship  between  the  frequency  of  the  oscillation  and  the  shear  rate. 

The  viscosity  measurement  in  the  current  study  was  measured  by  this  kind  of  viscometer. 


Summary  nf  Experimental  Investigations  of  Liquid  Viscosity  of  Refrigerants 


Table  4.8  is  a  summary  of  previous  experimental  investigations  for  the  viscosity  of  liquid  refrigerants  [65], 


TABLE  4.8:  LIST  OF  PREVIOUS  INVESTIGATIONS  OF  LIQUID 
REFRIGERANT  VISCOSITY  MEASUREMENTS 


v^aveiM^itov":' 

mmmm, 

B5m 

Kefrifi. 

im 

SHQI 

Benning  and  1 

1939 

CFC-11 

243-333 

sat. 

roll,  ball 

- 

[68] 

Markwood 

CFC-12 

243-332 

sat. 

roll,  ball 

* 

HCFC-22 

240-318 

sat. 

roll,  ball 

CFC-113 

243-333 

sat. 

roll,  ball 

CFC-114 

242-332 

sat. 

roll,  ball 

Lilios 

1957 

243-332 

sat. 

capillary 

2.2 

[69] 

1 

213-243 

sat. 

roll,  ball 

203?218 

sat. 

roll,  ball 

mwia;. 

243-293 

sat. 

roll,  ball 

Gordon  et  al. 

1969 

CFC-11 

245-347 

sat. 

capillary 

■gn 

CFC-12 

245-305 

sat. 

capillary 

HCFC-22 

246-312 

sat. 

capillary 

R-13bl 

246-301 

sat. 

capillary 

Phillips  and 

1970 

Viirssm 

209-352 

sat. 

capillary 

[64] 

Murphy 

IfiisalSi 

202-312 

sat. 

capillary 

HCFC-22 

201-299 

sat. 

capillary 

R-152a 

200-318 

sat. 

capillary 

- 

Kumagai  and 

1990 

11  kinds 

273-353 

sat. 

capillary 

±9 

[65} 

Takahashi 

of  refrig. 

sat. 

-13 

■MBS 

for  HFC-134a 

1992 

CFC-123 

233-418 

20  MPa 

Egggl 

[10] 

HFC-134a 

213-423 

30  MPa 

Assael  et  al. 

1992 

msmm 

270-370 

300  MPa 

vibrating 

±3 

Diller  et  al. 

1993 

HFC-134a 

175-320 

sat. 

torsional 

±3-8 

[13] 

crystal 

compared  to  others 

CFC-123 

170-320 

sat. 

torsional 

±7-12 

crystal 

compared  to  others 

R-141b 

175-320 

sat. 

torsional 

less  than  ±12 

crystal 

compared  to  others 

Note:  accuracy  not  reported,  "sat."  saturation  state 
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CHAPTER  5 


EXPERIMENTAL  APPARATUS  AND  OPERATION 

This  chapter  describes  the  details  and  designs  of  the  experimental  facility.  The  construction  of  the  test  rig 
is  described  and  shown.  Sensors  and  test  section  calibration  are  also  presented.  Finally,  the  experimental 
procedures  are  discussed  for  proper  operation  of  this  test  facility. 


EXPERIMENTAL  FACILITY  CONSTRUCTION 

The  experimental  test  rig  was  established  for  the  purpose  of  measuring  liquid  transport  properties  of 
refrigerant,  refrigerant  mixtures,  or  refrigerant/lubricant  mixtures.  Each  experiment  measures  the  single-phase 
heat  transfer  coefficients  of  liquids  by  a  surface  temperature  method. 

Heat  Transfer  Test  Facility 

The  test  section  is  a  0.9525-cm  (3/8-inch)  inside  diameter  by  2-m  long  smooth  copper  tube.  The  measured 
quantities  are  tube  wall  temperature,  inlet/outlet  fluid  temperatures,  absolute  and  differential  pressures,  dynamic  viscosity, 
and  the  mass  flow  rate.  A  schematic  of  this  test  rig  is  shown  in  Figure  5.1,  while  a  photograph  is  shown  in  Figure  5.2. 

Eleven  T-type  thermocouples  are  installed  on  the  outer  wall  of  the  tube  at  equal  distances  of  0.2  m, 
starting  from  the  inlet  point  and  ending  at  the  outlet  point,  along  the  2-m  long  test  section.  In  order  to  get  more 
average  temperature  measurements  at  the  inlet  and  outlet  locations,  two  additional  thermocouples  were  placed  0.1 
m  from  the  inlet  and  the  outlet  points  on  the  outer  tube  wall.  Moreover,  one  thermocouple  is  placed  on  the  outer 
insulated  wall  surface  for  measuring  the  temperature  there.  Also,  one  thermocouple  is  used  for  measuring  the 
room  temperature.  The  thermocouple  for  measuring  room  temperature  is  located  in  a  wooden  box  which  is  wrapped 
with  aluminum  foil  on  the  outer  surfaces  to  insure  the  measurement  of  a  radiation,  and  convection-free 
temperature.  Figure  5.3  shows  the  details  of  the  test  section  configuration.  A  picture  of  the  test  section  is  also 
shown  in  Figure  5.4.  Two  RTDs  were  placed  right  at  the  inlet  and  outlet  points  of  the  test  section  to  measure  the 
respective  fluid  temperatures.  All  thermocouples  and  RTDs  were  calibrated  to  ±0.05°C. 
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Figure  5.1:  Schematic  diagram  of  test  facility. 
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Figure  5.2:  Photograph  of  test  facility. 
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igure 
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Figure  5.4:  Photograph  of  test  section. 


In  order  to  get  more  accurate  inlet  and  outlet  temperatures,  three  plastic  mixers  were  installed  inside  the 
tube  for  the  purpose  of  determining  a  more  accurate  average  fluid  temperature.  The  first  one  was  located  vertically 
just  prior  to  the  tube  elbow  of  the  inlet  portion,  the  second  one  was  horizontally  placed  between  the  test  section  and 
the  tube  elbow,  and  the  third  one  was  installed  horizontally  immediately  after  the  outlet  of  the  test  section. 

Two  plastic  tube  connectors  which  served  as  heat  insulators  were  placed  at  both  the  inlet  and  outlet  of  the 
test  section  in  order  to  prevent  axial  heat  loss  from  the  ends  of  the  test-tube.  One  was  placed  right  before  the  inlet 
point  of  the  test-tube,  while  the  other  was  placed  just  after  the  outlet  point  of  the  test-tube.  The  test  section 
pressure  was  monitored  by  an  absolute  pressure  gauge  and  a  differential  pressure  transmitter  which  measured  the 
test  section  inlet  and  differential  pressures  across  the  test-tube,  respectively. 

The  system  pressure  was  controlled  by  an  accumulator  which  was  connected  to  a  nitrogen  vessel. 

Pressure,  therefore,  was  regulated  by  the  nitrogen  pressure.  The  system  mass  flow  rate  was  adjusted  by  a  variable- 
speed  controlled  motor  pump  assembly.  The  system  heating  control  was  achieved  by  two  parallel  power  connected 
heat  tapes  of  675  watts  each.  The  heating  capacity  was  controlled  by  adjusting  the  power  regulator  which  was 
connected  to  a  watt  transducer  with  an  accuracy  of  ±0.2%  of  reading.  The  heat  tapes  were  uniformly  wrapped 
along  the  test  section;  thus,  constant  heat  flux  was  expected  and  assumed. 

After  the  heated  test  section,  two  condensers  were  placed  between  the  outlet  of  the  test  section  and  the 
inlet  of  the  pumps  for  removing  the  heat  that  was  added  in  the  test  section.  Two  pumps  were  installed  in  a  parallel 
arrangement  with  each  of  the  condensers  in  order  to  enable  a  wide  operating  range  of  mass  flow  rates.  The  cooling 
capacity  could  also  be  controlled  by  adjusting  the  chiller  water  flow  rate  as  shown  in  Figure  5.1. 

All  of  the  measurement  signals  were  connected  to  a  HP3457 A  switch  and  control  unit  and  a  HP3488 A 
multimeter.  Data  acquisition  and  measuring  process  control  were  driven  by  a  personal  computer. 

Viscometer 


Viscosity  is  simultaneously  measured  with  other  measurements  in  the  test  system.  In  this  study,  the 
viscosity  is  measured  by  an  inline  viscometer  installed  at  the  inlet  portion  of  the  test  section.  Viscosity  can  be 
dynamically  measured  at  the  same  time  along  with  the  other  measurements.  This  viscometer  uses  the  torsional 
oscillation  principle  which  is  the  principle  of  "surface  load."  A  vibrating  surface  is  in  contact  with  a  liquid  which 
experiences  a  force  that  is  a  function  of  the  viscosity  [14].  This  type  of  viscometer  measures  the  product  of  dynamic 
viscosity  and  density.  The  range  of  this  instrument  can  vary  from  a  low  viscosity  range  of  0.1  cp  •  g  /  cm3  to  a 
high  viscosity  range  of  500  cp  •  g  /  cm3,  which  covers  most  of  the  pure  refrigerants,  refrigerant  mixtures,  and  even 
refrigerant/  lubricant  mixtures  with  low  lubricant  concentrations.  The  accuracy  of  the  viscosity  measurement  is 
±2%  of  reading.  The  fluid  is  controlled  by  a  bypass  valve  which  regulates  flow  through  a  sample  cell  in  which  the 
sensor  of  the  viscometer  is  mounted.  The  temperature  of  the  fluid  at  the  sample  cell  is  also  measured  by  a  calibrated 
RTD  with  an  accuracy  of  ±0.05°C.  To  determine  viscosity,  density  is  required.  In  this  test  rig,  density  is  directly 
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measured  by  a  densimeter  of  the  vibrating-wire  type  which  is  installed  in  series  with  the  viscometer.  A  schematic 
diagram  is  shown  in  Figure  5.5.  A  photograph  of  the  viscometer  sample  cell  is  shown  in  Figure  5.6. 

SENSOR  CALIBRATIONS 

The  physical  quantities  measured  directly  by  this  test  rig  include  temperature,  mass  flow  rate,  viscosity, 
pressure,  and  power  input  The  temperature  measurements  include  inlet  and  outlet  temperatures,  test-tube  wall 
temperatures,  insulated  wall  temperature,  and  the  ambient  temperature.  The  inlet  and  outlet  temperatures  are 
measured  by  RTD  sensors  while  the  other  temperatures  are  determined  with  thermocouples.  All  the  RTD  and 
thermocouple  sensors  were  calibrated  by  a  standard  thermometer  with  an  accuracy  of  ±0.05 °C. 

The  other  measurement  devices  such  as  for  mass  flow  rate,  viscosity,  pressure,  and  power  were  calibrated 
at  the  factory.  The  uncertainty  of  each  sensor  is  listed  in  Table  5.1. 


TABLE  5.1:  SENSOR  UNCERTAINTY 


measured  auantity 

1  uncertainty 

mass  flow  rate 

warn 

±0.15%  of  reading 

pressure 

psi 

±0.2%  of  reading 

viscosity  •  density 

cp  *  glen ? 

±2%  of  reading 

| 

density 

±0.15%  of  reading 

power 

w 

±0.2%  of  reading 

temperature  (RTD,  thermocouples) 

±0.05°C 

TEST  SECTION  HEAT  LOSS  CALIBRATION 

The  energy  balance  principle  can  not  only  be  applied  to  estimate  the  beat  loss  of  the  test  section  for  fluids 
with  known  specific  heats,  but  it  can  also  be  applied  in  reverse  to  obtain  the  heat  loss  for  fluids  of  unknown  specific 
hrats  However,  to  get  accurate  specific  heat  measurements,  heat  loss  calibrations  using  several  fluids  with  known 
specific  heats  are  required.  It  should  be  noted  that  the  heat  loss  estimation  will  be  dependent  on  the  specific  test 
section  with  its  own  geometric  and  insulating  characteristics.  In  other  words,  the  heat  loss  is  only  valid  for  the  test 
section  being  used  for  calibration  and  cannot  be  applied  to  other  test  sections  or  any  situations  changed  from  the  original. 

The  heat  loss  of  the  test  section  was  obtained  from  the  calibration  of  previously  tested  refrigerants  which 
areHCFC-22,  CFC-12,  and  CFC-113.  The  accuracy  of  the  heat  loss  calibration  was  tested  using  CFC-1 14.  Then,  the 
net  beat  input  of  the  test  section  was  used  to  calculate  the  specific  heat  of  the  test  fluids.  Figure  5.7  shows  the 
logarithmic  function  of  the  test  section  heat  loss  associated  with  the  logarithmic  function  of  the  temperature 
diffi»n»nn»  between  the  average  outer  insulated  surface  temperature,  Tv  and  average  ambient  temperature,  Tg. 
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Figure  5.5:  Viscometer  construction. 
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Figure  5.6:  Photograph  of  viscometer  sample  cell 


The  curve  fit  equation  of  the  heat  loss,  Qhss  (in  watts),  associated  with  the  In  (T}  -  Ta)  is  shown  below: 

In  (£to„)  -  0.612  +  1.33  In  (Ts  -  Ta)  (5.1) 

It  should  be  noted  that  this  equation  will  be  only  applicable  to  this  particular  test  section  as  constructed 
including  the  particular  insulation.  For  any  other  construction,  the  heat  loss  must  be  redetermined  and  the  heat 
loss  equation  must  also  be  reformulated. 

ACCURACY  VERIFICATION  OF  DENSIMETER 

Although  the  accuracy  of  the  densimeter  was  claimed  to  be  ±0.15%  within  1300  kg/m3,  it  was  necessary  to 
verify  the  accuracy  with  refrigerants  of  known  density.  The  verified  refrigerants  used  were  HCFC-22,  CFC-12, 
CFC-1 13,  and  CFC-1 14,  which  were  also  the  refrigerants  for  properties  measured  by  using  the  current 
methodology  and  facility.  Figure  5.8  shows  the  plot  of  measured  density  versus  the  ASHRAE  reported  density.  As 
indicated  in  the  deviation  band  in  this  figure,  the  deviation  was  within  ±2%.  However,  it  is  necessary  to  mention 
that  some  of  the  measured  densities  were  higher  than  1300  kg/m3,  which  is  considered  the  upper  level  with  highest 
accuracy  for  the  densimeter  used.  Moreover,  all  the  densities  measured  were  under  a  compressed  liquid  state 
rather  than  a  saturation  state.  As  it  can  be  seen,  most  densities  were  measured  within  ±1%. 

ACCURACY  VERIFICATION  OF  VISCOMETER 

Viscosity  is  measured  by  an  inline  torsional  oscillation  viscometer.  This  type  of  instrument  measures  the 
product  of  kinematic  viscosity  and  density.  The  accuracy  was  verified  using  CFC-12,  CFC-1 14,  CFC-1 13,  and 
pure  water.  Figure  5.9  shows  the  plot  of  the  measured  viscosity  of  CFC-12  versus  temperature  compared  with 
ASHRAE  data.  The  deviation  percentage  plot  is  shown  in  Figure  5.10.  As  the  plots  show,  the  measured  viscosity 
matches  the  ASHRAE  data  within  ±2%.  The  measured  viscosity  of  CFC-1 14  and  its  deviation  compared  with  the 
ASHRAE  data  are  shown  in  Figures  5.11  and  5.12,  respectively,  while  the  measured  viscosity  of  CFC-1 13  and  its 
deviation  compared  with  the  ASHRAE  data  are  shown  in  Figures  5.13  and  5.14,  respectively.  The  repeatability  of 
the  measurements  was  tested  by  CFC-1 14,  which  began  with  the  temperature  around  room  temperature  and  went 
down  to  nearly  0°C.  The  temperature  then  went  up  to  nearly  50°C  and  back  to  room  temperature.  Finally,  the 
measured  repeatability  of  viscosity  was  verified  with  pure  water.  This  test  measured  the  viscosity  of  pure  water  at 
room  temperature  for  nearly  17  hours.  The  measured  viscosity  of  pure  water  and  its  deviation  with  the  Chemistry/ 
Physics  Handbook  (CPH)  data  [70]  are  shown  in  Figures  5.15  and  5.16,  respectively.  As  can  be  seen,  the  measured 
accuracy  was  within  ±2%. 
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Figure  5.7:  Heat  loss  estimation  as  a  function  of  wall-fluid  temperature  difference. 
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Temperature,  T,  (°C) 

Figure  5.10:  Deviation  of  CFC-12  measured  viscosity  compared  with  ASHRAE  data. 


Legend 
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Figure  5.12:  Deviation  of  measured  and  ASHRAE-listed  viscosity  for  CFC-114. 
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Figure  5.13:  CFC-1 13  measured  viscosity  as  a  function  of  temperature. 
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Figure  5.14:  Deviation  of  measured  and  ASHRAE-listed  viscosity  for  CFC-1 1 3. 


Legend 
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Figure  5.16:  Pure  water  measured  viscosity  deviation  for  stability  test 


EXPERIMENTAL  OPERATION 


This  section  describes  the  operation  of  the  test  rig  and  the  experimental  procedures.  In  order  to  obtain 
reliable  test  data,  the  test  system  should  be  operated  under  specific  procedures.  Based  upon  operational 
experiences,  these  procedures  must  be  followed  by  the  operator  in  order  to  achieve  success.  The  test  procedures  are 
outlined  in  the  following  discussion. 

Test  System  Preparation 

Before  charging  the  facility  with  refrigerant,  the  test  system  should  be  evacuated.  A  system  pressure  less 
than  2.76  kPa  is  considered  to  be  evacuated  enough  and  ready  for  charging  refrigerant.  The  duller  can  be  turned  on 
and  the  chiller  working  fluid  can  be  circulated  through  the  heat  exchangers  while  the  pressure  approaches  2.76  kPa. 
Meanwhile,  the  refrigerant  to  be  used  can  be  pre-heated  in  order  to  increase  the  pressure  which  eases  transfer  into 
the  system.  During  this  time,  the  data  acquisition  system  is  monitoring  the  status  of  the  system. 

Some  potential  problems  can  occur  if  proper  preparations  are  not  made.  These  are: 

*  Check  the  material  compatibility  of  refrigerant  test  system  materials  with  refrigerants  in  all  wet  parts 
(seals).  Elastomer  material,  such  as  pressure  relief  valve  seals,  accumulator  bladder,  and  pump 
diaphragm  could  be  eroded  by  refrigerants  due  to  incompatibility. 

*  Leakage  of  the  system  could  cause  a  loss  of  refrigerant  or  an  intake  of  atmospheric  air.  This  could 
cause  problems  if  non-condensable  gases  were  to  exist  in  the  system,  and  an  adequate  system 
evacuation  could  be  difficult  to  achieve.  Leaks  should  be  detected  and  fixed  for  losses  of  14  kPa  (or 
more)  in  a  24-hour  period. 

*  In  order  to  store  more  refrigerant  in  the  accumulator,  chiller  glycol  must  be  circulated  around  the 
accumulator  in  order  to  keep  the  refrigerant  cool  and  ease  charging  the  accumulator.  In  addition,  the 
bladder  of  the  accumulator  should  be  pre-charged  to  a  pressure  up  to  276  kPa  in  order  to  control  the 
space  for  receiving  the  refrigerant  Failure  to  do  this  will  result  in  improper  functioning  of  the  accumulator. 

*  Check  all  the  valves  to  make  sure  they  are  appropriately  opened  or  closed. 

*  Monitor  the  system  status  to  see  if  its  properties  are  in  the  desired  range  (pressures  and  temperatures). 

*  Ensure  that  a  refrigerant  cylinder  is  used  as  a  buffer  for  overcharged  refrigerant,  this  is  required  by  the 
facility  for  reclaiming  extra  refrigerant  due  to  density  changes  (temperature  fluctuations).  This 
cylinder  should  be  properly  evacuated  when  hooked  up  to  the  system. 
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Viscometer  Preparation 


The  type  of  viscometer  used  measures  the  oscillation  frequent  of  a  cylinder  sensor  which  is  immersed  in 
the  refrigerant  in  a  sample  cell.  There  are  eight  bolts  surrounding  the  flange  to  connect  the  transmitter  portion  and 
the  sample  cell.  The  torque  balance  of  the  bolts  is  very  sensitive  to  the  frequency  output  which  is  converted  to  a  digital 
cunent  output  to  the  DAS  (data  acquisition  system).  The  following  procedures  are  required  for  viscometer  preparation: 

•  Carefully  lift  up  the  viscometer  sensor  and  carefully  clean  the  sensor  surface  with  a  soft  tissue  to  make 

sure  no  impurities  are  deposited  on  the  surface. 

•  Ch^rV  the  gasket  material  and  inspect  for  any  damage.  Replace  it  if  it  is  necessary. 

•  Carefully  place  the  sensor  into  the  sample  cell  while  lining  up  all  the  bolts. 

•  Tighten  each  of  the  eight  bolts  to  a  torque  of  45  Ibf-  ft . 

•  Eliminate  all  possible  vibration  sources  in  order  to  avoid  any  effects  on  the  measured  frequency. 

•  Adjust  the  potentiometer  of  the  viscometer  console  to  zero  voltage  output  at  a  current  at  4  mA  when 
the  sensor  is  subjected  to  complete  evacuation  (no  material).  If  this  is  not  properly  done,  an  inaccuracy 
could  be  introduced  in  the  measurement. 

Charring  Refrigerants 

The  following  steps  outline  the  procedures  for  charging  refrigerants: 

•  Make  sure  the  system  is  evacuated  at  a  proper  pressure  and  the  chiller  glycol-water  mixture  is 

circulating  in  the  heat  exchangers. 

•  Slowly  open  the  valve  and  let  refrigerant  flow  into  the  system. 

•  Watch  the  sight-glass  to  make  sure  refrigerant  is  flowing  into  the  system. 

•  Make  sure  the  liquid  refrigerant  fills  the  housings  of  the  refrigerant  pumps.  Then,  turn  on  the  pumps 
and  operate  at  a  low  speed. 

•  Release  the  nitrogen  gas  inside  the  accumulator  to  allow  proper  space  for  the  storage  of  liquid 
refrigerant. 

•  Continue  charging  refrigerant  until  the  refrigerant  is  seen  to  fully  fill  the  plastic  tubes  at  the  highest 
point  of  the  system.  Then  close  the  refrigerant  valve. 

•  When  the  refrigerant  is  properly  charged,  glycol  circulating  around  the  accumulator  could  be 

Charge  nitrogen  pressure  into  the  accumulator  to  the  desired  pressure. 

•  When  charging  a  non-azeotropic  refrigerant  mixture,  because  of  the  influence  of  composition  on  the 
mixture's  properties,  a  liquid  mixture  must  be  charged  to  the  test  system.  Vapor  charging  should  be 
avoided.  In  addition,  to  prevent  the  composition  from  changing,  the  use  of  distilled  refrigerants  should 
also  be  avoided. 
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System  Operation  and  Data  Acquisition 


Data  arc  taken  under  steady-state  conditions.  The  following  steps  control  steady-state  operation: 

•  Apply  power  to  the  heat  tape  and  adjust  power  to  desired  wattage  by  the  Variac. 

•  Apply  pressure  by  adjusting  the  nitrogen  regulation  valve  of  the  supply  tank. 

•  Adjust  the  pump  speed  controller  to  control  the  refrigerant  mass  flow  rate. 

•  Adjust  the  controlling  valve  of  the  glycol-waier  mixture  to  a  desired  flow  rate  as  indicated  by  the  rotameter. 

•  Steady-state  in  the  test  section  has  been  achieved  when  the  heat  input  equals  the  heat  removed.  Repeat 
any  of  the  above  steps  as  required  for  reaching  steady-state. 

•  Steady-state  is  considered  to  be  reached  by  monitoring  the  inlet  and  outlet  temperature  change  of  the 
test  section  within  ±0.05°C  between  two  scans  (approximately  one  minute).  Then,  a  data  point  is 
determined  as  the  average  value  of  20  consecutive  scans. 

•  Another  test  point  may  be  taken  by  repeating  the  above  procedures  for  a  different  inlet  temperature, 
mass  flow  rate  of  refrigerant,  or  mass  flow  rate  of  chiller  glycol. 

Injecting  Lubricant  and  Sampling  Lubricant  Concentration 

Lubricant  concentration  is  an  important  quantity  for  understanding  the  lubricant’s  effects  on  properties. 
There  are  a  number  of  methods  which  could  have  been  adopted  in  the  current  study,  such  as  the  light  absorption 
method  [71],  vibrating  U-tube  densimeter  sensor  and  ultrasonic  acoustic  sensor  [72,73],  viscometer  lubricant 
concentration  sensor  [74],  and  sampling  mass  method.  In  the  current  study,  the  last  method  was  adopted  over  the 
other  methods  due  to  their  limited  availability. 

The  injection  of  lubricant  into  the  test  system  is  achieved  by  an  oil-injection  piston  with  a  15.24-cm 
stroke.  The  lubricant  concentration  sample  cell  is  in  line  with  the  injection  cylinder.  This  assembly  is  placed  at 
the  inlet  portion  of  the  test  section.  The  schematic  diagram  is  shown  in  Figure  5.17.  A  picture  of  the  oil-injection 
assembly  and  oil-sample  device  is  shown  in  Figure  5. 18.  The  operational  procedures  are  described  as  follows. 

Oil  Injection— 

•  Evacuate  the  dead  space  of  the  piston  cylinder  assembly  and  fully  close  the  valve  at  the  oil  outlet  after 
completing  evacuatioa 

•  Evacuate  the  upper  volume  of  the  cylinder  above  the  piston  and  fully  close  the  upper  valve  after 
finishing  the  evacuation. 

•  Suck  oil  into  the  cylinder  from  the  lower  oil  outlet  by  slowly  opening  the  needle  valve  while  preventing 
air  from  leaking  into  the  cylinder. 

•  Close  the  lower  needle  valve  when  finished  charging  the  oil. 


52 


evacuating 

connector 


53 


Figure  5.17:  Schematic  diagram  of  oil  injection  and  sampling  devices. 
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•  Place  the  cylinder  up-side-down  and  evacuate  gases  from  the  cylinder  top  by  hooking  up  an  evacuating 
hose  to  remove  the  residual  air  remaining  or  dissolved  in  the  oil. 

•  Install  the  cylinder  at  the  position  for  oil  injection. 

•  Open  the  valves  connected  to  the  test  system  and  apply  high  pressure  nitrogen  at  the  top  of  the  piston. 

•  Inject  oil  into  the  system  until  the  desired  amount  of  oil  is  injected. 

•  Circulate  refrigerant  while  the  oil  is  being  injected  and  continuously  circulate  for  at  least  10  hours  to 
insure  the  oil  is  uniformly  mixed  with  the  refrigerant.  Finally,  sample  the  lubricant  concentration  in 
the  refrigerant. 

Lubricant  Concentration  Sampling- 

•  Close  the  valves  of  the  sample  cell  to  isolate  it  from  the  system. 

•  Remove  the  sample  cell  and  weigh  the  total  mass  of  the  cell  including  the  refrigerant  and  lubricant. 

•  Slowly  open  the  top  valve  of  the  sample  cell  to  allow  the  refrigerant  to  evaporate  from  the  cell.  An  oil 
absorber  hooked  up  to  the  cell  can  be  used  when  the  vacuum  pump  is  operated. 

•  Weigh  the  cell  and  oil  inside  the  cell  and  record  until  the  weight  remains  constant. 

•  Calculate  the  oil  concentration  by  dividing  the  net  oil  weight  (oil  in  the  cell  +  oil  in  the  oil-absorber) 
by  the  total  mixture  weight  (refrigerant  +  oil). 
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CHAPTER  6 


THEORY  OF  THERMAL  CONDUCTIVITY  MEASUREMENT 

In  this  chapter,  the  theory  of  thermal  conductivity  measurement  will  be  described.  Two  methods  were 
developed  in  this  study,  which  will  be  referred  to  as  Approach  1:  The  Nusselt  number  method,  and  Approach  2: 
The  Prandtl  number  method.  In  Approach  1,  the  heat  transfer  characteristics  of  the  test  section  and  Nusselt 
number  correlations  were  used  to  obtain  the  thermal  conductivity,  while  in  Approach  2,  the  thermal  conductivity 
was  directly  related  to  the  Reynolds  number,  ReD,  and  the  temperature  characteristics  of  the  test-tube.  The  detailed 
mathematical  approaches  are  described  below. 


APPROACH  1:  THE  NUSSELT  NUMBER  METHOD 

In  this  method,  the  thermal  conductivity  is  obtained  from  heat  transfer  characteristics  of  a  heated  test- 
tube.  From  experimental  measurements,  an  average  heat  transfer  coefficient  for  single-phase  conditions  in  the 
test-tube  can  be  determined.  Meanwhile,  the  heat  loss  from  the  test-tube  can  also  be  measured  for  a  specific  test- 
tube.  The  average  Nusselt  number  was  calculated  from  the  heat  transfer  measurements,  and  the  thermal  conductivity 
was  then  calculated  from  the  Nusselt  number.  Details  of  the  measuring  principles  are  described  as  follows. 

Heat  Transfer  Measurement 

The  average  heat  flux  is  calculated  from  the  following  equation: 


q"  =  h(Tw-Tf) 


(6.1) 


Thus,  the  average  heat  transfer  coefficient  can  be  written  as: 


Ti  = 


(Tw~Tf) 


(6.2) 


where  Tw  and  Tj  are  the  average  tube  wall  temperature  and  fluid  temperature,  respectively.  The  quantity  Tw-Tf 
can  be  calculated  from  either  the  individual  average  wall  temperature,  Tw ,  and  fluid  temperature,  Tf,  or  the 
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quantity  Tw  -  Tf,  which  is  described  in  a  later  section.  The  net  heat  transfer  rate,  qnet ,  to  the  fluid  can  be 
calculated  from  the  total  heat  input  from  the  power  applied,  'qM ,  and  heat  loss  to  the  environment,  q^ ,  as  follows. 


Qnet 


Qtot  Qloss 


(6.3) 


The  purpose  of  measuring  heat  loss  from  the  test  section,  is  to  obtain  the  net  heat  input  to  the  test 
section,  qner  The  heat  loss  is  transferred  to  the  environment  around  the  test  section  by  a  combination  of  natural 
convection  and  radiation.  For  calculating  the  heat  loss,  the  following  equation  is  used: 


where  Ts  and  Ta  are  the  average  outer  insulation  surface  temperature  over  area,  A0 ,  and  the  ambient  room  temperature, 
respectively,  which  can  be  obtained  from  experimental  measurements.  The  average  heat  transfer  coefficient  for  the 
heat  loss  of  the  test  section,  h0,  can  be  determined  and  calibrated  from  experimental  measurements  with  a  fluid  of 
known  specific  heat,  Cp.  This  equation  has  already  been  shown  in  Chapter  5  as  Equation  5.1.  Once  the  heat  loss 
characteristics  of  the  test  section  are  obtained,  they  can  be  applied  to  the  Cp  calculation  for  the  other  unknown  fluids. 

By  knowing  the  test  section  heat  loss  characteristics,  the  qneP  can  be  calculated  from  an  energy  balance  for 
a  working  fluid.  Then,  Cp  for  an  unknown  fluid  can  be  obtained  as  follows: 


Qnet 

Cp  " 


(6.5) 


where  Tf  and  T0  are  the  mean  inlet  and  outlet  fluid  temperatures,  respectively,  while  m  is  the  mass  flow  rate.  This 
specific  heat  determination  is  important  for  new  refrigerant  mixtures  and  blends  and  for  refrigerant/lubricant 
mixtures  when  non-ideal  mixing  occurs.  In  both  cases,  a  theoretical  calculation  of  Cp  for  the  pure  fluids  (either 
the  pure  refrigerant  or  pure  lubricant)  will  not  result  in  an  accurate  specific  heat  value. 


Nusselt  Number  Correlations 


Heat  transfer  data  are  usually  given  by  the  Nusselt  number  which  is  defined  a shD/k.  A  number  of 
single-phase  Nusselt  number  correlations  have  been  published  for  convection  heat  transfer.  Some  examples  are  the 
Dittus-Boelter  correlation  [75,76],  the  Petukhov  and  Popov  correlation  [77],  and  the  Gnielinski  correlation  [78], 
These  three  correlations  are  given  below: 


The  Dittus-Boelter  correlation  (for  heating): 


Nud  =  0.023 ReD0*Pr°4 


(6.6) 
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The  Petukhov  correlation: 


m)  RenPr _ 

NuD  =  1.07  +  12.7  (ft  8)1/2  (Pr273  -  1) 


(6.7) 


The  Gnielinski  correlation: 


-  078)  (Ren  -  1000)  Pr 

NuDa  1  +  12.7  (f!  8)172  (Pr273  -  1) 


(6.8) 


where  f  is  friction  factor.  For  smooth  tubes,  expressions  for/ are: 


Petukhov  correlation, 

/  =  (l.Z2logl0ReD  -  1.64)*2  (6.9) 

Gnielinski  correlation. 


/=  (0.19lnReD  -  1.64)*2 
In  the  above  equations,  ReD  and  Pr  are  defined  as: 


(6.10) 


(6.11) 

(6.12) 


As  mentioned  previously,  with  h  and  m  known,  the  thermal  conductivity  can  be  obtained  from  the  equality  by 
setting  the  Nusselt  number,  hD/k,  equal  to  the  above  correlations.  Equations  6.6, 6.7,  or  6.8. 


Calibration  Function 

The  method  used  to  determine  thermal  conductivity  involves  applying  the  Nusselt  number  correlations 
and  solving  for  the  thermal  conductivity.  Therefore,  an  accurate  Nusselt  correlation  is  necessary  in  order  to  obtain 
an  accurate  value  of  thermal  conductivity.  Although  some  published  Nusselt  correlations  have  claimed  good 
accuracy,  their  calibrations  from  the  heat  transfer  data  are  still  needed.  This  is  because  the  correlations  were 
obtained  by  curve  fitting  a  large  number  of  data  with  a  wide  range  of  Pr  and  Re^  which  might  not  apply  to  some 
specific  fluids  with  a  limited  range  of  Pr  and  ReD.  Moreover,  the  thermal  conductivity  was  originally  used  for 
obtaining  those  correlations.  Now  it  will  be  obtained  in  a  reversed  manner.  Therefore,  in  order  to  obtain  an 
acceptable  thermal  conductivity,  an  accurate  Nusselt  correlation  is  required  in  the  specific  range  of  Pr  and  ReD. 


58 


The  method  used  for  calibrating  Nusselt  number  correlations  assumes  a  calibrated  function  exists  between 
the  values  of  Nusselt  number  calculated  from  the  correlations  and  definition.  The  calibration  function  will  be  denoted 
as  CF.  With  this  assumption,  a  set  of  CF  values  can  be  obtained  from  experimental  data.  Then,  the  CF  correlations 
can  be  obtained  by  a  curve  fit  for  various  Pr  and  ReD.  The  mathematical  procedures  are  interpreted  as  follows. 


The  definition  of  the  average  Nusselt  number  is: 


The  CF,  therefore,  exists  between  NuD  and  correlations  like 

Nud  " 

CF  -  -  FVr*4 

D correlation 

The  calibrated  Dittus-Boelter  correlation  (for  heating): 


Nud  =  0.023 ReJ'P^CF 


The  calibrated  Petukhov  correlation: 


Nud 


(f/8)  RepPr 

1.07  +  12.7  (ft  8),/2  (Pr273  -  1) 


The  calibrated  Gnielinski  correlation: 


-  Q78)  (Ren  -  1000)  Pr 

N*l>m  1  +  12.7  ( ft  8)1/J  (Pr273  -  1) CF 


(6.13) 


(6.14) 


(6.15) 


(6.16) 


(6.17) 


Regression  of  the  Calibration  Function 


Because  the  ReD  and  Pr  are  not  of  the  same  order  of  magnitude,  a  logarithmic  transformation  function  is 
applied  for  convenience.  Other  types  of  transformations  are  possible.  In  the  following  equations,  the  superscript, 
e,  denotes  the  curve  fit  equations  obtained  from  measured  data  while  the  superscript,  c,  denotes  curve  fit  equations 
obtained  from  correlations.  Thus,  Equation  6. 14  becomes  in  natural  logarithmic  form: 


In  (Nu£)  -  In  Nucd  +  In  (CF) 

(6.18) 

In  (CF)  =  ln(Nu^-lnNucD 

(6.19) 
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A  curve  fit  for  In  ( CF )  can  then  be  obtained.  A  linear  relationship  between  In  (NuD  /  Pr")  and  In  (ReD)  is 
well  known  from  past  knowledge.  Therefore,  an  easy  method  to  correlate  these  two  logarithmic  quantities  is  of  the 
linear  form: 

In  (Nud  /  PH1)  -  A+B  In  (ReD)  (6.20) 


where  A,  B,  and  n  are  the  curve  fit  constants.  Different  constants  were  obtained  from  different  applied 
correlations.  The  calibration  function,  CF,  is  therefore  different  for  the  various  correlations  and  is  denoted  as: 


Nu*d 

CF  =  '=~~F 
Nucd 


(6.21) 


Once  the  calibration  function,  CF,  is  obtained  as  a  function  of  ReD  and  Pr,  the  modified  Nusselt 
correlations  can  be  used  to  determine  thermal  conductivity  as  shown  in  the  following  sections.  With  this 
assumption,  a  set  of  CF  values  can  be  obtained  from  experimental  data  and  correlation  calculations  for  different 
fluids.  CF  correlations  can  then  be  obtained  by  curve  fitting  CF  with  Pr  and  ReD.  The  mathematical  procedures 
are  described  below. 

For  heat  transfer  in  tubes  in  this  study,  the  following  curve  fit  equations  are  used  for  experimentally 
measured  Nusselt  numbers,  Nu^,  and  correlation  Nusselt  numbers,  M/£,  respectively. 

Nu*d  =  CeRe£Prl>e  (6.22) 

Nucd  -  CRetfr*  (6.23) 

It  should  be  noted  that  the  Dittus-Boelter  equation  is  already  in  this  form  while  the  other  two  in-tube 
equations  (Equations.  6.7  and  6.8)  must  be  curve  fit  into  this  form  for  consistent  treatment.  In  addition,  the 
magnitudes  of  Pr  and  Re  are  quite  different.  Again,  a  logarithmic  function  is  recommended  for  convenient  curve 
fitting.  By  dividing  Equation  6.22  by  Equation  6.23,  CF  can  be  expressed  as: 


CF  = 


C* 

C 


(6.24) 


Thermal  Conductivity  Calculations 

After  a  calibrated  is  obtained,  the  heat  transfer  coefficient  can  be  calculated  from  an  average  (fw  -  Tj), 

which  is  measured.  The  Nusselt  number  will  be  obtained  thereafter.  Then,  the  thermal  conductivity  can  be 
obtained  from  the  calibrated  Nusselt  number  correlations.  However,  the  calibration  function,  CF,  should  be 
obtained  from  known  property  fluids.  Once  the  calibration  function,  CF,  is  obtained  as  a  function  of  ReD  and  Pr, 
the  thermal  conductivity  can  be  expressed  in  general  form  as  follows. 
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k  -  [cF<MM(fyww»>r,pi> 

where  he  denotes  the  average  heat  transfer  coefficient  from  experimental  measurements. 

For  each  calibrated  Nusselt  number  correlation,  the  k  values  are  expressed  as  the  following: 


(6.25) 


for  Dittus-Boelter  correlation: 

kd  =  [  0.023  ReDMPi°ACF]‘l  hD 


for  Petukhov  correlation: 

_  1-07  +  12.7  (min(Pr2n  -  n  T 

kP  ~  CF  C/78)  RepPr 


for  Gnielinski  correlation: 


1  +  12.7  f/78)1/W/3-l)- 
ki  CF  C/78)  (ReD  -  1000)  Pr 


(6.26) 


(6.27) 


(6.28) 


For  each  of  the  three  calibrated  Nusselt  correlations,  the  thermal  conductivity  equation  is  of  a  different 
form.  The  final  equation  for  each  correlation  is  presented  below.  It  is  important  to  note  that  the  Prandtl  number 
used  in  each  theoretical  correlation  still  contains  the  thermal  conductivity,  k.  Each  of  the  three  correlations 
handles  this  problem  differently. 


The  kd  equation  from  the  Dittus-Boelter  correlation  is: 


,  r  h‘D  -4 
kd  =  *■  <?Rejf  (\iCp)bi  *  b 


The  Pr  equation  from  the  Petukhov  and  Popov  correlation  is: 


»«•»).  1.07 

~  12.7  CheD  127  078) 


ITT  —1=0 


(6.29) 


(6.30) 


However,  since  k  is  implicit  in  Pr,  k  must  be  solved  for  by  solving  for  Pr  first  from  the  above  equation.  In 
other  words,  k  is  divided  into  \iCp  to  obtain  Pr  which  is  then  treated  as  the  unknown  in  the  above  equation.  This 
equation  is  not  an  explicit  linear  type,  but  it  can  be  solved  by  a  numerical  iteration  such  as  the  Newton-Raphson 
method.  Once  Pr  is  determined,  the  thermal  conductivity,  *,  can  be  obtained  from  \iCp  /  Pr. 


*  The  Pr  equation  from  the  Gnielinski  correlation  is: 

As  before,  the  k  must  be  solved  implicitly.  The  resulting  Pr  equation  is: 
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100(0 pO ^  , - 1__1  =  0 

Pr^  "  12.7  CA*Z>  12.7  078)1/2 


(6.31) 


Again,  Pr  in  this  equation  must  be  determined  by  a  numerical  iteration  method. 


APPROACH  2:  THE  PRANDTL  NUMBER  METHOD 

In  the  previous  section,  the  Nusselt  number  method  was  discussed  for  calculating  thermal  conductivity  and 
the  concept  of  the  relative  accuracy  of  the  Nusselt  number  method  by  a  calibration  fiinctioa  This  section  will  follow 
simitar  concepts  but  bypass  the  Nusselt  number  to  discuss  the  Prandtl  method  for  thermal  conductivity  calculations. 
Unlike  the  hot-wire  transient  method  for  thermal  conductivity  measurements,  this  method  has  only  recently  been 
applied  to  thermal  conductivity  measurements.  This  method  measures  the  Prandtl  number,  and  the  thermal  conductivity 
is  calculated  by  knowing  the  Prandtl  number,  viscosity,  and  specific  heat.  The  mathematical  approach  follows. 


Prandtl  Number  and  Other  Related  Parameters 


Based  on  past  studies,  the  Nusselt  number,  Nu^  in  a  turbulent  channel  flow  was  known  as  function  of  the 
Reynolds  number,  ReD,  and  the  Prandtl  number,  Pr.  From  the  definition  of  Nusselt  number,  NueD  =  hD  Ik  and  the 
form  of  the  Dittus-Boelter  correlation,  the  following  equation  is  rewritten  to  correlate  heat  transfer  for  a  turbulent 


flow  heated  by  the  tube  wall: 


W-  =  C-Refi-Pr* 


(6.32) 


where  h  is  the  average  heat  transfer  coefficient  which  equals  q„et  lAUT^.  By  introducing  m,  Cp,  and  the 
temperature  difference  of  the  fluid  at  the  inlet  and  outlet,  the  above  equation  can  be  rearranged  as: 


mCp<T0-fjD 

kA(Tw-Tf) 


C  •  Re%  •  Pr * 


(6.33) 


After  substituting  the  definitions  of  the  Prandtl  and  Reynolds  numbers,  this  equation  can  also  be  rewritten  as: 


T*'T Y 


(6.34) 


In  other  words,  the  Prandtl  number  is  expressed  as  a  function  of  ReD  and  a  dimensionless  temperature  group, 
(To  -  TJ  /  (Tw  -  Tj).  If  this  dimensionless  temperature  group  is  defined  as  A 7*,  the  relationship  is  expressed  as: 


Pr  =  F  {Re D ,  AT*) 


(6.35) 
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This  is  the  relation  for  the  measured  Pr.  It  should  be  noted  that  all  the  quantities  on  the  right  hand  side  of 
Equation  6.34  can  be  obtained  from  experimental  measurements.  These  involve  m,  p,  7^  7),  and  Ta.  With  the 
measurement  of  these  quantities,  the  Prandtl  number  can  be  obtained.  However,  for  fluids  with  known  properties, 
the  Pr  can  also  be  calculated  from  \iCp  /  k.  Therefore,  the  thermal  conductivity  can  then  be  calculated  from  Pr. 
This  method  is  simpler  than  the  Nusselt  number  calibration  method. 

Prandtl  Number  Expression 


The  Pr  relationship  has  been  shown  in  Equation  6.35  as  a  function  of  ReD  and  A  7* .  An  expression  for  Pr 
in  terms  of  Re D  and  A7*  is  required  for  easy  use,  and  it  requires  calibration.  This  calibration  for  Pr  is  very 
important  to  the  accuracy  of  the  thermal  conductivity  measurement.  The  expression  of  Pr  with  Re^  and  AT*  in  a 
closed  form  is  very  difficult  to  obtain  because  it  requires  a  number  of  test  data.  However,  one  of  the  represented 
approaches  for  this  function  is  to  directly  express  Pr  as  a  function  of  ReD  and  AT*  in  a  polynomial  form.  An 
example  of  a  three-degree  polynomial  function  is  given  as: 

Pr  =  a0  +  often  +  a2AT*  +  afte*  +  a4(ReD  •  AT*)  (6.37) 

+  a5AT*2  +  aftel  +  a^Re*  •  AT*) 

+  at(ReD  •  AT*2)  +  a9(A7*3) 

where  the  coefficients,  at,  can  be  fitted  from  measured  data  for  fluids  of  known  properties. 

UNCERTAINTY  ANALYSIS 

This  section  describes  the  uncertainties  in  the  thermal  conductivity  measurements  by  the  methods  of 
Approach  1  and  Approach  2.  A  propagation-of-error  method  [79]  is  used  to  estimate  the  experimental 
uncertainties  associated  with  each  experimental  measurement.  The  sensor  and  equipment  uncertainties  are  listed 
in  Table  6.1.  The  derivation  of  all  related  uncertainty  equations  is  also  presented  in  Appendix  B,  which  includes 
the  thermal  conductivity  uncertainty  by  Approach  1  (three  correlations)  and  Approach  2,  heat  exchange  rate 
uncertainty,  ReD  uncertainty,  friction  factor  (/)  uncertainty,  heat  transfer  coefficient  ( h )  uncertainty,  AT* 
uncertainty,  and  Pr  uncertainty. 


Sensor  and  Geometry  Uncertainty 


Based  upon  the  uncertainties  in  the  measured  data  (both  sensor-  and  geometry-related),  which  is  used  to 
determine  thermal  conductivity,  the  uncertainty  in  k  calculated  from  the  three  different  correlations  can  be 
determined.  Therefore,  the  uncertainty  analysis  for  various  ranges  of  Pr  and  ReD  was  performed. 

Sensor  and  geometry  uncertainties  are  listed  in  Table  6.1. 


TABLE  6.1:  SENSOR  AND  GEOMETRY  UNCERTAINTY 


'  sources  .  . 

length,  L 

0.1  mm 

diameter,  D,  D0 

0.1  mm 

mass  flow  rate,  m 

0.15%  in  kg/s 

temperature,  T,,  T0,  Tw 

0.05°C 

viscosity,  p 

2.0%  of  measurement 

specific  heat  at  constant  pressure,  Cp 

2.0%  of  measurement 

Uncertainty  Estimation  and  Comparison  Between  Approaches  1  and  2 

Based  upon  the  existing  sensor  and  equipment  uncertainties,  the  thermal  conductivity  uncertainties  were 
estimated  and  compared.  In  Approach  1,  the  thermal  conductivity  calculated  from  the  three  different  correlations 
was  estimated  for  a  typical  operating  condition.  The  results  are  tabulated  in  Table  6.2.  Table  6.2  shows  the 
uncertainty,  presented  as  the  percentage  of  measured  thermal  conductivity,  by  Approach  1  for  the  three  Nusselt 
number  correlations  (see  Appendix  B).  The  uncertainty  analysis  was  based  upon  the  sensor  and  dimensional 
uncertainties  listed  in  Table  6.1  for  a  particular  operating  condition  listed  at  the  bottom  of  Table  6.2.  The 
uncertainty  varies  not  only  with  the  correlation  used  but  also  with  various  ReD  and  Pr. 


TABLE  6  2:  THERMAL  CONDUCTIVITY  UNCERTAINTY  PERCENTAGES 
FOR  THREE  NUSSELT  NUMBER  CORRELATIONS 
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From  the  uncertainty  results  shown  in  Table  6.2,  it  is  learned  that: 

l  At  a  low  Pr,  when  ReD  is  higher,  the  uncertainty  is  larger.  However,  the  uncertainty  of  k  through  the 
Dittus-Boelter  correlation  seems  less  affected  by  ReD  than  the  other  two  correlations. 

2.  At  a  moderate  ReD,  Pr  significantly  affects  the  uncertainty  of  k  for  the  Petukhov  and  Popov 
correlation  and  the  Gnielinski  correlation. 

3  At  a  low  ReD,  the  uncertainty  in  k  from  these  three  correlations  is  not  significant  for  a  change  in  Pr. 

4.  From  an  overview  of  the  uncertainty  in  k  from  these  three  correlations,  the  Dittus-Boelter  correlation 
seems  to  be  more  stable  than  the  other  two  correlations,  even  if  these  two  correlations  have  smaller 
uncertainty  values  at  some  particular  operating  conditions. 

Because  CF  is  obtained  from  two  Nusselt  number  ratios  as  mentioned  earlier,  the  only  uncertainty  source  of  CF  is 
from  the  Nusselt  number  curve  fit,  which  will  be  dependent  on  the  curve  fit  method  used  and  how  the  data  are  fit. 
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It  is  interesting  to  compare  the  uncertainty  differences  between  Approaches  1  and  2.  From  an  overview, 
the  variables  involved  in  Approach  1,  regardless  of  the  correlation,  are  more  complex  than  those  in  Approach  2. 
Therefore,  it  implies  that  the  uncertainty  of  Approach  2  should  be  less  than  that  in  Approach  1.  This  point  can  be 
verified  later.  The  thermal  conductivity  uncertainty  percentages  were  calculated  at  some  typical  conditions  of 
operation  for  both  approaches  (see  Appendix  B).  The  following  figures  plot  the  uncertainty  percentage  versus 
7 'W-Tj-.  Figure  6. 1  compares  the  results  at  Tt  -  5°C  and  A T  =  10°C  (A T  =  T0-  Tt),  and  Figure  6.2  compares  the 
results  at  Ti  -  5°C  and  A T= 20°C.  From  these  two  figures,  the  uncertainty  of  Approach  1  seems  twice  as  large  as 
the  uncertainty  in  Approach  2.  They  show  the  same  trends,  and  Tt  and  AT  do  not  seem  to  be  affected  significantly. 
For  other  operating  conditions.  Figure  6.3  compares  the  results  at  Tt  -  10°C  and  AT-  5°C,  and  Figure  6.4 
compares  the  results  at  Tt  =  10°C  and  AT  -  20°C.  Again,  from  these  two  figures,  the.trends  are  exactly  the  same 
as  with  the  first  two  conditions.  Therefore,  these  conclusions  can  be  drawn: 

*  A  7^  significantly  affects  the  uncertainty  of  the  measured  thermal  conductivity  in  both  approaches. 

*  Ti  and  AT  have  less  effect  on  the  uncertainty  of  the  measured  thermal  conductivity  in  both  approaches. 

*  Approach  2  has  less  uncertainty  than  Approach  1. 

SUMMARY  OF  APPROACHES 

Two  methods  were  discussed  for  the  calculation  of  thermal  conductivity:  the  Nusselt  number  calibration 
method  and  the  Prandtl  number  method.  In  the  Nusselt  number  calibration  method,  the  uncertainty  is  caused  by 
measured  h,  m,  p,  Cp,  and  CF  curve  fittings.  However,  in  the  Prandtl  number  method,  the  uncertainty  sources  are 
only  the  ReD  (which  depends  on  m  and  p),  the  AT*,  and  the  Pr  curve  fits.  The  accuracy  of  the  thermal  conductivity 
measurement  can  be  calculated  from  the  uncertainty  of  each  source.  However,  there  are  fewer  variables  used  for 
the  curve  fits  in  the  Prandtl  number  method  than  in  the  Nusselt  number  method.  Therefore,  the  Prandtl  number 
method  is  thought  to  be  simpler  and  more  convenient  for  obtaining  thermal  conductivity  values. 

Based  upon  the  uncertainty  analysis,  the  uncertainty  of  Approach  2  is  only  half  of  that  for  Approach  1. 
Approach  2  is  therefore  believed  to  be  more  accurate  than  Approach  1 . 
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Figure  6.2:  The  uncertainty  estimation  at  T,=  5'C  and  AT-  20*C. 


Temperature  Difference  Between  Wall  and  Fluid  (Tw— Tf),  (°C) 


CHAPTER  7 


EXPERIMENTAL  DATA  CALCULATIONS 

In  this  chapter,  the  experimental  analysis  of  the  raw  data  is  presented  and  the  heat  transfer  coefficient 
calculation  is  described.  Additionally,  the  thermal  conductivity  calculations  from  experimental  data  are  also 
described.  In  the  data  analysis,  each  measured  quantity  was  obtained  by  the  average  of  20  scans,  and  a  FORTRAN 
program  was  then  used  to  reduce  the  data.  Appendix  F  contains  a  copy  of  the  data  analysis  program.  For 
regression  of  multiple  variables,  a  SAS  [80]  program  was  also  used  for  aiding  in  the  correlation. 


DATA  REDUCTION 

The  data  analysis  procedure  for  the  calibration  refrigerants  (HCFC-22,  CFC-12,  CFC-113,  and/or  CFC-114) 
and  the  measured  refrigerant  (HFC-236ea)  is  different  for  some  treatments.  However,  for  the  raw  data  reduction,  it 
follows  the  same  procedure. 

Inner  Wall  Temperature  Corrections 


In  the  test-tube  used,  the  surface  temperature  is  not  measured  on  the  actual  inner  tube  wall  surface  because 
of  the  difficulty  of  installing  thermocouple  beads  there.  Instead  of  installment  on  the  inner  surface,  all 
thermocouple  beads  for  measuring  wall  surface  temperature  were  buried  midway  into  the  wall.  In  order  to  obtain  a 
more  realistic  surface  temperature,  a  correction  is  required  to  obtain  the  real  inner  wall  surface  temperature.  A 
correction  equation  was  developed  by  a  heat  conduction  model  and  is  described  as  follows: 


Qnel 


2nLks  (Tr2  -  Tr  l) 
In  (r2/r,) 


(7.1) 


where  qnet  is  net  heat  input  rate,  L  is  tube  length  (2  m  in  current  test-tube),  r2  is  the  tube  radius  at  the  point  where 
thermocouple  beads  are  buried,  is  the  inside  radius  of  the  tube,  ks  is  the  thermal  conductivity  of  solid  tube  (ks  = 
401  H7(/«  •  K )  for  pure  copper),  and  Tr  lt  and  Tr^  denote  the  temperature  at  r,  and  r2,  respectively.  A  correction 
temperature  is  then  expressed  as: 
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T  inethW  1) 

Arc  "  2 HkT 

Therefore,  the  real  inner  wall  temperature,  Tr  x  is  corrected  by: 

TrA  =  Trfi~*Tc 


Tra  in  this  case  is  the  measured  temperature  by  a  thermocouple. 


Average  Quantity  Calculations 

Average  quantities  include  a  timed  average  which  indicates  the  average  of  measured  quantities  in  a 
certain  scanning  time  interval,  and  a  location  average  which  indicates  the  average  of  measured  quantities  in 
different  positions.  These  average  measured  quantities  are  calculated  as  follows. 

•  Raw  Data  Average 

In  order  to  obtain  a  more  average  representation  of  the  experimental  data,  20  scans  were  taken  for  a 
steady-state  situation,  and  these  20  values  were  then  averaged  for  each  representative  quantity.  These  quantities 
include  the  following  variables:  Tt,T0,Twi,Ts,Ta,m  ,Q ,  and  p.  Therefore: 

average  quantity  =  ”[^1^  ( quantity ),  ]  (7.4) 

where  n  stands  for  the  total  number  of  scanning  times  (n  =  20),  and  /  represents  the  i-th  scan. 

•  Average  A  7^  Calculations 

The  average  fluid  temperature  in  the  test-tube,  7^  was  therefore  taken  as  the  average  temperature  of  the 
inlet  and  outlet  temperature,  ?,  and  Ta,  of  the  test-tube.  However,  in  the  average  A 7^  calculation,  which  is  Tw  -  Tf, 
the  finite  control  volume  method  was  used  which  considers  the  local  heat  transfer  characteristics  (thermal  entrance 
effect)  rather  than  the  all-points  average  method.  Equation  7.5.  The  mathematical  derivation  is  presented  as 
follows.  For  the  tube  wall  and  fluid  temperatures,  the  averages  are  calculated  as  the  following: 


<£><* 


w1  +  *’w2+,"  +  7’w13) 


Tf- 


Equations  7.5  and  7.6  are  used  for  obtaining  an  average  wall  temperature  which  assumes  a  linear  wall 
temperature  distribution.  It  is  applicable  for  fully  developed  regions  subjected  to  a  constant  heat  flux.  However, 


72 


the  entrance  effect,  in  some  situations,  is  significant  [29,76],  In  those  cases,  the  wall  temperature  distribution  in 
that  region  is  no  longer  linear.  Therefore,  if  a  linear  average  of  the  wall  temperature  is  taken  for  the  average  wall 
temperature,  that  will  affect  the  accuracy  of  the  heat  transfer  coefficient  calculation.  Another  approach,  the  finite 
control  volume  method  mentioned  above,  is  based  on  an  energy  balance  which  provides  a  better  representation 
than  the  linear  average  method.  This  method  is  applied  under  the  assumption  of  a  constant  heat  flux  situation.  In 
this  situation,  the  amount  of  heat  added  to  a  local  section  is  everywhere  equal  as  long  as  the  surface  area  is  equal. 
That that  an  equal  amount  of  heat  is  obtained  for  an  equal  length  of  heated  section  in  a  constant  cross 
station  area  tube.  An  energy  balance  equation  for  a  heated  section  can  be  expressed  as: 

Q  =  9j +  ?2  +  “* +  9, + ~  +  9«  (7’7) 

where  the  subscript  n  represents  the  total  number  of  subsections  and  i  represents  the  i-th  subsection. 

Fart,  q.  would  be  the  same  due  to  the  constant  heat  flux,  if  each  subsection  has  an  equal  surface  area  and 
the  same  property  values.  The  above  equation  can  also  be  written  as: 

Q  =  q"  (j4j  +A2  +  “*  +  •••  +An  )  (7.8) 

However,  since  Q  is  equal  to  hA&T^  the  above  equation  can  be  rewritten  as: 

t 

A T^f  hjdA  =  q"  (Al+A2  +  —  +Aj  +  —  +A„)  (7.9) 

or  rewritten  into  the  finite  control  volume  form, 

A7^  X  hlAi  =  q  (A}  +A2  +  ~  +At  +  —  +A„  )  (7.10) 

i-1 

Based  upon  the  assumption  of  a  constant  heat  flux  in  each  subsection.  A,  is  equal  to  qn  /  A 7^.  The  A T^. 
is  a  locally  averaged  quantity,  which  is  Tw.  -  7 Jr.  Substituting  A,  into  above  equation, 

A^r  ^  AT^l  =  q"  (■A1+A2  +  ~+Ai  +  ~+A„)  (7.11) 

L  <-l 
or, 

ATtfl  X  *  A1+A2  +  -+Ai  +  ~+A„)  *  A  (7.12) 

L 

where  A  is  the  total  heat  transfer  area. 

Finally  the  average  temperature  difference  between  the  average  wall  temperature  and  fluid  temperature 
can  be  expressed  as: 
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(7.13) 


Or,  in  another  form. 


(7.14) 


where  A/,  is  the  length  of  i-th  section. 

For  the  current  test  section,  the  wall  temperature  is  measured  at  13  locations.  Finite  control  volumes, 
therefore,  can  be  divided  into  13  subsections.  A  typical  wall  temperature  and  average  fluid  temperature 
distribution  is  shown  in  Figure  7.1.  Because  of  the  entrance  effect,  the  wall  temperature  at  the  inlet  portion  should 
reflect  a  non-linear  distribution.  For  a  better  representation,  smaller  subsections  were  taken.  Therefore,  the 
determination  of  subsections  at  the  inlet  portion  was  based  on  the  wall  temperature  at  the  measured  points.  Figure 
7.2  shows  the  finite  control  volume  configuration  of  the  current  test  section.  In  this  figure,  the  control  volumes  in 
the  inlet  and  outlet  portion  were  smaller  than  the  middle  portion  because  more  temperatures  were  measured  at  the 
ends  of  the  test  section.  This  method  for  obtaining  the  average  temperature  difference  was  thought  to  be  a  better 
representation  than  average  of  the  all-points  method  because  it  accounts  for  the  thermal  behavior  of  the  entrance 
length.  The  subsection  length  distribution  is  listed  in  Table  7.1. 


TABLE  1. 1 :  SUBSECTION  LENGTH  DISTRIBUTION 


S.v  subsection*#? 

m 

•  point  111 

Ip^ceftoro-; 
:;P::itilet.  Cnii||| 

:fj|  Aff,  cm 

i 

1 

o  ; 

5 

2 

2 

10 

10 

3 

3 

20 

15 

4 

4 

40 

20 

5 

5 
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140 
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10 

10 

160 

20 

11 

11 

180 

15 

12 

12 

190 

10 

13 

13 

200 

5 
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uniform  heat  flux  q" 


Figure  7.1:  Typical  temperature  distribution  along  a  heated  test  section. 
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Thermocouple  Position  Along  Tube 


Figure  7.2:  Finite  control  volume  configuration  of  the  test  section. 
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Heat  Transfer  Calculations 


One  of  the  thermodynamic  properties,  specific  heat  (Cp),  is  used  in  calculating  a  single-phase  heat 
transfer  quantity,  based  upon  the  energy  balance  principle.  This  property,  in  general,  is  function  of  temperature 
and  pressure.  However,  for  a  liquid  at  constant  pressure,  Cp  is  a  function  of  temperature  only.  But,  liquid  Cp  is 
not  a  strong  function  of  temperature  except  above  the  reduced  temperature,  Tr ,  of  0.7  to  0.8.  Detailed  information 
on  the  description  of  Cp  is  given  in  thermodynamic  and  property  textbooks  [28,81]. 

The  net  heat  input  rate  for  refrigerants  of  known  Cp  was  calculated  directly  from  the  energy  equation  for 
flow  in  the  test-tube.  That  is,  for  refrigerants  with  known  Cp: 

(7.15) 

For  refrigerants  with  unknown  Cp,  Qnet  is  calculated  from  the  energy  balance  of  the  test  section.  That  is,  for 
refrigerants  with  unknown  Cp,  Equation  6.5  is  applied.  Then  the  q„et  can  be  calculated  from  Equation  6.3. 

<716> 

In  Equation  6.3,  qhss  is  the  heat  loss  of  the  test  section.  A  correlation  for  calculating  qloss  was  described  in  Equation 
5.1  in  Chapter  5.  The  heat  transfer  coefficients  can  be  obtained  by  Equation  6.2.  The  parameters  in  this  study, 
including  Nup,  Pr,  Rep,  and  A7*.  were  then  calculated  by  their  definitions  presented  earlier  or  defined  in  Chapter  6. 

Viscosity  and  Thermal  Conductivity  Calculations 

The  product  of  viscosity  and  density  is  directly  measured  by  the  viscometer.  The  viscosity  is  then  obtained 
by  dividing  it  by  density  which  is  measured  by  a  densimeter  or  obtained  from  reliable  data  sources,  such  as  the 
ASHRAE  Handbook  [82]  or  REFPROP  [83],  which  are  shown  in  Appendix  A. 

There  are  two  approaches  proposed  in  this  study.  The  thermal  conductivity  can  be  obtained  using  either 
one.  However,  except  for  thermal  conductivity  being  directly  calculated  from  the  Dittus-Boelter  correlation  in 
Approach  1  (Equation  6.29)  the  other  thermal  conductivity  correlations  (Equations  6.30, 6.31,  or  6.34)  are 
developed  from  Pr.  Therefore,  Pr  must  be  obtained  before  the  thermal  conductivity  can  be  calculated. 

CORRECTION  OF  THE  PRESSURE  EFFECT  ON  THERMAL  CONDUCTIVITY  AND  VISCOSITY 

Pressure  has  some  effect  on  the  thermal  conductivity  and  viscosity  of  a  liquid.  However,  these  effects  are 
usually  significant  only  at  high  pressures.  Thermal  conductivity  and  viscosity  are  not  only  strong  functions  of 
temperature,  but  they  are  also  strong  functions  of  pressure  when  a  liquid  is  subjected  to  a  high  pressure.  The 
following  discussion  provides  corrections  for  the  pressure  effects  on  thermal  conductivity  and  viscosity. 
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Pressure  Effect  on  Liquid  Thermal  Conductivity 


In  general,  at  moderate  pressures,  up  to  5  or  6  MPa,  the  effect  of  pressure  on  the  thermal  conductivity  of 
liquids  is  usually  neglected,  except  near  the  critical  point.  Missenard  has  published  a  simple  correlation  for  k 
which  extends  to  high  pressures  [28],  The  thermal  conductivity  ratio  is  subjected  to  two  different  pressures  but  at 
the  gam*  temperature.  This  ratio  was  found  to  correlate  with  reduced  temperature  and  pressure  as  shown  in  the 
following  equation: 


k(Pr) 

k  ( low  pressure) 


1  +  Q-  P?1 


(7-17) 


where  Pr  is  reduced  pressure  and  Q  is  a  function  of  Prand  Tr  Q  is  given  in  Table  10-8  in  R.C.  Reid  et  al.  [28].  A 
plot  of  thermal  conductivity  ratio  against  Pr  is  given  by  Figure  7.3. 


Pressure  Effect  on  Liquid  Viscosity 


Liquid  viscosity  will  be  affected  by  pressure,  especially  at  high  pressures.  Increasing  the  pressure  over  a 
liquid  results  in  an  increased  viscosity.  However,  at  pressures  less  than  the  critical  pressure,  the  pressure  effect  on 
viscosity  is  very  weak.  Lucas  suggested  that  the  change  might  be  estimated  from  Equation  7.18  [28]: 


n  _  1  +  D  ( APr/2.118  y4 

”  1  +  Ca£J>r 


(7.18) 


where 

|x  «  viscosity  of  liquid  at  pressure  of  interest 
pMr  =  viscosity  of  liquid  at  saturated  state 
a  =  acentric  factor 
APr  =  (P-P^IP, 

P  -  liquid  pressure 
P^  ■  vapor  pressure 
Pe  *  critical  pressure 
Tr  =  reduced  temperature 

A  «  0.9991  -  [  4.674  •  10"4  /  (1.0523  T^  03877  -  1.0513  )  ] 

D  =  [  0.3257  /  (  1.0039  -  Tr2  m  )0  2906  -  0.2086 

C  =  -0.07921  +  2.1616  Tr-  13.4040  T2  +  44.1706  Tr3  -  84.8291  Tr4  +  96.1209  Tr5 
-59.812  Tr6  +  15.6719  T? 

A  plot  of  viscosity  ratio  versus  A Pr  is  shown  in  Figures  7.4  and  7.5. 
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Reduced  Differential  Pressure,  APr=(P-Pvp)/P( 


Thermal  Conductivity  and  Viscosity  Corrections  for  Pressure  Effects 


High  pressure  (higher  than  critical  pressure)  has  a  strong  effect  on  viscosity  and  thermal  conductivity,  while 
low  pressure  has  almost  no  effect  or  only  a  minor  effect  on  them,  as  discussed  in  the  above  sections.  In  the  current 
study,  refrigerants  were  circulated  in  the  low  pressure  range.  Therefore,  pressure  had  a  very  small  effect  on  viscosity 
and  thermal  conductivity.  However,  to  obtain  accurate  properties,  pressure  effects  on  both  viscosity  and  thermal 
conductivity  were  still  taken  into  account.  The  manner  in  which  pressure  effects  were  considered  is  outlined  in 
Equations  7.17  and  7.18  for  thermal  conductivity  and  viscosity,  respectively.  Because  the  operating  state  must  be 
controlled  with  limited  subcooling  at  the  outlet  of  the  test  section,  the  system  pressure  should  also  be  greater  than 
the  saturation  pressure  to  maintain  a  liquid  phase  throughout  the  test  section.  The  viscosity  and  thermal 
conductivity  at  a  selected  pressure  were  obtained  from  those  at  the  saturation  pressure  by  using  Equations  7.17  and 
7. 18.  By  reversing  the  usage  of  these  equations  and  by  knowing  viscosity  and  thermal  conductivity  at  a  specific 
pressure,  then  the  viscosity  and  thermal  conductivity  at  the  saturation  pressure  can  also  be  obtained. 

DATA  PRESENTATION 

This  section  describes  the  format  used  to  present  the  data  in  order  to  identify  trends  for  the  purpose  of  data 
comparison.  In  general,  the  thermal  conductivity  and  viscosity  are  strong  functions  of  temperature  and  weak 
functions  of  pressure  in  a  low  pressure  region.  Therefore,  the  thermal  conductivity  and  viscosity  were  presented 
only  as  a  function  of  temperature.  All  the  plots  shown  in  this  study  were  plotted  as  saturated  liquid  thermal 
conductivity  and  viscosity  versus  temperature.  Measured  data  compared  with  the  data  in  the  ASHRAE  Handbook 
and  from  REFPROP  were  also  shown  in  the  plots.  The  measured  specific  heat,  Cp,  and  density,  p,  of  interest 
versus  temperature,  T,  were  also  plotted  and  compared  with  ASHRAE  Handbook  [82]  or  REFPROP  [83]  values. 
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CHAPTER  8 


CALIBRATIONS  OF  THERMAL  CONDUCTIVITY  MEASUREMENT 

In  this  chapter,  calibration  of  the  thermal  conductivity  measuring  apparatus  by  using  refrigerants  of 
known  properties  will  be  shown  and  discussed  for  both  approaches.  The  calibration  function  (CF)  in  Approach  1 
and  the  correlation  obtained  by  Approach  2  will  be  presented.  In  addition,  the  accuracy  of  the  measured  data  by 
these  approaches  will  also  be  presented. 


CALIBRATION  OF  REFRIGERANTS  AND  OPERATING  RANGES 

Before  T"<re»»'™g  the  thermal  conductivity  of  HFC-236ea,  the  calibration  functions  should  be  obtained. 
The  calibration  functions  have  been  obtained  from  several  fluids  with  known  properties  by  the  method  proposed  in 
this  study.  Several  refrigerants  of  known  properties  were  selected  in  order  to  obtain  the  calibration  functions. 

They  were  HCFC-22,  CFC-12,  CFC-113,  and  CFC-1 14.  The  properties  used  for  these  refrigerants  are  from 
ASHRAE  data,  and  their  curves  are  presented  in  Appendix  A.  Table  8. 1  shows  the  operating  ranges  used  to  obtain 
the  calibration  functions  of  these  selected  refrigerants.  The  fluids  (HFC-236ea,  blends  A  and  B)  of  unknown 
thermal  conductivity  must  be  tested  within  the  Pr  number  ranges  of  the  basic  refrigerant  data  shown  below.  For 
HFC-236ea,  as  an  example,  with  operating  temperatures  between  0°C  and  50°C,  the  Pr  number  varies  from  5.60  to 
8.12,  while  ReD  can  be  controlled  by  regulating  the  mass  flow  rate  within  the  ranges  for  the  basic  refrigerants 
listed  above.  Under  most  conditions,  the  mass  flow  rate  varied  from  3  kg/min.  to  15  kg/min. 

TABLE  8.1:  PRANDTL  NUMBER,  REYNOLDS  NUMBER,  AND  A7*  RANGES  FOR  SELECTED  FLUIDS 


WmMJEMmM 

AT* 

HCFC-22 

2.28  -  2.56 

38360-184852 

0.81  -  1.42 

CFC-12 

2.82-3.08 

31031  -  168007 

0.97  -  1.30 

CFC-113 

7.00  - 10.58 

7970-  64354 

0.64  -  0.80 

CFC-1 14 

4.46  -  6.34 

16296  -  99517 

0.79  -  0.99 

overall 

2.28  - 10.58 

7970  -  184852 

0.64-1.42 
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CALIBRATION  FUNCTIONS 


The  calibration  functions  were  found  for  the  three  Nusselt  correlations  mentioned  earlier.  A  linear  regression 
between  ln(Nu  /  Pr")  and  In(ReD)  was  performed  for  the  experimental  curve  fits.  In  order  to  verify  the  current  method 
for  thermal  conductivity  measurement,  the  calibration  refrigerants  can  be  selected  such  that  the  calibration  functions 
were  obtained  only  for  HCFC-22,  CFC-12,  and  CFC-1 13,  which  covered  the  Pr  number  range  from  2.28  to  10.58. 
These  fiinrtinnc  could  then  be  applied  to  the  calculation  of  unknown  fluids  such  as  CFC-1 14.  Figures  8.1, 8.2,  and 
8.3  show  the  results  for  the  three  different  Nusselt  number  correlations  used  for  these  calibration  refrigerants. 

The  calibrated  Nusselt  number  can  be  obtained  from  the  curve-fitted  equations  such  that: 

-  CF  <*•» 

where  the  calibration  function,  CF  can  be  obtained  from  Equation  6.24. 

As  mentioned  in  Chapter  6,  the  Nusselt  number  can  be  obtained  from  the  Dittus-Boelter  correlation 
(Equation  6.6),  the  Petukhov  and  Popov  correlation  (Equation  6.7),  and  the  Gnielinski  correlation  (Equation  6.8). 
Then,  the  Nusselt  number  can  be  expressed  as  the  form  of  ReD  and  Pr  (Equations  6.22  and  6.23). 

The  coefficients  from  curves  of  the  calibration  functions  are  shown  in  Table  8.2  for  hD/k,  in  which  h  is 
experimentally  measured,  and  for  the  three  different  Nusselt  number  correlations. 


TABLE  8.2:  CURVE-FIT  COEFFICIENTS  FOR  NUSSELT  NUMBER 
-  THREE  CALIBRATION  REFRIGERANTS 


^iciefficlenf® 

%m 

bitfuihBoelter 

Petukhov  1 1 

*ftGn Ietiii&ff! 

c.c* 

0.006758 

0.023 

0.017266 

0.012093 

a,  <f 

0.903770 

0.8 

0.827298 

0.862783 

b,b* 

0.502199 

0.4 

0.455283 

0.441244 

It  should  be  noted  that  the  calibration  function,  CF,  is  still  expressed  as  function  ofReD  and  Pr. 

Therefore,  like  the  NuD  expression  in  Equation  6.20,  a  linear  relationship  between  IniCF/Pr**  ~b)  and  ln(ReD ) 
holds  true  and  can  be  described  as  follows: 

ln(CF/Pr **  -  b)  =  /n(C*/Q  +  (a9  -  a)  ln(ReD)  (8.2) 

The  plots  of  ln(CF)  versus  ln(ReD)  are  shown  in  Figures  8.4, 8.5,  and  8.6  for  the  CF  obtained  from  the 
Dittus-Boelter,  Petukhov  and  Popov,  and  Gnielinski  equations,  respectively.  As  indicated  in  these  figures,  the 
relationships  are  all  linear  with  the  various  Pr. 
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Figure  8.1:  Nusselt  number  calculated  from  hDIk  and  Dittus-Boelter 
correlation  for  the  calibration  refrigerants. 
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Figure  8.2:  Nusselt  number  calculated  from  hD/k  and  Petukhov  and  Popov 
correlation  for  the  calibration  refrigerants. 
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Figure  8.3:  Nusselt  number  calculated  from  hDIk  and  Gnielinski 
correlation  for  the  calibration  refrigerants. 
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Figure  8.5:  Calibration 


Figure  8.6:  Calibration  functions  for  the  Gnielinski  correlation. 


If  the  four  refrigerants,  HCFC-22,  CFC-12,  CFC-113,  and  CFC-114,  were  used  as  the  calibration 
refrigerants,  the  calibration  function,  CF,  would  be  slightly  different  but  very  close  to  that  obtained  for  the  three 
calibration  refrigerants.  The  CF  coefficients  are  listed  in  Table  8.3. 


TABLE  8.3:  CURVE-FIT  COEFFICIENTS  FOR  NUSSELT  NUMBER 
-  FOUR  CALIBRATION  REFRIGERANTS 


liicoefncient 

HD(k  1 

iillipiiiiel  l 

'Petukhov  1! 

GniefiMki 

c.cr 

0.006633 

0.023 

0.017372 

0.012381 

a, 

0.905027 

0.8 

0.826843 

0.860922 

b,  b* 

0.507640 

0.4 

0.457499 

0.442776 

THE  PRANDTL  NUMBER  REGRESSION 

As  mentioned  earlier.  Approach  2  for  calculating  the  thermal  conductivity  uses  the  Prandtl  number,  Pr. 
The  Pr  is  determined  as  function  of  ReD  and  AT*  as  pointed  out  previously.  For  the  three  basic  calibration 
refrigerants:  HCFC-22,  CFC-12,  and  CFC-1 13  and  the  four  basic  calibration  refrigerants:  HCFC-22,  CFC-12, 
CFC-113,  and  CFC-114,  the  full  term  polynomial  function  regressions,  based  on  Equation  8.3,  for  ln(Pr)  were 
fitted  as  shown  in  Table  8.4. 

ln(Pr)  =a0  +  axln(ReD)  +  a2ln(AT*)  +  a3ln(Rep  +  aAln(ReD)  •  ln(AT*)  +  a5ln(AT*2) 

+  a6ln(Rep  +  a-jln(Refy  •  ln(AT*)  +  agln(ReD)  •  ln(AT*2)  +  a9In(AT+2)  (8.3) 


TABLE  8.4:  CURVE-FIT  COEFFICIENTS  FOR  ln(Pr) 


•'coefficient 

3  refrigerants 

4  refrigerant* 

°o 

96.908792 

132.495075 

o, 

-26.452258 

-36.446721 

a2 

83.304684 

102.201609 

°3 

2.459228 

3.392931 

*4 

-16.020859 

-19.449363 

a3 

8.088965 

14.165714 

°6 

-0.076756 

-0.105753 

0.750949 

0.905506 

-0.668378 

-1.253316 

a9 

2.213377 

2.823725 
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Non-dimensional  Temperature,  AT* 


Plots  of  Pr  versus  AT*  are  shown  in  Figures  8.7  and  8.8,  for  the  three-refrigerant  and  four-refrigerant  bases.  As 
be  seen  from  these  figures,  the  smaller  the  AT*,  the  larger  the  Pr  for  a  fixed  Re^.  On  the  other  hand,  at  a  fixed 
AT*,  the  smaller  the  ReD ,  the  larger  is  the  Pr.  These  plots  furnish  the  trends  of  the  relationship  between  Pr,  Re^ 
and  AT*  developed  by  Approach  2. 

COMPARISON  OF  MEASURED  THERMAL  CONDUCTIVITY  WITH  REFERENCE  DATA 

It  is  of  interest  to  compare  the  measured  thermal  conductivity  with  the  ASHRAE  handbook  data  [82]  for 
the  calibration  refrigerants  (following  calibration).  To  do  this  comparison,  the  thermal  conductivity  of  the  four 
calibration  refrigerants  were  plotted  using  Approaches  1  and  2.  Figure  8.9  shows  the  measured  thermal 
conductivity  using  Approach  1  and  the  ASHRAE  handbook  thermal  conductivity  plot,  and  Figure  8.10  shows  the 
measured  thermal  conductivity  by  Approach  2  with  the  ASHRAE  handbook  thermal  conductivity.  Lines  of  ±5% 
were  nyd  to  present  an  indication  of  the  deviation  in  the  measured  data.  As  shown  in  these  figures,  most  of  the 
measured  points  were  located  inside  the  ±5%  band  except  for  some  points  scattered  outside  this  band  but  within 
±10%.  This  could  be  due  to  experimental  errors,  which  are  not  included  in  the  experimental  uncertainty,  such  as 
unsteady-state  or  inaccurate  property  data.  By  comparing  the  distribution  of  the  data  points  in  these  two  figures,  it 
is  found  that  there  seems  to  be  more  fluctuation  in  the  results  of  Approach  1  than  in  the  results  of  Approach  2. 
However,  these  results  show  that  the  experimental  uncertainties  were  within  ±5%  as  predicted  by  theoretical 
uncertainties  which  were  pointed  out  earlier.  These  results  provide  more  confidence  in  measuring  other 
refrigerants  by  this  method. 
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ASHRAE  data 


CHAPTER  9 


VERIFICATION  OF  METHODOLOGY 

The  purpose  of  this  study  is  to  determine  the  selected  properties  of  refrigerant  alternatives.  Because  HFC- 
236ea  is  a  proposed  alternative  for  CFC- 1 14,  properties  of  CFC- 1 14  were  measured  by  the  current  approach  and 
compared  with  the  ASHRAE  data  in  order  to  see  whether  this  approach  can  be  applied  to  fluids  of  unknown 
properties.  The  properties  of  CFC-1 14  being  verified  in  this  study  are  thermal  conductivity,  k,  viscosity,  p,  density, 
p,  and  specific  heat,  Cp.  However,  it  is  first  interesting  to  compare  these  properties  between  REFPROP  and  ASHRAE 
data.  Then,  the  measured  results  are  presented  and  compared  with  the  ASHRAE  data.  The  comparisons  are 
shown  in  the  following  sections. 


CFC-1 14  PROPERTIES  FROM  REFPROP  AND  ASHRAE 

REFPROP  is  a  computer  package  for  refrigerant  properties  recently  developed  by  the  National  Institute  of 
Standards  and  Technology  (NIST),  and  REFPROP-4.0  is  the  most  recent  version.  Because  of  the  limited  experimental 
data,  theoretical  predictions  were  thought  to  be  the  only  method  to  obtain  properties  for  this  heat  transfer  study.  In  this 
section,  interest  is  focused  on  a  comparison  of  density,  specific  heat,  viscosity,  and  thermal  conductivity  for  CFC-1 14. 
Figures  9.1, 9.2, 9.3,  and  9.4  show  the  plots  of  the  ASHRAE  data  and  REFPROP-4.0  data  for  these  properties. 

As  shown  in  these  comparison  plots,  some  deviations  exist.  The  deviations  presented  as  percentages  over 
a  temperature  range  of  -50°C  to  120°C  for  each  of  the  four  properties  discussed  here  are  listed  in  Table  9. 1 .  As  shown 
in  these  figures  (Figures  9. 1  through  9.4)  and  Table  9. 1,  these  properties  deviate  somewhat  from  each  other  over  the 
temperature  ranges  under  study.  Therefore,  the  results  obtained  might  be  different  if  other  data  sources  were  used. 

TABLE  9.1:  CFC-1 14  PROPERTY  COMPARISON  BETWEEN  ASHRAE  AND  REFPROP 


§§  property 

iliiiiiiilii 

dev.,%  !§ 
(high  temp.)- 

density,  p,  ( kglm 3) 

+5.2 

-0.02 

specific  heat,  Cp,  tlfikg  *  K)] 

+6.8 

-2.2 

viscosity,  p,  (Pa  •  s) 

+4.2 

-7.1 

thermal  conductivity,  k,  [H7<»» '  Ml 

-2.7 

-11.7 

°dev.,%  =  [(REFPROP  value  -  ASHRAE  value )/ ASHRAE  value  ]  100 
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Comparison  plot  between  ASHRAE  and  REFPROP  data 
for  CFC-114  thermal  conductivity. 


VERIFICATION  OF  CFC-114  PROPERTIES 


In  this  section,  the  measured  properties  of  CFC-1 14  which  include  thermal  conductivity,  viscosity,  specific 
heat,  and  density  will  be  presented  and  compared  with  ASHRAE  data. 

CFC-114  Thermal  Conductivity 

The  CFC-1 14  thermal  conductivities  were  obtained  from  experimental  data  operating  from  near  0°C  to 
50°C.  Figure  9.5  shows  the  thermal  conductivity  of  CFC-114  by  Approach  1  as  a  function  of  temperature  along 
with  the  ASHRAE  data  [82],  The  deviations  of  thermal  conductivity  between  measured  and  ASHRAE  data  are 
shown  to  be  within  ±5%  over  the  measured  temperature  range.  Figure  9.6  shows  the  thermal  conductivity  of 
CFC-1 14  lty  Approach  2  as  a  function  of  temperature  along  with  the  ASHRAE  data.  Again,  the  deviation  is  within 
±5%.  It  is  important  to  mention  that  the  properties  required  in  the  thermal  conductivity  calculations,  such  as 
viscosity  and  specific  heat,  were  all  based  on  ASHRAE  data.  Deviations  of  thermal  conductivity  for  measured  and 
ASHRAE  data  are  shown  in  Figures  9.7  and  9.8  for  Approaches  1  and  2,  respectively. 

CFC-114  Viscosity 

The  CFC-1 14  viscosity  was  obtained  from  the  product  of  viscosity  and  density,  directly  measured  by  the 

inline;  viscometer.  The  accuracy  of  the  viscometer  had  already  been  verified  for  several  fluids.  Figure  5.11  shows 
the  viscosity  as  a  function  of  temperature  for  CFC-1 14.  As  can  be  seen  from  the  figure,  the  measured  viscosity  and 
the  ASHRAE  viscosity  were  in  veiy  close  agreement  with  the  deviation  being  ±2%.  Figure  5.12  shows  the 
deviation  plot  over  the  measured  temperature  range. 

CFC-114  Specific  Heat 

The  measured  specific  heat  was  calculated  from  an  energy  balance  of  the  test  section.  With  the  net  heat 
transfer  calculated,  the  specific  heat  was  calculated  from  Equation  6.5.  Figure  9.9  shows  the  specific  heat  of  CFC- 
1 14  over  a  temperature  range  of  0°C  to  50°C.  As  shown  in  this  figure,  the  measured  Cp  closely  matches  the 
ASHRAE  data.  This  suggests  that  the  temperature  and  the  net  heat  transfer  to  the  refrigerant  values  are  quite  good 
and  provide  much  confidence  in  measuring  unknown  fluids  such  as  HFC-236ea. 

It  should  be  noted  that  the  specific  heat  of  CFC-1 14,  Cp,  was  calculated  by  using  the  test  section  heat  loss 
estimation  described  above.  The  resulting  heat  loss  is  shown  in  Figure  5.7,  where  Equation  5. 1  was  used.  The 
measured  Cp  compared  with  the  ASHRAE  Cp,  shown  in  Figure  9.9,  is  within  ±3%. 
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Comparison  of  measured  thermal  conductivity  by  Approach  1 
and  ASHRAE  data  for  CFC-1 14. 
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Figure  9.6:  Comparison  of  measured  thermal  conductivity  by  Approach  2 
and  ASHRAE  data  for  CFC-1 14. 
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ASHRAE  Thermal  Conductivity,  (  WHm  •  A' )  1 
Fiaure  9  T  Comparison  of  measured  thermal  conductivity  by  Approach  1 
and  ASHRAE  data  for  CFC-1 14. 


CFC-114  Density 


The  density  of  CFC-1 14  was  measured  by  a  densimeter.  The  densimeter  was  verified  using  known 
refrigerant  data,  and  the  results  were  compared  with  the  ASHRAE  data  as  shown  in  Figure  5.8  (Chapter  5).  The 
measured  data  are  shown  in  Figure  9. 10.  The  calculated  deviation  percentage  was  calculated  to  be  within  ±1.9% 
of  the  ASHRAE  data  for  CFC-1 14. 

SUMMARY 

In  this  chapter,  the  measured  properties  of  CFC-1 14  have  been  presented  and  compared  with  ASHRAE 
data.  The  purpose  of  the  CFC-1 14  property  measurements  was  to  verify  the  methodologies  used  here.  The  main 
properties  of  interest  are  thermal  conductivity  and  viscosity.  However,  other  properties,  such  as  specific  heat  and 
density,  were  also  measured  and  compared.  Since  specific  heat  is  a  required  property  for  determining  the  thermal 
conductivity  for  the  methodologies  of  this  study,  the  validity  of  the  Cp  data  impacted  the  accuracy  of  the  calculated  k. 
In  general,  the  measured  properties  closely  matched  the  ASHRAE  data.  This  supports  the  validity  of  the 
methodology  developed  in  this  study.  A  summary  follows  of  this  verification.  Table  9.2  shows  this  summary  of 
verified  results  for  each  property'. 

TABLE  9.2:  SUMMARY  OF  MAXIMUM  DEVIATION  FOR  CFC-1 14  PROPERTIES 


property 

density,  p 

±5.2% 

+1.9% 

m msarmm 

±6.8% 

±3% 

viscosity,  p 

±7.1% 

±2% 

thermal  conductivity*,  k 

±11.7% 

±5% 

*Dev., :  [  (REFPROP  value  -  ASHRAE  value)  /  ASHRAE  value )  •  100 
^Dev.,:  [  (Experiment  value  -  ASHRAE  value)  /  ASHRAE  value  ]  *  100 
‘Tor  both  Approaches  1  and  2 


The  deviations  shown  in  Table  9.2  indicate  that  larger  deviations  exist  between  REFPROP  and  ASHRAE 
data  for  transport  properties  (viscosity  and  thermal  conductivity),  while  smaller  deviations  were  detected  for 
thermodynamic  properties  (density  and  specific  heat).  For  measured  properties,  it  was  shown  that  they  deviate  less 
from  the  ASHRAE  data  than  from  the  REFPROP  predictions.  This  indicates  that  the  methodology  worked 
adequately.  However,  to  improve  accuracy,  more  calibrations  might  be  necessary.  Other  possible  error  sources, 
excluding  experimental  uncertainties  such  as  operation  errors  (e.g.  steady-state  requirement),  must  be  eliminated 
as  much  as  possible  during  the  experimental  operation. 


107 


Legend 
ASHRAE  data 


KX  •**>//*] ‘<b 


108 


Figure  9.9:  Comparison  of  measured  specific  heat  for  CFC-114. 
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Figure  9.10:  Comparison  of  measured  density  for  CFC-114. 


CHAPTER  10 


PROPERTIES  OF  HFC-236ea 

In  the  previous  chapter,  the  measured  properties  of  CFC-1 14  were  compared  with  ASHRAE  data.  The 
small  differences  seen  gave  confidence  for  applying  the  current  approach  to  measure  refrigerant  properties.  Based 
on  the  thermal  conductivity  calibrations  (both  Approaches  1  and  2),  same  test  section  conditions  and  viscometer, 
the  properties  of  HFC-236ea  were  measured.  The  results  were  also  compared  with  REFPROP-4.0  data  [83]  and  are 
presented  below. 

SPECIFIC  HEAT 

Based  on  the  same  test  section  heat  loss  estimation  mentioned  previously,  the  Cp  of  liquid  HFC-236ea  was 
determined.  Figure  10.1  shows  the  Cp  of  HFC-236ea  and  compares  it  with  the  REFPROP-4.0  data.  Because  of  the 
lack  of  the  published  experimental  data,  the  REFPROP-4.0  data  only  are  used  for  comparison.  The  measured  data 

show  the  deviations  from  REFPROP-4.0  values  are  within +5%. 

A  linear  curve  fit  equation  is  determined  for  the  measured  Cp  [  k.H(k%  ■  C  )l  versus  temperature  (°C)  and  is 
shown  below: 

Cp  =  1.2048  +  0.001925  ■  T  (10.1) 


VISCOSITY 

HFC-236ea  viscosity  was  measured  by  a  viscometer  with  its  accuracy  verified  with  a  number  of  fluids  as 
shown  in  Chapter  5  to  be  ±2%.  Figure  10.2  shows  the  measured  viscosity  of  liquid  HFC-236ea  versus  temperature 
over  a  temperature  range  from  -5°C  to  nearly  60°C.  The  deviations  from  REFPROP  values  in  Figure  10.3  indicate 
the  deviations  are  within  +5%. 

In  order  to  conveniently  use  the  results,  a  curve  fit  equation  expressed  as  a  polynomial  is  given  below: 

p  =  i  a.-T1  (10.2) 

i  =  0 

where  p  is  in  cp  and  T  is  in  °C.  The  curve  fit  coefficients  are  listed  in  Table  10.1. 
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Figure  10.1:  Comparison  of  measured  specific  heat  and  REFPROP  data  for  HFC-236ea. 
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TABLE  10.1:  POLYNOMIAL  COEFFICIENTS  FORHFC-236ea  VISCOSITY  (Equation.  10.2) 


second  degree  polynomial 

aQ 

5.65237 164E-01 

5.62097609E-01 

a\ 

-1.06925368E-02 

-8.41619726E-03 

a2 

2.31479004E-04 

5.20047470E-05 

°3 

-4.579 16803E-06 

- 

fl4 

3.800203 17E-08 

- 

°5 

-1.66399134E-11 

- 

DENSITY 

The  density  of  HFC-236ea  was  also  measured  by  a  densimeter.  The  accuracy  of  the  densimeter  has  been 
verified  with  a  number  of  refrigerants  as  discussed  in  Chapter  5.  The  resulting  plot  of  density  versus  temperature 
is  shown  in  Figure  10.4.  The  deviations  were  calculated  to  be  within  +1%  compared  with  REFPROP-4.0  values. 
The  linear  curve  fit  equation  for  the  measured  data  is: 

p  =  1514.30969  -  3.13807  •  T  (10.3) 

where  p  is  in  kg/m3  and  T  is  in  °C. 

THERMAL  CONDUCTIVITY 

The  thermal  conductivity  of  HFC-236ea  over  a  temperature  range  of  0°C  to  60°C  was  measured  by 
Approaches  1  and  2  in  this  study.  The  resulting  plot  is  shown  in  Figure  10.5  using  Approach  1  and  in  Figure  10.6 
using  Approach  2.  As  shown  in  these  two  figures,  the  measured  results  average  15%  higher  than  the  REFPROP- 
4.0  values.  Based  on  the  verification  data  for  CFC-114,  the  measured  results  imply  that  the  REFPROP  data  are 
low  for  HFC-236ea  for  the  measured  temperature  range.  A  linear  curve  fit  equation  is  given  for  easy  use  below: 


0.097013  -  0.000277  •  T 

Approach  1  (Dittus-Boelter  CF) 

(10.4) 

0.094944  -  0.000285  *  T 

Approach  1  (Petukhov-Popov  CF) 

(10.5) 

0.093838  -  0.000282  •  T 

Approach  1  (Gnielinski  CF) 

(10.6) 

0.0959604  -  0.000268  •  T 

Approach  2 

(10.7) 

where  k  is  in  \VI{nt  •  K )  and  T is  in  °C. 
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Figure  10.4:  Comparison  of  measured  density  and  REFPROP  data  for  HFC-236ea. 
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Figure  10.5:  HFC-236ea  thermal  conductivity  by  Approach  1 . 
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Figure  10.6:  HFC-236ea  thermal  conductivity  by  Approach  2. 


OTHER  PROPERTIES 


In  the  current  study,  there  are  liquid  refrigerant  properties  which  were  measured  or  calculated.  These 
inriiiA»  thermal  conductivity,  viscosity,  specific  heat,  and  density.  However,  other  properties  based  on  these  four 
properties  also  can  be  calculated.  These  include  thermal  diffusivity  and  the  Prandtl  number.  Thermal  diffusivity  is 
defined  as  k  /  pCp,  and  Prandtl  number  is  \iCp  /  k.  Figures  10.7  and  10.8  show  the  thermal  diffusivity  and  Pr  of 
HFC-236ea  versus  temperature,  respectively. 

To  produce  the  values  of  thermal  diffusivity  and  Pr,  Equation  10.7  was  used  for  k  and  Equations  10.1  and 
10.3  were  used  for  Cp  and  p,  respectively.  For  viscosity,  p.  Equation  10.2  and  the  coefficients  associated  with  the 
fifth  degree  in  Table  10.1  were  used. 

SUMMARY 

In  this  Chapter,  HFC-236ea  properties  were  measured  by  the  current  facility  and  methodology  which  were 
verified  by  the  results  of  CFC-1 14  as  discussed  in  the  previous  chapter.  The  thermodynamic  properties  ofHFC- 
236ea  (density  and  specific  heat)  showed  good  accuracy  when  compared  with  REFPROP.  However,  the  transport 
properties  (viscosity  and  thermal  conductivity)  showed  larger  deviations  than  the  thermodynamic  properties.  A 
summary  table  of  the  deviations  is  shown  in  Table  10.2.  As  shown  in  this  table,  the  deviations  of  density,  specific 
heat,  and  viscosity  are  quite  reasonable,  while  the  thermal  conductivity  deviation  is  quite  large.  However,  as 
shown  in  Figures  10.5  and  10.6,  the  trends  of  measured  and  REFPROP  data  were  quite  consistent.  A  modification 
between  these  two  measures  is  necessary  to  achieve  closer  agreement. 

TABLE  10.2:  SUMMARY  OF  MAXIMUM  DEVIATION  FOR  HFC-236ea  PROPERTIES 


properties 

■  deviation* 

<0-50*0  • 

density,  p 

±1% 

specific  heat,  Cp 

+4.8% 

viscosity,  |i 

-5% 

thermal  conductivity*,  k 

+15% 

'deviation  =  [  (measured  value  -  REFPROP  value)  /  REFPROP  value  ]  100 
*for  both  Approaches  1  and  2 
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Temperature,  T,  (°C) 

Figure  10.7:  Variation  of  thermal  diffusivity  of  HFC-236ea  with  temperature. 
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Figure  10.8:  Variation  of  Prandtl  number  of  HFC-236ea  with  temperature. 


CHAPTER  11 


PROPERTIES  OF  HFC-236ea  AND  LUBRICANT  MIXTURES 

Property  characteristics  of  a  HFC-236ea  mixed  with  a  polyol  ester  (POE)  lubricant  (Icematic  Castrol  oil 
SW68)  were  investigated.  Properties  studied  included  thermal  conductivity  and  viscosity.  However,  specific  heat 
and  density  were  also  determined  in  this  study. 

Procedures  for  lubricant  injection  and  mass  fraction  sampling  were  described  in  Chapter  5.  Four  lubricant 
concentrations  were  studied  in  this  work.  Each  concentration  determination  was  taken  after  10  hours  of  refrigerant 
circulation  through  the  test  loop  following  the  injection  of  the  lubricant.  A  summary  table  for  lubricant 
concentrations  is  shown  in  Table  11.1. 

TABLE  11.1:  LIST  OF  HFC-236ea  AND  LUBRICANT  CONCENTRATION  SAMPLES 


number 

r.fof 

11*1111 

total  cell 

l|||fH 

oil  weight  ! 

tiltH®! 

net  oil  * 

f«liwe||hli 

net  all 

!•••  weight  H 

■  1 1 

net  ref. 

■  *w%ht--4 

1  <2  ■ :: 

<41  mass 
.-^percentages 

•  m  f  i 

l 

521.91 

628.60 

106.69 

523.81 

1.90 

104.79 

1.78 

2 

521.91 

627.12 

105.21 

525.52 

3.61 

101.60 

3.43 

3 

521.91 

626.58 

104.67 

527.69 

5.78 

98.89 

5.52 

4 

521.91 

625.90 

103.99 

529.43 

7.72 

96.45 

7.42 

EFFECTS  OF  LUBRICANT  CONCENTRATION  ON  THERMAL  CONDUCTIVITY  OF  MIXTURE 

Because  of  limited  lubricant  data  sources,  the  thermal  conductivity  of  the  pure  lubricant  could  not  be 
obtained.  So  far,  few  papers  have  been  published  in  this  area.  For  refrigerant/lubricant  mixtures,  there  have  also 
been  very  limited  data  available  for  use  in  real  refrigeration  systems.  In  the  current  study,  the  thermal  conductivity 
of  HFC-236ea  with  a  selected  lubricant  (Icematic  Castrol  oil  SW68)  mixture  was  determined.  In  this  study,  the 
thermal  conductivity  calculation  was  based  on  Approach  2  (Chapter  7)  because  it  showed  less  uncertainty  and  was 
easier  to  apply  (without  CF)  than  Approach  1  as  shown  Ity  previous  results. 

Figure  11.1  shows  the  HFC-236ea  and  the  lubricant  thermal  conductivity  versus  temperature  for  various 
lubricant  concentrations.  As  shown  in  this  figure,  at  lower  temperatures  the  thermal  conductivity  increases 
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Dependence  of  liquid  thermal  conductivity  of  HFC-236ea 
and  lubricant  mixture  on  temperature. 


somewhat  with  increasing  lubricant  concentration.  It  also  shows  the  lubricant  concentration's  effect  on  the  thermal 
conductivity.  At  higher  temperatures,  a  minor  effect  of  the  lubricant  concentration  on  the  mixture's  thermal  conductivity 
was  found.  A  slight  decrease  of  thermal  conductivity  of  the  mixtures  at  higher  temperatures  and  concentrations 
was  detected  The  curve  fit  equations  are  provided  later  in  this  discussion.  It  was  found  that  the  thermal 
conductivity  of  the  mixture  increased  with  lubricant  concentration  at  the  higher  concentrations  and  lower  temperatures. 

A  nearly  40%  thermal  conductivity  increase  was  found  at  -108C  and  at  the  highest  lubricant  concentration  tested 
(7.42%).  At  a  high  temperature  (60°C),  the  thermal  conductivity  of  the  mixtures  was  found  not  to  significantly 
change  with  lubricant  concentration,  at  least  for  concentrations  of  less  than  8%  as  used  in  the  current  study. 

An  alternative  expression  for  the  mixture's  thermal  conductivity  plotted  for  thermal  conductivity  versus 
lubricant  concentration  is  shown  in  Figure  11.2.  In  this  figure,  four  constant  temperature  lines  were  plotted  versus 
lubricant  concentration.  As  shown,  thermal  conductivity  increases  with  increasing  lubricant  concentrations  at 
lower  temperatures  and  decreases  as  the  temperature  increases.  This  is  because  the  thermal  conductivity  of  the 
lubricant  is  generally  higher  than  that  of  the  refrigerant,  and  thermal  conductivity  of  liquids  generally  decreases 
with  increasing  temperature. 

A  two-variable  correlation  of  temperature  (7)  and  lubricant  concentration  (Coil)  for  thermal  conductivity 
of  the  lubricant  and  HFC-236ea  mixture  is  developed  by  a  curve-fitting  scheme  as  shown  below: 

=  a0  +  aj  +  a2Coil  +  a3  P  +  aATCoil  +  o5C^;  (11.1) 

+  a6P  +  a7PCoil  +  agTC*,  +  a9C*, 

where  *  is  in  W/m  •  AT,  f  is  in  °C,  and  is  the  mass  fraction  of  oil  in  the  mixture.  The  a,  values  are  listed  in  T able  11.2. 

TABLE  1 1.2:  COEFFICIENTS  OF  THERMAL  CONDUCnVITY 
FOR  LUBRICANT  AND  HFC-236ea  MIXTURES 


1 

a0 

0.095728 

a\ 

-0.000659 

°2 

0.004704 

°3 

3.4517E-05 

O4 

•0.000218 

ai 

-0.000273 

Of 

-5.67E-07 

8.44E-07 

1.4516E-05 

°9 

-7.671E-06 
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ctivity  of  lubricant  and  HFC-236ea  mixture 
four  different  temperatures. 


A  comparison  of  the  deviation  between  experimental  data  and  predicted  values  by  the  currently  developed 
equation  (Equation  11.1)  is  shown  in  Figure  11.3.  As  shown  in  this  figure,  the  predicted  values  agree  with  the 
experimental  data  within  ±5%. 

An  individual  equation  for  the  thermal  conductivity  of  the  lubricant  and  HFC-236ea  mixture  as  a  function 
of  temperature  is  also  provided  for  various  lubricant  mass  concentrations  over  the  temperature  range  from  -10°C 
to  60°C  and  is  shown  below: 


0.0%  lubricant:  km  =  (9.596040 10£  -  02)  +  (-2.68429925E  -  04 )7  (1 1.2) 

1.78%  lubricant:  km  =  (1.073 18327E  -  01)  +  (-7.3558 1270E  -  04)7  (1 1.3) 

3.43%  lubricant:  km  =  (1.10977404E-01)  +  (-8.13782855E-04)7  (11.4) 

5.52%  lubricant:  km  =  (1.13845766E  -  01)  +  (-9.09929338E  -  04)7  (1 1.5) 

7.42%  lubricant:  km  =  (1.15142770E  -  01)  +  (-9.53352253E  -  04)7  (11.6) 


where  km  is  in  H7(m  •  K )  and  7  is  in  °C. 

EFFECTS  OF  LUBRICANT  CONCENTRATION  ON  THE  VISCOSITY  OF  MIXTURE 

A  lubricant  (oil)  usually  has  a  very  high  viscosity  compared  with  a  refrigerant.  In  general,  the  lubricant's 
viscosity  is  several  thousand  times  greater  than  most  refrigerants  at  low  temperatures.  In  general,  the  viscosity  of  a 
lubricant  changes  rapidly  with  temperature,  especially  in  the  low  temperature  range.  However,  the  viscosity  in  the 
high  temperature  range,  although  not  significantly  changing  with  temperature,  is  still  very  high  and  could  be 
several  hundred  times  that  of  the  refrigerant  viscosity.  As  a  lubricant  is  added  to  a  refrigerant,  the  viscosity  of  the 
mixture  rapidly  increases  due  to  the  high  viscosity  of  the  lubricant  added. 

In  the  present  study,  the  viscosity  of  the  lubricant  and  HFC-236ea  mixture  was  measured.  Figure  1 1 .4 
shows  the  mixture  viscosity  versus  temperature  for  various  lubricant  concentrations.  As  shown  in  this  figure,  the 
viscosity  increases  with  a  lubricant's  increasing  concentration  and  with  a  decreasing  temperature.  An  interesting 
result  was  found  that  in  the  mixture  viscosity  rapidly  increased  with  higher  lubricant  concentrations  at  lower 
temperatures.  At  low  temperatures  (below  -10°C)  and  lubricant  concentrations  under  8%,  it  was  found  that  the 
effect  of  lubricant  concentration  on  the  mixture  viscosity  was  significant.  At  high  temperature  situations,  the 
lubricant  concentration  was  shown,  not  so  obviously,  to  affect  the  mixture  viscosity  just  as  markedly  as  at  low 
temperatures.  An  alternative  plot  for  viscosity  versus  concentration  for  various  temperatures  is  shown  in  Figure  1 1.5. 

Equations  developed  for  applying  these  results  are  useful  and  they  are  provided  in  the  following 
discussion.  Viscosity  equations  for  various  lubricant  concentrations  are  listed  below: 
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Figure  11.4:  Viscosity  versus  temperature  for  lubricant  and  HFC-236ea 
mixture  under  various  lubricant  concentrations. 
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Figure  11.5:  Viscosity  versus  lubricant  mass  fraction  for  lubricant 
and  HFC-236ea  mixture  for  various  temperatures. 


0.0%  lubricant:  In  (p) 

=  -52.355901  +  5796.69002o|;+ 0.1633357-  0.0001887* 

(11.7) 

1.78%  lubricant:  In  (p) 

*  -0.457921  +  0.0063 19^- 0.0184637+  0.0000809067* 

(11.8) 

3.43%  lubricant:  In  (p) 

=  -0.365710  -  0.008252^:-  0.0189907+ 0.0000752787* 

(11.9) 

5.52%  lubricant:  In  (p) 

*  -0.265842  -  0.000275^;  -  0.0182487  +  0.00005568672 

(11.10) 

7.42%  lubricant:  In  (p) 

=  -0.190103  -  0.0 10952^;- 0.0193397+  0.0000719047* 

(11.11) 

where  7 is  in  AT  and  p  is  in  cp. 

A  two-variable  correlation,  temperature  (7)  and  lubricant  concentration  (Coj/),  for  the  viscosity  of 
lubricant  (Icematic  Castrol  oil  SW68)  and  HFC-236ea  mixture  (pm)  was  also  developed  by  the  curve  fit  scheme  as 
shown  in  Equation  1 1. 12: 

In  (pm)  =  a0  +  axT  +  a2Coil  +  a3 7*  +  a4TCoil  +  a5C*n  (11.12) 

+  a6P  +  a7T*Coil  +  asTC*t  +  a9C*u 


where  \im  is  in  cp,  7  is  the  temperature  in  °C,  and  Ca/  is  the  mass  fraction  of  lubricant  in  the  mixture.  The  o,  are  listed  in 
Table  1 1.3.  An  accuracy  test  for  using  this  equation  was  verified  within  ±2%  with  the  experimental  data  covering  a 
viscosity  range  of  0  to  1.2  cp.  Figure  11.6  shows  the  predicted  viscosity  using  Equation  11.12  versus  experimental 
data.  It  shows  that  for  the  viscosity  range  of  the  experiments  (Oto  1.2  cp),  the  values  predicted  by  using  Equation  11.12 
agree  with  the  experimental  data. 


TABLE  1 1.3:  COEFFICIENTS  OF  VISCOSITY  FOR  LUBRICANT  AND  HFC-236ea  MIXTURE 


coefficient^ 

°o 

0.567099 

ai 

-0.010817 

a2 

0.074007 

a3 

0.000168 

a4 

-0.001519 

ai 

-0.005964 

a6 

-1.4E-06 

°7 

1.1996E-05 

aS 

4.380 IE-05 

a9 

0.000302 
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Experimental  Viscosity,  (cp) 

Figure  11  8-  Predicted  viscosity  versus  experimental  data  of  the  lubricant 
(Icematic  Castrol  SW68)  and  HFC-236ea  mixture. 


OTHER  PROPERTIES  OF  LUBRICANT  AND  HFC-236ea  MIXTURE 


In  this  section,  specific  heat  and  density  of  HFC-236ea  and  Icematic  Castrol  SW68  oil  are  calculated  and 
presented  as  curves.  Figure  1 1.7  shows  the  specific  heat  of  these  mixtures.  As  shown,  the  specific  heat  is  not 
significantly  affected  by  lubricant  concentrations  except  at  high  concentrations  of  lubricant  in  the  lower 
temperature  range.  As  the  specific  heat  of  a  pure  lubricant  is  difficult  to  obtain,  it  is  assumed  that  the  specific  heat 
for  both  materials  (HFC-236ea  and  lubricant)  should  be  fairly  close.  However,  further  verification  will  be  required 

when  specific  heat  data  of  pure  lubricants  become  available. 

The  density  of  the  mixture  is  shown  in  Figure  11.8.  As  indicated,  the  density  of  the  mixture  decreases  as 
the  temperature  increases.  However,  at  a  fixed  temperature,  the  densities  of  the  mixture  with  a  low  lubricant 
concentration  are  higher  than  those  of  high  lubricant  concentration.  This  implies  that  pure  HFC_236ea  has  a 
higher  density  value  than  pure  lubricant  currently  used. 

SUMMARY 

In  this  chapter,  the  transport  properties  of  lubricant  and  HFC-236ea  mixtures  are  presented.  Lubricant 
concentration  effects  on  both  thermal  conductivity  and  viscosity  are  discussed.  However,  viscosity  is  more 
significantly  affected  by  lubricant  concentration,  especially  in  a  lower  temperature  range.  Thermal  conductivity  is 
affected  by  lubricant  concentrations  by  8  to  40%  at  -10°C  and  only  by  1  to  10%  at  40°C  for  a  lubricant 
concentration  varying  from  2  to  7.4%. 

Other  properties,  such  as  specific  heat  and  density  of  the  mixture  of  HFC-236ea  and  Icematic  Castrol 
SW68  oil,  were  also  calculated  and  discussed.  The  results  showed  that  temperature  and  lubricant  concentration 
have  limited  influence  on  these  properties. 
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Figure  1 1 .8:  Density  versus  temperature  of  lubricant  and  HFC-230ea  mixture. 
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APPENDIX  A.  THERMOPHYSICAL  PROPERTIES 


This  appendix  discusses  the  use  of  thennophysical  properties  such  as  saturated  pressure  (Pui)>  density  (p),  specific 
heat  (Cp),  viscosity  (p),  and  thermal  conductivity  (*).  Data  sources  of  these  properties  are  found  in  a  number  of 
references,  such  as  ASHRAE,  NIST  (REFPROP),  and  Chemistry  and  Physics  Handbook.  However,  for  purposes  of 
refrigerant  study,  ASHRAE  data  and  REFPROP  data  are  usually  the  preferred  sources.  It  should  be  noted  that 
REFPROP  is  a  recently  developed  package  which  uses  theoretical  prediction  methods.  Because  of  limited 
experimental  data  for  new  refrigerants,  such  as  HFC-236ea,  REFP4.0  (REFPROP  version  4.0)  was  employed  as  the 
reference  data  source  for  comparison.  Although  REFPROP  data  includes  almost  all  of  the  refrigerants  and  their 
mixtures  in  any  combination,  it  was  discovered  that  the  transport  properties  between  these  two  data  sources  were 
sometimes  quite  different  from  each  other.  Therefore,  proper  care  should  be  taken  before  using  them.  In  this  study, 
those  properties  used  for  thermal  conductivity  measurements  are  all  based  on  the  ASHRAE  data  source  [82]  while 
comparisons  between  measured  properties  and  REFPROP  [83]  were  made  for  the  fluids  of  unknown  properties. 

In  order  to  use  property  data  conveniently,  based  on  the  ASHRAE  or  REFPROP  data  source,  curve  fit 
equations  were  developed  by  using  polynomial  functions  described  by: 

5 

property  =  .  a,  '  T‘  (A.1) 

where  temperature,  T  is  in  °C  here  and  in  the  following  discussions  unless  otherwise  stated.  The  applied 
temperature  range  of  regression  for  all  properties  is  at  the  saturation  temperature,  and  the  pressure  is  also  at  the 
saturated  condition.  The  curve  fit  temperature  range,  in  °C,  is  listed  in  Table  A.l. 
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TABLE  A.1:  CURVE-FIT  TEMPERATURE  RANGE  IN  °C  FOR  TEST  REFRIGERANTS 


m 

PM 

IBlIlBBilllll 

■Mill 

WHii 

mmm 

CFC-12 

-50  - 100 

-50  - 100 

-40-80 

-70  -  68 

-60-68 

HCFC-22 

-130-96.14 

-130-96.14 

-90-70 

o 

1 

•u 

o 

00 

1 

o 

*T 

CFC-113 

-30-214 

-30  -  214 

■n 

0-175 

0-170 

CFC-114 

-40-115 

-40-115 

-40-115 

-40-115 

-40-115 

HFC-236ea 

-40-100 

-40-100 

-40-100 

-40-100 

-40-100 

blend  A 

-40-85 

-40-85 

00 

1 

© 

"f 

o 

1 

00 

u» 

-40-85 

blend  B 

© 

1 

n| 

o 

1 

o 

-40-70 

-40-70 

-40-70 

SATURATION  PRESSURE 

The  P^j  versus  temperature  curve  fit  equations  are  provided  in  the  following  sections  for  the  tested 
refrigerants.  The  unit  for  saturation  pressure,  is  MPa  for  ASHRAE  and  psia  for  REFPROP  data,  while  the 
temperature  unit  is  °C. 

PMt  for  Calibration  Refrigerants  -  ASHRAE 

Table  A.2  lists  the  curve  fit  coefficients  of  the  saturation  pressures  of  refrigerants  based  on  the  ASHRAE 
data.  The  form  of  the  equation  for  each  curve  fit  is  a  polynomial  function  as  described  by  Equation  A.  1.  These 
data  were  used  to  verify  the  subcooled  conditions  and  to  verify  that  all  operating  conditions  were  in  the  liquid  phase. 


TABLE  A.2:  SATURATION  PRESSURE  CURVE-FIT  COEFFICIENTS 
FOR  HCFC-22,  CFC-12,  CFC-113,  AND  CFC-114 


.  polynomial 

CFC-1J  ! 

1  CFC-113  i 

CFC-114 

a®,  MPa 

5.00005E-01 

2.91020E-01 

1.58554E-02 

8.76786E-02 

aj  ,MPa/°C 

1.61909E-02 

1.00845E-02 

6.77690E-04 

3.47903E-03 

aj.MPa/rC)2 

1.95017E-04 

1.49920E-04 

1.25848E-05 

5.45875E-05 

a3 ,  MPa/(°C)3 

1.03242E-06 

5.6405  IE-07 

2.282 13E-07 

3.74883E-07 

a4 ,  MPa/(°C)4 

2.23206E-09 

5.91877E-09 

1.47081E-10 

4.43447E-10 

a5 ,  MPa/(°C)5 

1.71047E-12 

5.79095E-11 

1.60668E-12 

1.88150E-12 

determinant 

9.99996E-01 

9.97749E-01 

9.99998E-01 

9.99999E-01 
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P„,  for  HFC-236ea  -  REFPROP 


Table  A.3  lists  the  saturation  pressure  (psia)  of  refrigerants  based  on  the  REFPROP  data. 
TABLE  A.3:  SATURATION  PRESSURE  CURVE-FIT  COEFFICIENTS  FOR  HFC-236ea 


polynomial  ; 
lllll  coefficient^ 

|  HFC«236ea  '  ' 

ao ,  psia 

1.13349352E+01 

a} ,  psia/°C 

4.8487079 IE-01 

*2 »  psia/(°C)2 

8.3043659  IE-03 

a3 ,  psia/(°C)3 

6.785 101 14E-05 

a4 ,  psia/(°C)4 

2.04624 143E-07 

a5 ,  psia/(°C)5 

-2. 1903540  IE-10 

determinant 

1.00000000E+00 

LIQUID  DENSITY 

The  liquid  density  versus  temperature  curve  fit  equations  are  provided  in  the  following  sections  for  the 
tested  refrigerants.  The  unit  for  density,  p,  is  (k g/mi3)  while  the  unit  for  temperature  is  °C. 

Density  for  Calibration  Refrigerants  -  ASHRAE 

Because  the  type  of  viscometer  used  in  this  study  measures  the  product  of  density  (g/cm 3)  and  dynamic 
viscosity  (cp),  therefore,  density  is  a  property  needed  to  be  measured  in  the  present  study  in  order  to  obtain  the 
viscosity.  For  the  calibration  of  refrigerants  of  known  properties,  the  ASHRAE  data  for  density  was  used  to  verify 
both  measured  density  and  viscosity.  Table  A.4  lists  the  density  curve  fit  equation  coefficients  at  the  saturated 
liquid  state  based  on  ASHRAE  data. 
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TABLE  A.4:  LIQUID  DENSITY  CURVE-FIT  COEFFICIENTS  FOR  HCFC-22,  CFC-12,  CFC-1 13,  CFC-1 14 


erant 

jf  j  eeeflfcitmt | § '  j| 

HI 

a« ,  kg/m3 

1.27397E+03 

1.39594E+03 

1.613 16E+03 

1.50572E+03 

MSiHM 1 

-3.39602E+00 

-3.22976E+00 

-2.00082E+00 

-1.76217E+00 

a 2 ,  kg/m3  •  (°C)2 

5.30874E-03 

-7.53629E-03 

1.26301E-02 

-1.78123E-02 

a3 ,  kg/m3  •  (°C)3 

1.21421E-05 

-7.65376E-06 

-4.930 12E-04 

-5.05889E-05 

a, ,  kg/m 3  *  (°C)4 

-2.76909E-06 

3.27685E-07 

4.29266E-06 

1.49907E-06 

MUEBSSai 

-1.82476E-08 

-1.14387E-08 

9.94815E-01 

-8.35448E-09 

determinant 

9.93922E-01 

9.99988E-01 

9.94815E-01 

9.99996E-01 

Liquid  Density  for  HFC-236ea  -  REFPROP 

The  curve  fitting  coefficients  of  liquid  density  based  on  REFPROP  are  listed  in  Table  A.  5. 
TABLE  A.5:  LIQUID  DENSITY  CURVE-FIT  COEFFICIENTS  FOR  HFC-236ea 


§|  polynomial 

coeIDclcii^^®i 

:HFC-236ea 

a Q,kg/m3 

1.50280E+03 

a, ,  kg/m3  •  °C 

-2.84734E+00 

a2  ,  kg/m3  •  (°C)2 

-5.44754E-03 

a3  ,  kg/m3  •  (°C)3 

-2.65424E-05 

a4  ,  kg/m 3  •  (°C)4 

1.76196E-08 

a5  ,  kg/m3  •  (°C)5 

-1.89195E-09 

determinant 

9.99996E-01 

LIQUID  SPECIFIC  HEAT 


The  liquid  specific  heat  versus  temperature  curve  fit  equations  are  provided  in  the  following  sections  for 
the  tested  refrigerants.  The  unit  for  specific  heat,  Cp,  is  U/ikg  •  K>,  while  the  unit  for  temperature  is  °C. 
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Specific  Heat  for  Calibration  Refrigerants  -  ASHRAE 

Table  A.6  lists  the  specific  heat  curve  fit  equation  coefficients  for  the  saturated  liquid  state  of  HCFC-22, 
CFC-12,  CFC-113,  and  CFC-114  based  on  the  ASHRAE  data. 

TABLE  A.6:  LIQUID  SPECIFIC  HEAT  CURVE-FIT  COEFFICIENTS 
FOR  HCFC-22,  CFC-12,  CFC-113,  CFC-114 


jxdy&omiil 

refrigerant  1 

coefficient 

ao ,  kJ/(kg  •  K) 

1.17392E+00 

9.27994E-01 

9.21305E-01 

9.3086  IE-01 

aj ,  kJ/(kg  •  K)-°C 

2.72938E-03 

1.87657E-03 

1.70157E-03 

1.13601E-03 

a 2,kJ/(kg’K)‘CC? 

1.46720E-05 

1.52956E-05 

-1.40780E-05 

-3.84778E-06 

a3 ,  kJIQcg •  K)  *  (°C)3 

1.18154E-07 

-1.12761E-08 

9.15992E-08 

1.96439E-07 

a4 ,  kJ/Qcg  •  K)  •  (°C)4 

6.23786E-09 

-5.42684E-10 

6.44030E-1 1 

8.34917E-11 

a}  ,  kJ/(kg  •  K)  •  (°C)5 

5.16438E-11 

4.6642  IE-11 

-1.804 50E-13 

-4.94558E-13 

determinant 

9.99485E-01 

9.99954E-01 

9.99994E-01 

9.99993E-01 

Specific  Heat  for  HFC-236ea  -  REFPROP 


The  coefficients  of  specific  heat  curve  fit  based  on  REFPROP  are  listed  in  Table  A.7. 

TABLE  A.7:  LIQUID  SPECIFIC  HEAT  CURVE-FIT  COEFFICIENTS  FOR  HFC-236ea 


§|§|||||  polynomial  pllllll 
coefficient 

HFC-236ea 

a© ,  kJ/(kg  ’K) 

1.15099E+03 

a, ,  kJ/(kg  'K)‘°C 

2.45936E+00 

aj ,  kJKkg  •  K)  •  CQ? 

4.03757E-03 

a3 ,  kJ/(kg  ‘  K)  •  ("C)3 

1.64306E-05 

i4,kJ/(kg‘K)-(0C)4 

-7.58220E-07 

a3 ,  kJ/(kg  •  K)  •  (°C)5 

1.53071E-08 

determinant 

9.99988E-01 

liquid  viscosity 


The  liquid  viscosity  versus  temperature  curve  fit  equations  are  provided  in  the  following  sections  for  the 
tfstfd  refrigerants.  The  units  for  viscosity,  p,  are  Pa  *  s  while  the  unit  for  temperature  is  °C. 

Viscosity  for  Calibration  Refrigerants  -  ASHRAE 

Table  A.8  lists  the  dynamic  viscosity  (in  Pa  •  s)  curve  fit  equation  coefficients  at  the  saturated  liquid  state 
based  on  ASHRAE  data. 

TABLE  A.8:  LIQUID  VISCOSITY  CURVE  FIT  COEFFICIENTS  FOR  HCFC-22, 

CFC-12,  CFC-113,  AND  CFC-114 


polynomial 

...  ||  - 

refrieerant  111 

■ran 

ao  ,Pa’s 

2.10126E-04 

2.57691E-04 

9.545 13E-04 

4.91168E-04 

&y,Pa'  sf°C 

-2.278 12E-06 

-2.77653E-06 

-1.31618E-05 

-6.55940E-06 

a 2,Pa-s/(°C? 

1.14319E-08 

1.89625E-08 

1.11619E-07 

4.48919E-08 

Zy,  Pa'  s/(°C)3 

3.80621E-12 

-8.65530E-11 

-6.34269E-10 

-1.26949E-10 

a4 ,  Pa •  s/(°C)4 

-7.75142E-13 

-2.19818E-14 

2.15999E-12 

-4.40069E-13 

a 5,Pa-s/(°C)5 

7.72755E-15 

2.9647  IE-15 

-3.2360  IE-15 

3.18456E-15 

determinant 

9.99999E-01 

9.99999E-01 

9.99999E-01 

9.99999E-01 

Liquid  Viscosity  for  HFC-236ea  -  REFPROP 

The  coefficients  of  dynamic  viscosity  curve  fit  equation  based  on  REFPROP  are  listed  in  Table  A.9. 
TABLE  A.9:  LIQUID  VISCOSITY  CURVE-FIT  COEFFICIENTS  FOR  HFC-236ea 


coefficient  ill 

BfC-236e | 

ao ,  kJ/(kg '  K) 

5.83657E-04 

aj ,  kJI(kg '  K)‘°C 

-9.32850E-06 

Zi'kJ/ikg-V-CQ? 

1.03792E-07 

ZytkJI{kg'K)'(°Cf 

-1.08120E-09 

a4 ,  kJ/(kg '  K) '  (°C)4 

8.20849E-12 

a}  ,  kJ/(kg  •  K) '  (°C)5 

-2.863 13E-14 

determinant 

9.99997E-01 
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liquid  thermal  conductivity 


The  liquid  thermal  conductivity  versus  temperature  curve  fit  equations  are  provided  in  the  following 
sections  for  the  tested  refrigerants.  The  units  for  thermal  conductivity,  k,  are  Wl(m  *  K)  while  the  unit  for 
temperature  is  °C. 

Thermal  Conductivity  for  Calibration  Refrigerants  -  ASHRAE 

Thermal  conductivity  quoted  for  the  calibration  refrigerants  was  based  on  ASHRAE  data.  Table  A.  10  lists 
the  thermal  conductivity  (in  lV/(m  •  K))  curve  fit  equation  coefficients  at  the  saturated  liquid  state. 

TABLE  A.  10:  LIQUID  THERMAL  CONDUCTIVITY  CURVE-FIT  COEFFICIENTS 
FOR  HCFC-22,  CFC-12,  CFC-113,  AND  CFC-114 


polynomial 

erant 

|  coefficient 

■  ,  HCFC22 . 

a0 ,  W[(jn •  K) 

9.61623E-02 

7.71084E-02 

8.19249E-02 

7.04566E-02 

a, ,  Wl{m  -K)-°C 

-4.245 13E-04 

-3.66836E-04 

-1.74554E-04 

-2.41089E-04 

2i2 ,  W/(m  •  K)  *  (°C)2 

4.25538E-07 

1.77984E-07 

-6.073 17E-08 

2.54213E-07  . 

a3  ,  WKm  •  K)  •  (°C)3 

4.88343E-10 

-4.95688E-10 

1.01692E-09 

-9.79673E-10 

a4  ,  WHjn  •  K)  •  (°C)4 

-2.09592E-I1 

6.91404E-13 

-6.47263E-12 

-2.41232E-11 

a5  ,  Wl(m  ■  K)  ■  (°C)5 

-4.18749E-13 

1.10221E-13 

1.55361E-14 

1.19913E-13 

determinant 

9.99996E-01 

9.99996E-01 

9.99989E-01 

9.99994E-01 

Thermal  Conductivity  for  HFC-236ea  -  REFPROP 


The  coefficients  of  thermal  conductivity  curve  fit  equation  based  on  REFPROP  are  shown  in  Table  A.  1 1 . 
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TABLE  A.ll:  LIQUID  THERMAL  CONDUCTIVITY  CURVE-FIT  COEFFICIENTS  FOR  HFC-236ea 


coefficient  fM  . 

Si  ^fi®C436ea:>;i 

ao ,W/(m-K) 

7.98508E-02 

at ,  W/(m  ‘K)’°C 

-3.22129E-04 

,  W/(m  *  K) '  (°C)2 

-5.87303E-07 

a3 ,  WHm  •  K)  •  fC)3 

8.10775E-09 

a4 ,  Wl(m  •  K)  •  (°C)4 

2.08164E-1I 

-4.47036E-13 

determinant 

9.99743E-01 
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APPENDIX  B.  UNCERTAINTY  ANALYSIS  FOR  THERMAL  CONDUCTIVITY  MEASUREMENTS 


This  appendix  analyzes  the  uncertainties  of  the  thermal  conductivity  measurement  for  the  methods  of 
Approach  1:  the  Nusselt  number  method,  and  Approach  2:  the  Prandtl  number  method. 


UNCERTAINTY  ANALYSIS  FOR  APPROACH  I 

Three  different  single-phase  correlations  employed  in  the  thermal  conductivity  calculations  will  be  analyzed 
for  their  The  uncertainties  of  the  measured  thermal  conductivity  arise  due  to  both  sensor  (or  instrument) 

uncertainties  and  the  Nusselt  number  correlation  uncertainties.  However,  in  this  analysis,  only  the  sensor  uncertainties 
are  included  for  the  uncertainty  calculation,  while  the  Nusselt  number  uncertainty  is  discussed  in  this  current  analysis. 


Uncertainty  Analysis  for  Thermal  Conductivity  Calculated  from  the  Dittus-Boelter  Correlation 


The  thermal  conductivity  obtained  from  the  Dittus-Boelter  correlation  is: 


u  -  r  J£- »,-08  ,|-0.4  r„4A 11«7 
kd  -  1-0023**  ^  Cp  ■» 


=  FlhD.Re&^Cp)  =  F(h,D, m,  p,  Cp) 


(B.l) 

(B.2) 


The  thermal  conductivity  presented  here  is  a  function  of  A,  D,  Re^,  p,  and  Cp.  Therefore,  the  uncertainty  of  the 
sum  of  ukd  I  kd  can  be  expressed  as  follows. 


U, 


/  4  **Pv 
v"3  Rer  > 


-  U  «  Uc 

'  3  3  r> 


(B.3) 


'  SSkd  (B.4) 

where  Uj,  /  h,  and  uRep)  /  Re^  are  discussed  in  the  following  sections.  Then,  the  uncertainty  ratio  UrA  (=  ukd  /  kd) 
is  obtained  by  taking  square  root  of  the  sum  of  squares  of  each  uncertainty  source,  ^kd  That  is. 
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The  thermal  conductivity  from  the  Petukhov  and  Popov  correlation  is  found  to  be. 


kp  =  pQ? 


[ 


mmuCpRen  MZ/ZvW 
101.6 hD  "  12.7  v  8  * 


(B.6) 


=  F(  p,  CpJ, Rep,  h,Z>)  “  F(  p,  CpJ,  m,h,D)  (B.7) 


In  this  equation,  kp  is  function  of  p,  CpJ,  Re^  h,  and  D.  Therefore,  the  uncertainty  of  sum  of  square  of  \  /  kp 
can  be  expressed  as  follows: 


if? =  <a£**)  ^scPX)  kP} 

■AUxt+fStEzr 

v  ar  1 k-J 


=  sst 


(B.8) 

(B.9) 


where  each  item  in  the  above  equation  is  derived  as  follows. 


_  5l_2  r  J*£lL_L21  (£yin  +  i  T1 .  [  frcPRen  5*1 
ap  kp  p  2 1 J01.6SD  12.7  J  L101.67iDm  kPJ 

dk^uez  ucvj  \  mm  m  (Lx™ + ! I1 .  r  frcPRe,n—  S 

acp  Cp  2  L 101 67f2>  12.7  v8'  J  L101.6SD(/78)^  ^ 


3  r  (8/)1/2  107  rrvl/2  .  .  V» .  [  t'Wh&Ren  .  81/2  (1.07)  M  -I  S 

+lJ  L  203.2SD  25.4  '  J/ 

3kp  UX*D  3  [•  f8Am  1.07  ,  /•  yl/a  1-1  r  Sm/'a»CpRen  1 

V.  E  ^  I _ --*r- _  117  \  8  /  J  L  i  A1  /Tn  J  R& n 


(B.10) 


(B.ll) 


(B.12) 


(B.13) 


^rEH-  3  r  (8rt1/2  L07  r  l  yl/2  T-l  r  _  8»«/«pCpRen  1  jfr 
afc  kp  "2l.joi.6SD  12.7  '8'  J  L  101.6SD  J  h 


CB.14) 


(B.15) 


U°  _  _  3  T  (8y)1/2  -  107  (  W2  4  1  T1  r  -  8l/2^1/2^gn  1  UD 

dD  kp  2  k  101 .6HD  12  7  ' 8 '  ^  101.6 HD  *  D 

Then,  the  uncertainty  of  ufpk  (=  ukp  tkp)  is  obtained  by  taking  square  root  of  the  sum  of  squares  of  each 
uncertainty  source,  SSk^.  That  is 

^L=(SS>cp)in  0.16) 

P 

Uncertainty  Analysis  for  Thermal  Conductivity  Calculated  from  the  Gnielinski  Correlation 


The  thermal  conductivity  from  the  Gnielinski  correlation  is  found  to  be: 

m.6hD 


*  =  mq,  [  ™m»Cp  (Rej>:  10°0)  -  ^({yw  4 1  Ya 

g  r  r  L  im  «fcn  12.7  v »  '  J 


F(  m  Cp,f, Rep,  h,D)  *  F( n,  Cp,f,  m,h,D) 


(B.  17) 
(B.  18) 


In  this  equation,  kg  is  function  of  m  Cp,f  Rep,  h,  and  D.  Therefore,  the  uncertainty  of  the  sum  of  the  squares  of 
Ukg  /  kg  can  be  expressed  as  follows: 


£4 


dkB  U„ 


3*. 


.8k. 


g 


(^r> 

dh  Kg 


scp  kk 


£/c  St  c/, 

T£)2  +  (  3/  *  )2  +  ( 


3tg 
dRep 


V, 


Rer 


CB.  19) 


(B.20) 


where  each  term  in  the  above  equation  is  derived  as  follows: 


8kg  Un  Un  3  r  (.V)inmCp  (Ren  -  1000)  1  ,  fym  .  l-l  \foCp  (Ren  - 1000)  i  Uu 
Btkkg-p-n  JOUSD  +1J  L  101.6^8)^  if 


(B.21) 


^Cp  =  ^Cp  _  3  r  (8/)1/2mQ>  (fen  -  1000)  _  1  /Xyl«  .  1*» .  r/nQ>  (flen  -  1000)  i  ^Cp 
3C/>  Q>  2  l  101.6AD  12.7  ^8^  J  L  101.6AD  (/78) 1/2  *"  CP 

dkzUf  Wm™Cp  (Ren  -  lOOO)  1  //•via  1-1 

dfkg  ~ll  101.6 hD  “12.7^  8^  J 


(B.22) 


(B.23) 
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_s KU*‘d_  3  r  (8J)1  V/>  (Reg  -  1000)  JL  (£ym  +  i  I'1  •  f  1  (B  24) 

mtD  kg  2  L  101.6AD  127  ^  J  L  101 .6hD  1  ReD 

dj^u-h  arCs^Vp^p-iooo)  J_rA-i*  -i-i .  r .  miHA  <&d  z  im  1 S  ^ 

MTg=-2l  l0l6-hD  -UJU)  +1J  L  101.6AD  Ja 

a^C^_  3  r  (8»1/2^Q>  (Ren  -  1000)  V1  T  (Reg  -  1000)  1  ^  _ 

dD  kg  ~  2*-  101.6AD  12-7  8  101.6AZ? 

Then,  the  uncertainty  ratio  of  ^rgjk  (=  ^  /  *^ )  is  obtained  by  taking  square  root  of  sum  of  square  of  each 

uncertainty  source,  That  is 

^,(SSis)‘n  (B.27) 


Uncertainty  Analysis  for  Heat  and  Flow  Measurements 

The  uncertainty  calculation  in  heat  and  flow  measurements  would  contribute  to  those  uncertainty 
quantities  of  the  thermal  conductivity  uncertainties  mentioned  in  the  above  three  sections.  In  this  section,  analysis 
of  these  heat  and  flow  uncertainties  will  be  discussed.  The  uncertainties  of  the  heat  transfer  measurements  include 
those  for  the  heat  transfer  rate,  Q,  and  heat  transfer  coefficient,  h,  while  the  uncertainties  of  fluid  flow 
measurements  include  those  from  the  calculations  of  the  Reynolds  number,  Re^,  and  friction  factor,  f 

Uncertainty  Analysis  of  Heat  Transfer  Measurement 


The  Uncertainty  of  Q— 


The  net  heat  transfer  rate  is  calculated  from  the  following  equation: 


-  hA{Tw-Tf) 

(B.28) 

=  mCp  (T0-Tj) 

(B.29) 

=  F(m,  Cp,  Ta,  Tt ) 

(B.30) 

Therefore,  the  uncertainty  of  the  sum  of  the  squares  of  Uq  /  Q  is  given  as: 
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(B.31) 


(f  )2  =  (^)2  +  (^)2  +  (f)2  +  (f)2 


The  resulting  Uq  /  Q  is  given  by: 


(B.32) 


(B.33) 


The  Uncertainty  of  h- 


The  average  heat  transfer  coefficient  is  obtained  from  the  following  equation: 


A(Tw~Tf) 


(B.34) 


nDL(Tw~Tf) 


(B.35) 


=  F(  Q,D,  L,  Ty) 


(B.36) 


The  uncertainty  of  the  sum  of  the  squares  of  I  h  is  given  by: 


&r -(%r+(-%r+(-%r+ 


uT  u7r 

(-rry)2H-=rV): 

lw  Iw  If 


(B.37) 


-3% 


(B.38) 


Then,  the  uncertainty  of  Uj;  /  h  is  given  by: 


T-tatf 


(B.39) 
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Or,  the  resulting  uncertainty  of  Uj-  //is  derived  as 


1.58 


URer 


f  0.79 In  (ReD)  -  1.64  ReD 


(B.49) 


UNCERTAINTY  ANALYSIS  FOR  APPROACH  2 

The  Approach  2  for  thermal  conductivity  is  based  on  the  Prandtl  number,  Pr,  measurement.  The  Pr  is 
expressed  in  Equation  6.34,  which  means 


Pr  =  F(D,L,ReD,  AT*) 


(B.50) 


Because  Pr  —  \iCp  /  k,k  is  function  of  p,  Cp,  and  Pr.  In  other  words, 

k  =  ^  =  F(vi,Cp,Pr ) 


(B.51) 


The  uncertainty  of  k  is  expressed  as  follows. 


.  iy*3v>+/’  + 

*  "  LvSm,  *  ;  +\ecP  k  >  +{<dPr  k  >  J 


(B.52) 


After  taking  the  partial  derivative  of  each  term  and  rearranging,  this  equation  can  be  written  as 


ui  -  r 


t/.. 


U, 


UD 


(B.53) 


Uncertainty  of  Pr 


The  measured  Pr  is  a  function  of  D,  L,  ReD,  and  AT*.  Therefore,  the  uncertainty  of  Pr  is  expressed  as 

=  T(^.^y  +  (^.Hky  +  (JZ!L.¥x£n)2  +  (J£!L.E*IlyV'2  &54) 

Pr  l-V5Z)  D  )  +  \  dL  L  t  +  V  dReD  ReD  >  ^  DAT*  AT*  >  J  (  ) 

where  UReD  /  ReD  is  found  from  Equation  B.42  and  /  AT*  is  given  in  the  following  section. 
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Uncertainty  of  AT* 


The  AT*  is  defined  as  (Ta  -  Ty)  /  (Tw  -  Tj).  The  uncertainty  of  AT*  is  expressed  as  the  following  equation: 


APPENDIX  C.  SOME  PHYSICAL  PROPERTIES  OF  TEST  REFRIGERANTS 


TABLE  C.  1 :  SOME  PHYSICAL  PROPERTIES  OF  TEST  REFRIGERANTS 


KHfll 

1*  'fatiwvai^ 

CFC-12 

120.91 

0.1814 

-29.76 

111.80 

4113 

HCFC-22 

86.47 

0.2211 

-40.86 

96.15 

4974 

HFC-32 

52.02 

0.2671 

-51.75 

78.21 

5830 

CFC-113 

187.38 

0.2563 

47.65 

214.35 

3437 

CFC-114 

170.92 

0.2511 

3.65 

145.65 

3259 

CFC-125 

120.03 

0.2953 

-48.57 

66.18 

3571 

HFC-134a 

102.03 

0.3235 

-26.15 

101.15 

4067 

HFC-143a 

84.04 

0.2566 

-47.35 

73.10 

3811 

HFC-152a 

66.05 

0.2573 

-24.65 

113.55 

4492 

HFC-236ea 

152.05 

0.3860 

6.50 

139.29 

3533 

In  Table  C.l,  the  acentric  factor  is  a  component  constant,  which  is  defined  as: 

<D  =  -  logP^  (at  Tr  -  0.7)  -  1.000 


(C.l) 


To  obtain  values  of  to,  the  reduced  vapor  pressure  (P vpr  Pvp  /  Pc)  at  rr  =  77  Tc  =  0.7  is  required.  One  of  the 
calculation  methods  was  developed  from  the  Lee-Kesler  vapor  pressure  relation  which  was 

co  =  J  (C.2) 

a  =  -In  (Pe)  -  5.92714  +  6.096489'1  +  1.28862/n  (9)  -  0. 16934796  (C.3) 

p  =  -  15.2518  -  15.68759*1  -  13.4721/n  (9)  +  0.4357796  (C.4) 

9  =  Tb/Tc  (C.5) 


In  the  equation  for  a,  the  critical  pressure,  Pc ,  must  be  expressed  in  atmospheres. 
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APPENDIX  D.  EQUIPMENT  AND  INSTRUMENTATION  SPECIFICATIONS 


The  manufacturer  and  model  number  of  each  component  in  the  test  apparatus  are  shown  in  Table  D.l. 
Electronic  instrumentation  is  listed  in  Table  D.2. 


TABLE  D.l:  COMPONENTS  OF  THE  TEST  RIG 


£ :  component 

refrigerant  pump 

Wanner  Eng. 

D-10 

motor 

GE 

1  hp,  dc 

58PB56SAA42A 

motor  controller 

Dart  Inc. 

variable  current 

250G 

filter  drier 

Sporlan 

desiccant,  262  cm3 

C-164-S 

auxiliary  heater 

Omega 

heat  tape,  627  W,  120V, 

1.27  cm  W  x  244  cm  L 

STH05 1-080 

teat-section  heater 

Omega 

heat  tape,  2  @  627  W 

120V,  1.27  cm  W  x  244  cm  L 

STH05 1-080 

mixer 

Omega 

FMX  7301 

sight  glasses 

Sporlan 

1.27  cm,  solder  fit 

- 

condenser 

Doucete 

coaxial  coil 

CX-H050 

accumulator 

Oil- Air  Ind. 

Bladder  (butyl,  Buna-N) 

1-100-1 

charging  valve 

packed  angle 

9271 

ores,  relief  valve 

SS-4R3A-A 

chiller 

Heat-X  Inc. 

CFC-12,  3  HP 

PC-300 

line  conditioner 

Tripp  Lite  Co. 

spike,  surge  suppresser 

LC1800 

1MM 

pneumatic  cylinder  type 

UDR-SS-32-6 

sample  cylinder 

Swagelok  Co. 

closed-ends  cylinder 

304L-HDF4-75 
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TABLE  D.2:  DATA  ACQUISITION  AND  INSTRUMENTATION  COMPONENTS 


■■■HI 

model 

Zenith 

386  PC 

Z-386-20 

multimeter 

HP  Inc. 

digital  multimeter 

3457A 

scanner 

HP  Inc. 

digital  control  unit 

3488A 

meter  board 

HP  Inc. 

armature  relay 
multiplexer,  10  channel 

44491A 

scanner  board 

Hewlett-Packard 

multiplexer,  10  channel 

44470A 

interface  card 

National  Instr. 

776207-01 

MiH 

bare  bead,  teflon  shield 

T-type 

m 

electric  ice  point 

CJ-T 

RID  probe 

Hy-Cal  Eng. 

100  ohm  platinum 
o.oo385/n/n/°c 

RTS-36-T-100-A 

-5-2-36-X55 

RTD  transmitter 

tSBSSMB 

2-wire,  4-20  mA,  100  ohm 

CT-801-A-S-X5 

watts  transducer 

Jemtec  Co. 

2- element,  3-wire, 

3- phase,  accuracy:  ±0.2% 

XL31K5A2 

mass  flow  meter 

Micro  Motion 

0-1800  kg/hr 
accuracy:  ±0.15% 

RFT9739 

densitometer 

Micro  Motion 

RFT9739 

viscometer 

Nametre  Co. 

■Hi 

1810-LV 

Setra  Co. 

0-250  psia,±0.11%FS 

C280E 

Setra  Co. 

0-1  psi,  ±0.15%FS 

C228-1 

voltage  regulator 

120  V,  1.4  kVa 

1010 
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APPENDIX  E.  SAS  REGRESSION  OUTPUTS 


This  appendix  lists  the  output  of  the  SAS  program  for  three  correlation  regressions  in  the  Nusselt  number 
method  (Approach  1)  and  two  variable  (ReD  and  63*)  regressions  in  the  Prandtl  number  method  (Approach  2)  for 
three  calibration  refrigerants  (HCFC-22,  CFC-12,  and  CFC-113)  and  four  calibration  refrigerants  (HCFC-22, 
CFC-12,  CFC-1 13,  and  CFC-1 14).  The  SAS  program  (version  6.07)  was  accessed  from  the  campus  Vincent 
computer  woric  station  network.  The  first  section  in  this  appendix  is  the  SAS  output  for  three  correlations 
discussed  in  Approach  1  whereas  the  second  section  is  the  SAS  output  for  In  (. Pr )  regression  with  variables  of  In 
(ReD)  and  AT*  discussed  in  Approach  2.  Both  sections  contain  three-  and  four-calibration  refrigerants  cases. 
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NUSSELT  NUMBER  REGRESSION  IN  APPROACH  1 


THREE  CALIBRATION  REFRIGERANTS 


The  SAS  System  1 

22:35  Saturday,  September  10. 1994 


Model:  MODEL  1 
Dependent  Variable:  NUE 

Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

2 

31.00300 

15.50150 

6419.126 

0.0001 

Error 

120 

0.28979 

0.00241 

C  Total 

122 

31.29279 

Root  MSE 

0.04914 

R-square 

0.9907 

Dep 

Mean 

5.53266 

Adj  R-sq 

0.9906 

C.V. 

0.88821 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  |T| 

INTERCEP 

1 

-4.997045 

0.10783630 

-46.339 

0.0001 

LNRE 

1 

0.903770 

0.00876542 

103.106 

0.0001 

LNPR 

1 

0.502199 

0.01131934 

44.366 

0.0001 

The  SAS  System 

22:35  Saturday,  September 

Model:  MODEL  2 

Dependent  Variable: 

NUD 

Analysis  of  Variance 

Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

2 

25.5452 2 

12.77261 

• 

• 

Error 

120 

0 

0 

C  Total 

122 

25.54522 
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Root  MSE  0.00000  R-square  1.0000 

DepMean  5.48384  AdjR-sq  1.0000 

C.V.  0.00000 

Parameter  Estimates 
Parameter  Standard  T  for  HO: 

Variable  DF  Estimate  Error  Parameter  =0  Prob>|T| 

INTERCEP  1  -3.772261  0.00000000 

UMRE  1  0.800000  0.00000000 

LNPR  1  0.400000  0.00000000 

The  SAS  System  3 

22:35  Saturday,  September  10, 1994 

Model:  MODEL  3 
Dependent  Variable:  NUP 

Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square  , 

F  Value 

Prob>F 

Model 

2 

26.09945 

13.04973 

189914.474 

0.0001 

Error 

120 

0.00825 

0.00007 

C  Total 

122 

26.10770 

Root  MSE 

0;00829 

R-square 

0.9997 

Dep  Mean 

5.57329 

Adj  R-sq 

0.9997 

C.V. 

0.14873 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  |T| 

INTERCEP 

1 

-4.059037 

0.01819021 

-223.144 

0.0001 

LNRE 

1 

0.827298 

0.00147858 

559.521 

0.0001 

LNPR 

1 

0.455283 

0.00190939 

238.445 

0.0001 

The  SAS  System  4 

22:35  Saturday,  September  10, 1994 

Model:  MODEL  4 
Dependent  Variable:  NUG 

Analysis  of  Variance 
Sum  of  Mean 

Source  DF  Squares  Square  F  Value  Prob>F 
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Model 

2 

29.39758 

14.69879 

51780.573 

0.0001 

Error 

120 

0.03406 

0.00028 

C  Total 

122 

29.43164 

Root  MSE 

0.01685 

R-square 

0.9988 

Dep  Mean 

5.58170 

Adj  R-sq 

0.9988 

C.V. 

0.30185 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  > 

INTERCEP 

1 

-4.415096 

0.03697204 

-119.417 

0.0001 

LNRE 

1 

0.862783 

0.00300526 

287.091 

0.0001 

LNPR 

1 

0.441244 

0.00388087 

113.697 

0.0001 
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FOUR  CALIBRATION  REFRIGERANTS 

The  SAS  System  1 

23:37  Saturday.  September  10. 1994 


Model:  MODEL  1 
Dependent  Variable:  NUE 

Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

2 

42.89045 

21.44522 

10119.359 

0.0001 

Error 

166 

0.35179 

0.00212 

C  Total 

168 

43.24224 

Root  MSE 

0.04604 

R-square 

0.9919 

Dep  Mean 

5.48795 

Adj  R*sq 

0.9918 

C.V. 

0.83884 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  |T| 

INTERCEP 

1 

-5.015673 

0.08491757 

-59.065 

0.0001 

LNRE 

1 

0.905027 

0.00687972 

131.550 

0.0001 

LNPR 

1 

0.507640 

0.00970171 

52.325 

0.0001 

The  SAS  System 

23:37  Saturday,  September  10, 1994 


Model:  MODEL  2 
Dependent  Variable:  NUD 

Analysis  of  Variance 
Sum  of  Mean 

Source  DF  Squares  Square  F  Value  Prob>F 

Model  2  35.03234  17.51617 

Error  166  0  0 

C  Total  168  35.03234 

Root  MSE  0.00000  R-square  1.0000 

DepMean  5.43887  AdjR-sq  1.0000 

C.V.  0.00000 

Parameter  Estimates 
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Parameter 

Standard 

Variable 

DF 

Estimate 

Error 

INTERCEP 

1 

-3.772261 

0.00000000 

LNRE 

1 

0.800000 

0.00000000 

LNPR 

1 

0.400000 

0.00000000 

The  SAS  System 

T  for  HO: 

Parameter  =  0  Prob  >  |T| 


23:37  Saturday,  September  10, 1994 


Model:  MODEL  3 
Dependent  Variable:  NUP 


Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

2 

35.99052 

17.99526 

167377.653 

0.0001 

Error 

166 

0.01785 

0.00011 

C  Total 

168 

36.00837 

Root 

MSE 

0.01037 

R-square 

0.9995 

Dep  Mean 

5.53387 

Adj  R-sq 

0.9995 

C.V. 

0.18737 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  | 

INTERCEP 

1 

-4.052869 

0.01912666 

-211.896 

0.0001 

LNRE 

1 

0.826843 

6.00154957 

533.594 

0.0001 

LNPR 

1 

0.457499 

0.00218519 

209.363 

0.0001 

The  SAS  System 


23:37  Saturday,  September  10, 1994 


Model:  MODEL  4  • 
Dependent  Variable:  NUG 


Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

2 

40.13652 

20.06826 

67855.395 

0.0001 

Error 

166 

0.04909 

0.00030 

C  Total 

168 

40.18562 

Root  MSE 

0.01720 

R-square 

0.9988 
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Dep  Mean 
C.V. 


5.53955 

0.31045 


Adj  R-sq 


0.9988 


Parameter  Estimates 
Parameter  Standard  T  for  HO: 

Variable  DF  Estimate  Error  Parameter  =  0  Prob>|T| 

INTERCEP  1  -4.391624  0.03172280  -138.437  0.0001 

LNRE  1  0.860922  0.00257007  334.980  0.0001 

LNPR  1  0.442776  0.00362428  122.169  0.0001 
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PRANDTL  NUMBER  REGRESSION  IN  APPROACH  2 


THREE  CALIBRATION  REFRIGERANTS 


The  SAS  System  1 

12:56  Sunday,  August  21, 1994 

Model:  MODEL  1 
Dependent  Variable:  PR 

Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

9 

40.84252 

4.53806 

969.099 

0.0001 

Error 

113 

0.52915 

0.00468 

C  Total 

122 

41.37167 

Root  MSE 

0.06843 

R-square 

0.9872 

Dep  Mean 

1.45033 

Adj  R-sq 

0.9862 

C.V. 

4.71827 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  |T| 

INTERCEP 

1 

96.908792 

38.77678932 

2.499 

0.0139 

X 

1 

-26.452258 

10.79544008 

-2.450 

0.0158 

Y 

1 

83.304684 

26.64656684 

3.126 

0.0023 

X2 

1 

2.459228 

0.99959736 

2.460 

0.0154 

XY 

1 

-16.020859 

4.93940149 

-3.243 

0.0016 

Y2 

1 

8.088965 

7.10238635 

1.139 

0.2572 

X3 

1 

-0.076756 

0.03078419 

-2.493 

0.0141 

X2Y 

1 

0.750949 

0.22852190 

3.286 

0.0014 

XY2 

1 

-0.668378 

0.65967019 

-1.013 

0.3131 

Y3 

1 

2.213377 

1.01672374 

2.177 

0.0316 
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FOUR  CALIBRATION  REFRIGERANTS 


The  SAS  System  1 

23:58  Saturday,  September  10, 1994 

Model:  MODEL  1 
Dependent  Variable:  PR 

Analysis  of  Variance 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value 

Prob>F 

Model 

9 

42.70508 

4.74501 

1012.626 

0.0001 

Error 

159 

0.74505 

0.00469 

C  Total 

168 

43.45013 

Root  MSE 

0.06845 

R-square 

0.9829 

Dep  Mean 

1.50972 

Adj  R-sq 

0.9819 

C.V. 

4.53416 

Parameter  Estimates 

Parameter 

Standard 

T  for  HO: 

Variable 

DF 

Estimate 

Error 

Parameter  =  0 

Prob  >  |T| 

INTERCEP 

1 

132.495075 

34.02680322 

3.894 

0.0001 

X 

1 

-36.446721 

9.44835075 

-3.857 

0.0002 

Y 

1 

102.201609 

23.81121972 

4.292 

0.0001 

X2 

1 

3.392931 

0.87286660 

3.887 

0.0001 

XY 

1 

-19.449363 

4.40931003 

-4.411 

0.0001 

Y2 

1 

14.165714 

6.12753638 

2.312 

0.0221 

X3 

1 

-0.105753 

0.02682882 

-3.942 

0.0001 

X2Y 

1 

0.905506 

0.20379577 

4.443 

0.0001 

XY2 

1 

-1.253316 

0.57134480 

-2.194 

0.0297 

Y3 

1 

2.823725 

0.81305846 

3.473 

0.0007 
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APPENDIX  F.  DATA  ANALYSIS  PROGRAM 


This  appendix  contains  a  copy  of  the  FORTRAN  program  used  to  reduce  the  raw  data  from  the  transport 
property  tests.  This  program  opened  the  raw  data  files  which  were  collected  from  the  data  acquisition  system  used 
to  collect  all  of  the  property  and  heat  transfer  information.  Two  approaches,  the  Nusselt  number  method  and 
Prandtl  number  method,  were  also  used  to  analyze  the  thermal  conductivity  in  this  program.  All  the  equations 
referred  to  in  this  program  were  mentioned  and  discussed  in  the  previous  chapters. 

Properties  needed  were  obtained  from  the  ASHRAE  Handbook  or  REFPROP  (version  4.0),  which  were 
mentioned  in  the  Appendix  A,  and  written  in  subroutine  form.  Calculated  or  measured  properties  (including 
thermal  conductivity,  viscosity,  specific  heat,  and  density,  etc.)  are  available  to  execute  this  computer  program. 
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program  analysis 
implicit  real*8(a-h,o-z) 

oommon/acp/acpOtacpl,acp2,acp3,acp4,acp5 

common/av/av0,avlfav2,av3  ,av4,av5 

common/ak/akO,ak  1  ,ak2,ak3  ,ak4,ak5 

common/gpda/ga,pa,da 

cornmon/gpdb/gb,pb,db 

oommon/gpdc/gc,pc,dc 

common/cc/cc^ccg,ccp,ccd 

dimension  tw(15yttw(13Xtwo(13),twi(13) 
real*8  kUdMiu 

data  di,dl,pai/0.009525,2.0,3.141592654/ 
c  toanv»(ti+81.485443iyi.00110283e+04  !to  an^current) 
c  toamp=(ti+8 1. 485443  iyi. 001 10283e+04 

c . open  output  files 

open(unit=9/ile=,pr.dat',sUtiis=lunknownf) 

open(unit=l  O^^lognu.dat'^tatus^nknown’) 

open(unit=l  l,file=Vnu.datVstalus=Wcnownf) 
open(unit=  1 2411  e=,cfiihVstatus=\inkno\vnT) 
open(unit=  1 3/Ile^yynu.dat'rstAtus=WaK)wn,) 
open(iuiit=  M^le^Vawl  .dat\sUtus=‘unknownT) 
open(unit^l54ile=,raw2.dat,^Utus=,unknown') 
oper^unit3  1 6,file*^refincd.dat^status=Hinknownr) 
ope^um^^jfile^sf.datVstatus^Hinknown') 
open(unit=  1 84»le=loss.dat,^Utus-lunknownr) 
open(unit=  1 94He-cp.dat,rstatus=,unknovvnr) 
open(unit=2  l4ile^VcT.dat,^Jdais=HinkrK>vvnT) 
open(llnit=8/Ile=,ypgnu-dat,^tatus=Wcnown,) 
open(unit=7Tfile=‘modifiedk.dat'4tatus^lunknownr) 
open(imrt=23^Ue=,new.k.dat,,status=,unkiK)vvnr) 
open(umt=2241e=WLdat,^tatus^Wknown’) 
open(unrt=24/ile=Viscp.dat,^tatus=,unkiK)wn*) 

c  —  - 

c . Cp,  viscosity,  K  curve  fitting  coefficient  from  ASHRAE  or  REFPROP 

c . nr=l  for  r22;  nr=2  for  rl2;  nr=3  for  rl  13;  m=4  for  rl  14; 

c . nr=5  for  r236ea(REFPROP); 

c 

write(*,*),input  refrigerant  number,  nr=7, 
&’nr=El(r22)^ir=2(rl2)mr=3(rl  13);nr=4(rl  14);’, 
&'nrc5(r236ea)*4ir=s6(90%  E.G.);nr=7(80%  E.G.);', 

&'nr=8(70%  E.G.)^irss9(60%  E.G.)*,nr=  10(50%  E.G.)’ 
read(*,*)nr0 

write(V)'input  ending  refrigerant  number,  nend=7 

read(*,*)nend 

nr=nrO 

write(*,*)'inpiit  number  of  data  set,  NTDATA=7 

read(*,*)ntdata 

ndata=l 

90  if(nr.eq.l)then 
c . r22  transport  properties  data 

c . r22  Cp,  viscosity,  K  curve  fitting  coefficient  from  ASHRAE 

1  tc=96.15+273.33 
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tb— 40.86+273.33 
pc-733.02 
wm»86.47 
w-0.2211 

call  prop(nr)  lcall  prop(nr)  to  properties  coefficients 

open(uiift*531e«W2.dat,^Utus=,ol<f) 

goto  30 

c 

c . rl2  transport  properties  data 

c . rl2  Cp,  viscoaity,  K  curve  fitting  coefficient  from  ASHRAE 

c 

else  iflnr.eq.2)then 

2  tc=l  11.8+273.33 
tb-29.76+273.33 
pc-606.26 
wm=120.91 
w-0.1814 

call  prop(nr)  !call  prop(nr)  to  properties  coefficients 
open(unit=5rfile=,rl2.dat'rstatus='ol£f) 
goto  50 
c 

c....t1  13  transport  properties  data 

c . rl  13  Cp,  viscosity,  K  curve  fitting  coefficient  from  ASHRAE 

c 

else  if{nr.eq.3)then 

3  tc=2 14.35+273.33 
tb=47. 65+273 .33 
wm=187.38 
w=0.2563 

call  prop(nr)  {call  prop(nr)  to  properties  coefficients 
open(unit=5tfile=,rl  B.dat’^tatus- ol<f) 

go  to  50 
c 

c . rl  14  transport  properties  data 

c . rl  14  Cp,  viscosity,  K  curve  fitting  coefficient  from  ASHRAE 

c 

else  iffar.eq.4)then 

4  tc=145.65+273.33 
tb=3.65+273.33 
pc=471.08 

wm=  170.92 
w-0.2511 

call  prop(nr)  {call  prop(nr)  to  properties  coefficients 
open(unit=5^1e-,rl  1 4.dat\status-  ohf) 
goto  50 
c 

c . r236ea  transport  properties  data 

c . r236ea  Cp,  viscosity,  K  curve  fitting  coefficient  from  REFPROP-4.0 

c 

else  if(nr.eq.5)then 

5  tc= 139.29+273.33 
tb=6.5+273.33 
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pc=5 12.42 
wm=l  52.05 
w*0.3860 

call  prop(nr)  !call  prop(nr)  to  properties  coefficients 

opcn(unit=5tfile=,r236ca.dat,^tatus=,ol<r) 
c  open(unit*5 ,file^/home/yuan/viscometer/r236^  dat\ 

c  Aststus-olcT) 
goto  50 
end  if 


c . output  bead  lines 

50  write(V)Tave(C)  mass(kg/s)  heal(W)  h(W/mA2.C)  Nu' 
wr T(C)  Kg(W/m.C)  Kp(W/m.C)  Kd(W/m.C)  Kt(W/m.C)  E(% 
ScyG  E(%)-P  E(%>D* 

write^,*)’ - — - - - 

& - ’ 

write(V) 

areasspai*di*dl 

nexp=20 

do  1000  i^nditijiitdata 
c 

c . reading  data  over  nexp(=20)  times 

c 

do  200  k=l,nexp 
reid(V,end=1001Xtw(j)j=l,15) 

read(5Acnd=1001)power,ti,to,pd,pifden,tsefMmass 

c . ltd  anti-calibration  for  checking  accuracy  purpose 

c  tiamp=(ti+81. 223861 7)/ 1.001 84971 e+04 

c  ti=-80.0+(l  .0e4O4)*tiamp 
c  toamp=(ti+8 1. 485443 1 yi. 00 110283e+O4 

c  to==*80.0+(l  .0e+04)*toamp 
tf=0.5*(ti+to) 
c 

c . fluid  transport  properties  at  Tsat:  Cp,  vis,  and  K 

c . (ASHRAE  data  from  curve  fitting) 

c . coefficients  obtained  from  calling  prop(nr)  previously 

c 

fcp=acp0+acp  1  *tf+acp2  *(tf**2)+acp3*(tfM3)+acp4*(tf>M)+ 
&acp5*(tf**5)  !Cp 

fvis=av0+avl*tf+av2*(tf**2)+av3*(tf**3)+av4*(tf**4)+ 

&av5*(tf**5)  !  viscosity 

flct=ak0+akl  V+ak2*(tf**2)+ak3*(tf**3)+ak4*(tf**4)+ 

&ak5*(tf**5)  {thermal  conductivity 

c . net  heat  input  to  the  test  section 

heat*roiass*fcp*(to-ti)  !real  heat  input  to  the  test  section 

hflux=heat/area  !heat  flux  of  the  test  section 

c 

c . inner  wall  temperature  correction  by  heat  conduction  through  wall 

c 

rl=(3.0/16.0)*0.0254  !tube  inside  diameter 

r2=(4.0/16.0)*0.0254  !tube  outside  diameter 

twc=hflux*rl  *log(r2/riy40 1 .0  {wall  temp.  correction(wall  conduction) 
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do  300 11=1,13 

twi(n)=tw(n>twc  '.inner  wall  temp. 

300  continue 

c  twi(10)-0-5*(twi{9>+twi(n))  !b*d tw(10),  so  use  9,1 1  for  average 
call  dtwf(twi,ti,to,dt)  {call  dtwf  to  get  dt 

h-heal/(area*dt)  Ibeat  transfer  coefficient 

c 

c... -total  quantities  for  wall  temp.,  fluid,  power,  insulation  surface, 
c _ and  air  temperature 


ht~bt+h 

tti-tti+ti 


tto=tto+to 
do  180  n=l,13 

ttw(n>=ttw(n>Hw(n)  leach  outer  wall  temp,  total  quantities 

180  continue 

ttwa-ttwa+twa  !twa-average  wall  temp. 


c 


c 


ttfa=ttfa+tf 

tts=tts+tw(14) 

tta=tu+tw(15) 

tt*tt+t 


!tfa=average  fluid  temp. 

!tw(  1 4)-insul  ati  on  surfeace  temp. 
!tw(15)*=ambient  temp. 
!t«viscometer  temperature(C) 


tpi*=tpi+pi  !pi«inlet  pressure(psta) 

tpower^tpower+power  !  power  •power  input(l*V) 

tvismm=tvismm+vismm  !vismm-=viscosity  measured  by  viscometer(cp) 


trmass=trmass+nnass  !imass=mass  flow  rate 
tden=tden+den  !den*=density 

cpt~cpt+cp  !cp*specific  heat  capacity 

vist=vist+vis  !vis=viscosity 

ktt^^ktl+kt  !kt=thermal  conductivity 


200  continue 


c 

c . average  quantities  for  nexp(=20)  times 

c 


do  190  n=l,13 
two(n)=ttw(nynexp 
twi(n)=ctwo(n)-twc 
190  continue 
tiave*tti/nexp 
toave*=tto/nexp 
twa~ttwa/nexp 
tfa-ttfa/nexp 
c  t-tl/nexp 
pi*tpi/nexp 
tsavr-tta/nexp 
taave=tta/nexp 
powcra*5tpower/nexp 
c  vismm3StvUmm/nexp 
nnass-trmass/nexp 
den*tden/nexp 
have-ht/nexp 
c  epave-ept/nexp 
c  visave-vist/nexp 


leach  average  outer  wall  temp,  of  nexp  times 


♦average  inlet  temp. 

{average  outlet  temp. 

'average  wall  temp. 

{average  fluid  temp. 

{average  sample  cell  tempcrature( viscometer) 
{average  inlet  pressure 
{average  insulation  surface  temp. 

{average  ambient  temp. 

'.average  power  input 

{average  viscosity(cp)  measured 

{average  mass  flow  rate 

{average  density 

{average  heal  transfer  coefficient 

{average  cp 

{average  visoostty 
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{average  thermal  conductivity 


c  ktave=ktt/nexp 
c 

c.....P*at 

if(nr.eq.l)then 
w-0.2211 
pc“733.02 

paat=((5.00004649E-01)+(1.61909014E-02)*tfa+ 
ft  (1.95017230E-04)*tfa**2+ 

&  ( 1.03241 905 E-06)*tfa**3 +{2.23205 832E-09)*tfa**4+ 

ft  (1.71046532E-12)*tfa**5)*1.0E+06*1.45038E-04 
else  if(nr.eq.2)then 
w=0.1814 
pc=606.26 

p«at=((2.91020423E-01>+<100844698E-02)nfa+ 
ft  (1.49919608E-04)*tfa**2+ 

ft  (5.64050538E-07)*tfa**3+(-5.91876770E-09)*tfa*M+ 
ft  (5.79094585E-ll)*tfa**5)*1.0E+06*1.45038E-04 

else  if(nr.eq.3)then 

w-0.2563 

pc*501.25 

paat~((1.58554465E-02)+(6.77690492E*O4)*tfa+ 

ft  ( 1 .2 5 84830 1 E-05)*tfa**2+ 

ft  (2.28212969E-07)*tfa**3+(-1.47080736E-10)*tfa**4+ 
ft  ( 1 .60667465E-1 2)*tfa**5)*  1 .0E+06* 1.45038E-04 

else  if(nr.eq.4)then 
w*=0.25 1 1 
pc=471.08 

psats=((8.76785517E-02>+(3.47902835E-03)*tfa+ 
ft  (5.45874973E-05)*tfa*#2+ 

ft  (3.74883 1 15E-07)*tfa**3+(4.43446752E-10)*tfa**4+ 

ft  (1.8815001 6E-1 2)*tfa**5)*  1 .OE+06* 1 .4503 8E-04 

else  il(nr.eq.5)then 

w=0.3860 

pc-512.42 

p«at=(1.13349352E+01)+(4.84870791E-01)*tfa+ 
ft  (8.3043659  lE-03)*tfa**2+ 

ft  (6.785 1011 4E-05)*tfa**3+{2.04624 143E-07)*tfa**4+ 
ft  (-2.19035401E-10)*tfa**5 
end  if 

c....iiuid  transport  properties  based  on  temperature  at  tfa:Cpf  vis,  and  K 

c . (ASHRAE  data  from  curve  fitting) 

c . coefficients  obtained  from  calling  prop(nr)  previously 

c 

cp~acp0+acp  1  *tfa+acp2*(tfa**2)+acp3  *(tfa**3)+acp4*(tfa**4)+ 
&acp5*(tfa«5)  !Cp 

vis-av0+avl*tfa+av2*(tfa**2>+av3<(tfa**3)+av4*(tfa**4)+ 

&av5*(tfa**5)  {viscosity 

kl=ak0+akl*tfa+ak2*(tfa**2)+ak3*(tfa**3)+ak4*(tfa**4)+ 

&ak5*(tfa**5)  {thermal  conductivity 

c . non-dimensionized  temperature  and  pressure 


tr~(tfa+273.33ytc 


!tc  is  critical  Temp. 


ppr-pVpc 

dpr-(jK-p«*typc 

call  kfceto»<lr4ipr^»liok) 
call  vfcctoc(w,tr,dpf^atiov) 
wriu(*t*yraliok“,j»liok,^aliov='  .raliov 


!pc  is  critical  pressure 

Iptal  it  saturated  pressure 

{k  midified  factor  based  on  current  pressure  pi 

{viscosity  modified  factor  based  on  pi 


vis~vii*raliov  {modifying  viscosity  from  saturated  stale  to  pi  stale 
viscp-vis*cp 

kt~kt*ratiok  {modifying  k  from  saturated  state  to  pi  stale 
call  dtwf(twi,tUveftoavetdt)  {call  dtwf  to  get  dt 

dtn*<toave4iaveydt  {non-dimensionized  dt 

twave-tfa+dt 

h*nnass*cp*(toave^tiavey(area*dt)  {beat  transfer  coefficient 

c _ write  aveaged  data 

write(143000)ittfa4mass,tiave,toavettwave,dt,taavettsave,powera 
3000  fonnat(lx»D 


c 

c . heat  loss  calculations  and  calibrations 


tak-Uave+273.33 

call  aiiprop(tak,dena,cpa,vUa,tka)  {call  airprop  to  get  prop. 


c . heat  loss  calculations 

c 


!air  Prantdl  number 
lair  expansion  coefficient 
{gravitional  acclcration 

lair  kinennatic  viscosity 
!air  thermal  diffiisivity 
{outer  insulation  diameter 
{outer  insulation  surface  area 
!temp.  difference  between  surf/air 
{temp,  difference  between  wall/air 
{temp.  diff.  ratio  between  suf-wall-atr 
!Ri  number 


pra-visa*cpaAka 
beta=1.0/taave 
g=9.8 

dvisa^visa/dena 
alpha*stka/(dena*cpa) 
dia=2.5 *2.54/100.0 
ao*pai*dia*dl 
dua^tsave-taave 
dtwa  ^twave-Uave 
rdtswa^dtsa/dtwa 
rmda»=(g*beU*dtsa*dia*  *3 .0y(dvisa*alpha) 


c ch urchill  free  convection's  Nusseh  number  calculations 


beat“rmass*cp*(toave-tiave)  {real  heat  input  to  the  test  section 

hflux~heat/area  {heat  flux  of  the  test  section 

pcir=(l  .0+((0.5  59/pra)**(9.0/l  6.0)))*  *(8.0/27.0) 
anuic-<0.6+(0.3r7*(rada**(1.0/6.0)ypm))**2.0  lair  Nu(unmodify) 


yanuic-log(anuic) 

anu~1.2*anuic 

ha~anu*tka/dia 

qloss~ha*ao*(Uave-taavc) 

qlossx-powera-heal 

hac-qlostx/ao/(tsave-Uave) 

anuc«hac*dia/tka 

yanuc-log(anuc) 

cfrqlostx/qloas 

ybcf^logCcf) 


!log(anufc) 

{modified  Nu  number,  update  1.2 

{air  heat  loss  cofficient,  ha 

{heat  loss  calculated  from  C.C.  eq. 

{true  heat  loss 

{true  ha 

{true  Nuair 

!log(anuc) 

{heal  loss  ratio 
Hog  of  heal  loss  ratio 
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xbn>log(rada)  HogofR* 

write(18^250)log(dU»Xlog(qlossx)  ! write  to  loss.dit 

2250  formal  x3(n5-6^x)) 
c 

beat  input 

c....  .heat  lots  prediction  by  measured  Ts-Ta 

qlo«s*(exp(0.612001419))*(dtsa)**1.33022535  !beat  loss  prediction 

cpiet~powera-qlo«8  !net  heat  to  the  test  section 

qbatio*(qloas/powera)*100.0  Ibeat  loss  percentage  ratio 

cpc°=qnet/(rmass*(toave4iave))  !  calculated  Cp 

rcp=((cpc<pycp)*  100.0  !Cp  deviation  percentage 

write(19,2260)tfa,cp,cpcjrcp  twrite  to  cp.dat 

2260  format(lx^l0.5,2x£fl2.6,2x47.3) 
c.....Re  and  Pr  number 

rt=(4.0*nnassy(pai*di  *  vis)  !  Reynolds  numbers 

xre=log(re)  !  logarithm  Re  transformation 

pr=vis*cp/kt  Ifluid  Prantdl  number 

ec=rmass**2/(den**2*di**2*cpc*dt) 

write(24,10  !3)tfaJcMniass,vis,cpc 
write<23,*)log(pr),log(re),log(dtn) 
c  write(23,1013)kl^inass,vis#1.0E+03,cp/1000.0,dtn 
c  write(*,*)lvis=^vis,,cp=,,cpfrinass-,rmass 
call  newk(prtre,dtn) 
if(nr.eq.5)then 
hk=vis*cpc/pr 
else 

hk=vis*cp/pr 
end  if  . 
hk=hk/ratiok 

rhkK(bk-ktykt)*100.0 

write(22,1014)tfa4}kjct4frk 
1014  format(lx49-6,2(fl2.6,lx)^10.6) 
c 

c . parameter  logarithm  transfonnation  for  corrected  Nusselt  number 

c 

pr=vis*cp/kt  {fluid  Prantdl  number 

xpr=log(pr)  !log  Prandtl 

nu=h*di/kt  {corrected  Nu  number 

ynu=log(nu)  SlogNu 

yynu=log(nu/(pr**0.507056)) 
write( 1 5,300 1 )i,tfa,tr,re,pr,cp 

3001  format(lx,i3,lx,2(f6.3,lx),lx^l0.2,lxJ8.4,lxfL0.5) 
c 

c . correlation  Nusselt  number  calibration  functions 

c . call  subroutine  cnucf(re,gc£pcf,dcf) 

c 

call  cfiiu(pr,rc,enu,cfg,cip,cfd,rfh) 
c 
c 

c . modified  h  values  from  enu(curve  fitting  equation)  values 

c 

hstar=(enu*ktydi 
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cfh=hstar/h 

ybfUr*log(h*Ur) 

yh»log(h) 

hl=h*cfh 

yyb=log(h/pr*  *(-0.477820)) 
yyh*Ur*log(hiUr/pr#*(-0.477820)) 
c 

c . Dittus  Boeher  correlation 

c 

c  h=have 
*1=0.905027 
b 1=0.507640 
cl  =exp(-5.0 1 5673) 

*2=0.800013840 

b2=0.4 

c2=exp(-3.772261) 

eml=al-*2+0.8 

em2=bl-b2+0.4 

en=1.0/(0.6-bl+b2) 

dks=((c2*h*diy(0.023*cl*re**eml*(vis*cp)**e!n2))**en 

dks=dks/ratiok 

dnust*r=h*di/kt 

dnu=0.023*re**0.8*pr**0.4  ’Dittus-Boelter  eq.  Nu 
ydnu=log(dnu)  !ln(Nud) 

ydcf=ynu-ydnu 
dcf=nu/dnu 

dk=((h*diy(0.023*ie**0.8*(vis*cp)**0.4))**(5.0/3.0) 

dk=dk/ratiok 

1980fomut(lx49.5,lx49.5,lx,4(n2.6flx)) 

prd=vis*cp/dk 

c  prstar=7.4 13721 08e-2+9. 8 5 503972e-l  *prd 

c  dk=vis*cp/prrt*r 

nik={(dk5-ktykt)*  100.0  !  deviation  ratio 

if(abs(rdk).le.5.0)then 

write(16,1800)powera4mass,den, 

&tfa,tiavc,toavc,twavc,dt,taave,tsave 

else  if(*bs(rdk).gt5.0)then 

goto  800 

end  if 

1800  form*t(lx49.44x48.54x^l0.3,lx,7(f8.5,lx)) 

800  write(*,*)re,pr,xre,xpr4iu,dnu 
write(  10,201 0)re,pr,xre,xp!\nu,dnu 
2010  form*t(lx/15.6,lx48.4,lx,4(fl2.6,2x)) 
whe(*,900)tr,tfa,pr3ctTdk4tlk 

write(9,900)tr,tf*,prjct,dk4tik,dens,den  ’write  to  pr.dat 
900form*t(lx^3.2,2x^5.2,2x48.4,2x^(f8.5,2xXf8.4,2x,2(f83,2x)) 

c . Petukhov  Popov  correlation 

f=(1.82*logl0(re)-1.64)**-2 

pnul={f78)*re*pr 

pnu2=1.07+12.7*((f78)**(0.5))*(pr**(2./3.)-1.0) 

pnu=cpnul/pnu2 

ypnu=log(pnu) 
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ypcf-ynu-ypou 

pcf*nu/pou 

cstar-=cl 

astar*al 

b*Ur*bl 

cpe«exp(-4.052869) 

ape=0.826843 

bpe*0.457499 

prpold=pr 

850  ipl^crtar*(fr8.)**0.5*vis*q>>re**(astar-apc+l))/ 
&<12.7*cpe*h*di) 

4>2=(1.07/(12.7^8.)^*0.5)) 

^>==prpold**(273.)-fpl*pfpold**(bstar*bpcy»-^)2-l 

&prpold**(b*Ur-bpc-l) 
pqmew=prpotd^lp/ipp) 
epcp=prpnew-prpold 
if(aba(epsp).le.l.0e-12)goto  920 
prpold^prpocw 

goto  850 

920  pk=vis*cp/pfpncw 
pk=pk/ratiok 

ipk=«pk-iciyki)*ioo.o 

c . Gnietinski  correlation 

fg=(0.79*log(re>1.64)"-2 

gnul^fg/8)*(re-1000.0)*pr 

gnu2=1.0+12.7*((fg/8)**0.5)*(pr**(2./3,)-1.0) 

gnu=gnul/gnu2  IGnielinski  Nu 

gkc=vis*cp*((gnu  l/(  I2.7*(fg/8)**(0.5)*h*di))- 

&(l/(12.7*(fg/8)**0.5))+l)**(-1.5) 

ygnui=log(gnu)  !ln(Nug) 

ygcf=ynu-ygnu 

gcf=nu/gnu 

cge=exp(-4.391624) 

age=0.860922 

bge=0.442776 

prgold=pr 

860  fgl=(c8Ur*(fg/8.)**0*5*vis*cp*(re-1000.)**(astar-age+l))/ 
&(12.7*cge*h*di) 
fg2=(1.0/(12.7*(f^8.)**0.5)) 
ffg==Prgold**(2^3  .)-fgl  *prgold*  *(bstar-bge)+fg2-l 

73  .)-(bstar-bge)*fgl  *prgold** 

&(bsUr-bge*l) 
prgr»ew=prgold-{£fg/ffgp) 
q*g=prgnew-prgold 
if(abi(epag).le.l.0e-12)go  to  930 
prgold=prgnew 
goto  860 

930  gk=vi**cp/prgnew 
gk^gk/ratiol k 
rgkK(g^yfc)M00.0 
c . write  to  output 


c  gkcc^(gko-ktykt)*100.0  !k  deviation  from  Gnielinski  cq. 
c  if(gkcc.gtlO.O)go  to  990  {based  ongkcto  throw  away  bad  points 
write(8,2500)i^,xpr,ynu,ydnu,ypnu,ygmi  ‘write  to  file  ypgnu.dat 
990  write(7^510)tf3UcUdkstpk,gk^dk^kti^k  ‘write  to  file  modified.k.dat 
write(2  lr2520)xre,xpr,ydcf,ypcf,ygcf  Swritetofileycf.dat 
2500  formal  lx,i4,lx,6{fL  5.6,  lx)) 

2510  format(lx415.3t2x,4(f8.6,2xX2x3(n.3^x)) 

2520  format(  1  xv2x,5(fl  2.6,2x)) 
c . Nusselt  number  curve  fittings 

c . n  for  Pr^n  in  3  and  4  calibration  fluids(R22,R12,Rl  13,  add  R1 14)  are: 

c . expriment:0.498730;0.504359 

c.....Dittus:0.4;0.4 

c . Petukhov:0.455420;0.457575 

c . Gnielinki:0.441340;0.442798 

yenu=log(nu/pr**0. 5230 86) 
ydnu=log(dnu/pr**0.4) 
ypnu=log(pnu/pr**0.457575) 
ygnua=log(gnu/pr**0.442798) 

write(  1  l,2600)xre,yenu,ydnu,ypnu,ygnu  !  write  to  ynu.dat 

2600  format(2x,5(fl2.8,2x)) 
c 

c..  .read  another  set  of  data  at  different  temp,  level 
c . set  total  quantity  equal  to  0.0 


600  do  601  n=l,13 
ttw(n)*0.0 

601  continue 
ht=0.0 
tti=0.0 
tto=0.0 
ttwa=*0.0 
ttfa=0.0 
tts=0.0 
tta=0.0 

c  tt=0.0 
tpi=0.0 

c  tvismm^O.O 
tpower=0.0 
trmass=0.0 
tden=O.0 
ndata=i+l 


! total  quantity  of  each  wall  temp. 

!  total  quantity  of  heat  transfer  coef. 

‘total  quantity  of  inlet  temp. 

‘total  quantity  of  outlet  temp, 
total  quantity  of  average  wall  temp, 
total  quantity  of  average  fluid  temp. 

!  total  quantity  of  insulation  surface  temp, 
{total  quantity  of  ambient  air  temp. 

{total  quantity  of  sample  cell  temperature 

{total  quantity  of  inlet  pressure 

{total  quantity  of  viscosity 

{total  quantity  of  power  input 

{total  quantity  of  mass  flow  rate 

{total  quantity  of  density 

{number  of  data  set(each  set  has  nexp  times) 


1000  continue 

c . read  another  working  fluid  data  set 

1001  if(nr.eq.nend)go  to  2000 
close(unit-5) 

ndala=i 


goto  90 

1010  formal(  1  x,i2,2x,6(fl  5.8)) 

1011  format(lx46.3,7(fl0.5,lx)) 

1012  format(  1  x»i2,2x,6(fl  5.8)) 

1013  format(lx,5(fl5.8,lx)) 


‘end  of  data  file  detected  come  here  and  check  nr 
{temporary  close  5  and  ready  open  5  for  another  fluid 
{set  data  number  continuously 
{number  of  refrigerant  increased  by  1 
{back  to  very  beginning  of  data  reading 
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2000  clo«e(iinit=5) 
stop 
end 

£••*••*•*••*•**••*•***•*****•*♦*•**•*•*****•**•*•*♦*•**************••***♦**** 
C  This  subroutine  calculates  the  refrigerant  properties,  such  as  specific  * 

C  heat  capacity,  viscosity,  and  thermal  conductivity,  by  curve  fittings.  * 

C  The  source  data  are  based  upon  ASHRAE  except  R236ea  based  on  REFPROP4.0  * 

subroutine  prop(nr)  !  input  number  of  refrigerant 

implicit  real*8(a-h,o-z) 

cooimoa/acp/acp0,acpl,acp2,acp3,acp4,acp5  !Cp  coefficients 
common/av/avO,avl,av2,av3,av4,av5  !  Viscosity  coef. 

convnorv/ak/akO,ak  1  ,*k2,ak3tak4,alc5  !K  coef. 
if(nr.eq.l)then 

open(unit=20,file==,prop.22Vstatus-oldr)  !coef.  file  for  HCFC-22 

c . coefficients  for  HCFC-22  Cp 

read(20,*)acp0,acpl,acp2,acp3,acp4,acp5 

c . coefficients  for  HCFC-22  viscosity 

read(20,*)av0,avl,av2,av3,av4,av5 

c . coefficients  for  HCFC-22  K 

rcad(20,*)ak0,akl,ak2tak3fak4,ak5 

clo«e(unit=20) 

return 

else  if(nr.eq.2)then 

open(unit=204ile=*prop.  12',status^'olcT)  !coef.  file  for  CFC-12 

c . coefficients  for  CFC-12  Cp 

read(20,*)acp0,acpl,acp2,acp3,acp4,acp5 
c..... coefficients  for  CFC-12  viscosity 
read(20,*)av0,avl,av2,av3,av4,av5 

c . coefficients  for  CFC-12  K 

read(20,*)ak0,ak  1  ,ak2,ak3,ak4,ak5 

close(umt=20) 

return 

else  if(nr.eq.3)then 

open(unit=204ile=,prop.l  13',status=’ol<f)  !coef.  file  for  CFC-1 13 

c . coefficients  for  CFC-1 13  Cp 

read(20,*)acp0,acp  I,acp2,acp3,acp4,acp5 

c coefficients  for  CFC-1 13  viscosity 

read(20,*)av0,av  1  ,av2,av3,av4,av5 
c..., .coefficients  for  CFC-1 13  K 

read(20,*)ak0,alcl^k2,ak3,ak4,ak5 

close(unit-20) 

return 

else  if(nr.eq.4)then 

open(unit=20,file=,prop.  1 14’^tatus=’oIcT)  !coef.  file  for  CFC-1 14 
c.....coefficients  for  CFC-1 14  Cp 

read(20,*)acp0,acpl,acp2,acp3,acp4,acp5 

c . coefficients  for  CFC-1 14  viscosity 

read(20,*)av0,avl,av2,av3,av4,av5 
c..... coefficients  for  CFC-1 14  K 

read(20,*)ak0,ak  1  ,ak2,ak3  ,ak4,ak5 
close(unit=20) 
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return 

else  if{nr.eq.5)then 

opcn(unit*20ffile*,prop.236,^Utus-olcr)  looef.  file  for  HFC-236ea 

c . coefficients  for  HFC-236ea  Cp 

read(20,*)acp0,acpl,acp2,acp3,acp4,acp5 
acp0=  1.1 93407 10E+03 
acpl=  6.33325 863 E+00 
acp2=-4.13577229E-01 
acp3=  1.47026302E-02 
acp4=-2.220295  86E-04 
acp5=  1.21 074845 E-06 
c  acp0=1.20477688E+03 
c  acpl=1.92495738E4)0 
c  acp2*=0.0 
c  acp3=0.0 
c  acp4*0.0 
c  acp5=0.0 

c . coefficients  for  HFC-236ea  viscosity 

read(20,*)av0,avl,av2,av3,av4,av5 
av0=(  5.63759208E-01)*1.0e-3 
avl  =(- 1 .03053898E-02)*  1 .0e-3 
av2=(  2. 1 434 1 562E-04)*  1 .2e-3 
av3={-5. 0634444 1 E-06)*  1 .0e-3 
av4=(  6.690402 19E-08)*1.0e-3 
av5*=(-3.19661297E-10)*1.0e-3 
^....coefficients  for  HFC-236ea  K 
read(20,*)ak0,ak  1  ,ak2,ak3  ,ak4,ak5 
close(unit=20) 
return 
endif 
end 

C  This  subroutine  calculates  the  air  properties,  such  as  density,  Cp,  * 

C  viscosity,  and  thermal  conductivity,  by  curve  fittings.  The  source  data  • 

C  are  based  on  back  tables  in  Appendix  of  the  book  by  Incoporea  and  DeWitt* 
£*•«**********»*»****»*****••*******•***» ***************************** 
subroutine  airprop(tak,dena,cpa,visa,tka) 
implicit  real*8(a-h,o-2) 

data  d0,d  1  ,d2,d3,d4,d5/8.402 1 902 1 e+00, -7.3 1001 496e-02, 
&3.01471766e-04,-6.34908758e-07,6.58567756e-10,-2.66433143e-13/ 
data  c0,cl,c2,c3,c4,c5/l.l  15073 12e+03 ,-1.3 15441 13e+00, 
&5.91507275e-03,-1.29040736e-05, 1.44722 101e-08,-6.33937087e-12/ 
data  v0,vl  ,v2,v3,v4,v5/-3 .79 5 14205e-07,+8.36939549e-08, 
&-9.95046476e-ll,1.33037244e-13,*1.27004888e-16,5.52250075e«20/ 
daU  t0,tl,t2,t3,t4,t5/-4. 1 9377378e-04,l  .038693 52e-04, 

&-7.4703 5287e-08, 1 .22 800797e- 1 0,- 1 .5 5 56836 1 e- 1 3 ,7. 1 02 1 5 3 5 5e- 1 7/ 

c . getting  density,  specific  heat,  viscosity,  and  thermal 

c..... conductivity  of  air  at  temperature  tak(in  K) 

dena=d0+d  1  *tak+d2  *tak**2+d3  *tak**3  +d4*tak**4+d5  *tak**  5 
cpa^cO+cl  •tak+c2*tak**2+c3*tak**3+c4*tak**4+c5*tak**5 
visa=v0+vl*tak+v2*tak**2+v3*tak**3+v4*tak**4+v5*tak**5 
tka=st0+tl*tak+t2*tak**2+t3*tak**3+t4*tak**4+t5*tak**5 
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return 

cad 


c  This  subroutine  calculates  the  average  wall  and  fluid  temperature  * 

C  difference  by  using  finite  control  volume  method  * 

subroutine  dtwf(twi,tiave,toave,dt) 
implicit  real*8(a-h,o-z) 
dimension  twi(13),tfi(13),dx(13),x(13) 
c 

c.... .finite  control  volume  grid  increament  and  coodinates 
c 

dx(l)=O.05 
x(l)=0.0 
dx(2)=0.1 
x(2)=0.1 
dx(3)=0.13 
x(3)=0.2 
do  10  i-4,10 
dx(i)=0.2 
x(i)=x(i-l>K).2 
10  continue 
dx(ll)=0.15 
x(ll)-0.18 
dx(12)=0.1 
x(12)=0.19 
dx(13)=0.05 
x(13)=2.0 
dl=2.0 
c 

c... .fluid  temperature  distribution(assume  linear) 
c 

tfi(l>=tiave . 
tfi(13)=toave 
do  20  i*=2,12 

t^i)^(i-i>+<x(i)-x<:i-i))*<(tfi(i3>tfi(i)ydi> 

20  continue 
c 

c . dt  calculations 

c 

dtt-0.0 
do  30  i-1,13 

dtt=dtt+<dx(iy(twi(i>tf!(i))) 

30  continue 
dt=dl/dtt 
return 
end 

C*«**«MM*iM««#*««M«M**«M«M»t*****«*«*W**«*«*«« 

C  This  subroutine  calculates  the  calibration  functions  for  the  Nusseh  * 

C  number1!  equations,  such  as  Gnielinski,  Petukhov,  and  Dittus-Boeher  * 
C  equations.  Also,  the  experimental  data  forms  the  Nusseh  number  * 

C  correlation,  Nue.  Those  are  obtained  by  curve  fittings  from  SAS.  * 
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subroutine  cfiiu(pr^e,enu,cf&cfp,cfd,cfii) 

implicit  real*8(a-h,o-z) 

common/gpda/ga,pa,da 

oommon/gpdb/gb,pb,db 

common/gpdc/gc,pc,dc 

common/cc/ccs,ccg,ccp,ccd 

data  ae,be,ce/0.85054475 1,0.507056, -4.39778757/ 

data  ag,bg,cg/0.835702658,0.495804,-4. 16347694/ 

dataap,bp,cp/0.819511116, 0.516636,-4.03165674/ 

daU  ad,bd,cd/0.799977124,0.400000,-3.77200365/ 

data  haJib,hc/0. 74666 1,-0.477820,-0. 159677/ 

data  hal,hbl,hcl/0.746638,-0.477820,-0. 1594 15/ 

ga=ae-ag 

pa-ae-ap 

da=ae-ad 

gb=be-bg 

pb=be-bp 

db=be-bd 

ccs=exp(ce) 

ccg-exp(cg) 

ccp=exp(cp) 

ccd=exp(cd) 

gc=ce-cg 

pc=ce-cp 

dc=ce-cd 

enu=pr**ae*re**be*exp(c€) 

rfg^P^<(^g))*(re**(1^ag))*(cM>(ce-cg)) 

cfp=^pr**(be-bp))*(re**(ae-ap))*(exp(ce-cp)) 

cfd=(pr**(be-bd))*(re**(ae-ad))*(exp(ce-cd)) 

cfti={pr<*(hb-hb  1  ))*(re**(ha-ha  1  ))*(cxp(hc-hcl )) 

return 

end 

c  The  following  subroutine  calculates  the  thermal  conductivity  * 
c  ratio  at  different  pressure  but  at  the  same  temperature.  This  is* 
c  so  called  pressure  effects  on  thermal  conductivity.  The  equation* 
c  used  here  is  Missenard  correlation 1970)  * 

subroutine  kfactor(tr,pr,ratiok) 
implicit  real*8(a-h,o-z) 
dimension  a(21) 

data  a(lXa(2Xa(3Xa(4Xa(5),a(6),a(7),a(8X 
&  a(9Xa(10),a(llXa(12),a(13),a(14),a(15),a(16), 

&  a(17Xa(18Xa(19Xa(20),a(21) 

ft  /-0.376480, 1 .998078,0.007849,-3.4 1 9245,-0.053279, 

ft  0.000131,1.956588,0.115588,-0.000281,-0.000001554, 

&  0.0,-0.075644,  -0.00029 1,0.000005608,-5.59027E- 10, 

ft  0.0,0.0,0.0005 17,-0.000004349,9.446296E-1 1,2.460594E-13/ 

q=a(l)+a(2)*ti^a(3)*piH-a(4)*ti^*2+a(5)*ti^pr+a(6)*pr**2+ 
ft  a(7)*tr**3+a(8)*tr**2*pr+a(9)*tr*pr**2+a(10)*pr**3+ 
ft  a(ll)*tr**4+a(12)*tr**3*pr+a(13)*tr**2*pr**2+a(14)*tr*pr**3+ 
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&  «(15)*pr**4+«(16)*tr**5+a(17)*tr**4*pr+»(18)*tr**3*pr**2+ 

&  »(19)*tr**2*pr**3+»(20)*tr*pr**4+*(21)*pr**S 

nliok=1.0+<j*pf**0.7 

return 

cod 

c  The  following  lubroutine  calculates  the  liquid  viscosity  * 
c  ratio  at  different  pressure  but  at  the  same  temperature.  This  is* 
c  to  called  pressure  effects  on  liquid  viscosity.  The  equation  * 
c  used  here  is  Lucas  correlation  1981)  * 

c************************************************************ 

subroutine  vfactor(w,tr,dpr,ratiov) 
implicit  real*8(a-h,o-z) 

a-0.9991-(4.674E-04/(1.0523*tr**(-0.03877>1.0513)) 

d^0.3257/(1.00394r**2.573)**0.2960>0.2086 

c=-0.07921+2.1616*tr-13.4040*tr**2+44.1706*tr**3- 

&  84.8291*tr**4+96.1209*tr**5-59.8127*tr**6+15.6719*ti**7 

ratiovKlO+d' ‘(dpr/2.1 18)**a)/(1.0+c^v*#r) 

return 

end 

c*************w********m**********************  **•* 
c  The  following  program  is  used  to  calculate  k  by  directly  curve  * 
c  fitting  by  function  log(h),  log(vis*cp),  and  log(re)  * 

^•••••♦•••••♦•••••♦••••••♦♦•♦♦••♦♦♦•••********************** 

subroutine  newfc(pr,re,dtn) 
implicit  real*8(a-h,o-z) 
data  a0,al,a2,a3,a4,a5,a6,a7,a8,a9 
&  /96.908792,-26.452258,83.304684,2.459228, 

&  -16.020859,8.088965,-0.076756,0.750949, 

&  -0.668378,2.213377/ 

x=log(re) 

y=log(dtn) 

pr*a0+al  *x+a2*y+a3  *x**2+a4*x*y+a5  *y*  *2+ 

&  a6*x**3+a7*x**2*y+a8*x*y**2+a9*y**3 

pr=exp(pr) 

return 

end 
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